m 


m 


A 003  ;  FINAL  REPORT 

COMPUTER  ALGORITHMS  AND  ARCHITECTURES 
FOR  THREE-DIMENSIONAL 
EDDY-CURRENT  NONDESTRUCTIVE  EVALUATION 


Contract  No.  N  00019-86-C-0219 


Sabbagh  Amoaatea,  Inc. 
4639  Morningaide  Drive 
Bloonangton,  IN  4T401 


Naval  Air  Syttema  Command 
20  January,  1989 


VolUBM  II 


CHAPTERS  I-V 


DT1C 


4  apr  nat 


r 


*  89  4  04  08 


REPORT  DOCUMENTATION  PAGE 

U  REPORT  SECUR.TV  CuASS  F<AT)ON 

UNCLASSIFIED 

•b  RESTRICTIVE  MARKINGS  1 

2a  SECuR-TY  C-ASSiFhCATJON  AoThORiTY 

3  DiSTR.Bur»ON'APA.LA8iUrr  Of  REPORT 

Approved  for  Public  Release; 

Distribution  is  unlimited. 

2b  OEC_ASS*f<CAT>ON  '  DOWNGRADING  SCHEDULE 

SA/TR-2/89  . 


6a.  NAME  Of  PERFORM.NG  ORGANIZATION 

SABBA6H  ASSOCIATES,  INC. 


6c  ADDRESS  (Oty.  Sun.  end  ZtP  Code) 

4639  Morningside  Drive 
Bloomington,  IN  47401 


Sa.  NAME  Of  FUNDING/ SPONSORING 
ORGANIZATION 

Naval  Air  Systems  Command 


Sc  ADDRESS  (Oty.  Stite,  end  ZIP  Coot) 

AIR-931A 

Washington,  DC  20361-9310 


'1  TlTwfi  (Include  Security  CliSSifidtion) 


60  OFFICE  SYMBOL 
(If  ififiliCtbl*) 


7*  NAME  OF  MONITORING  ORGANIZATION 

Department  of  The  Navy 
Naval  Air  Systems  Command 


7b  ADDRESS  (Oty.  Staff,  end  HP  Coo*) 

ATTN:  AIR-931A 
Washington,  DC  20361-9310 


80  OFFICE  SYMBOL  9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
*  (If  ipptkibie) 

N00019  N0001 9-86-C-0219 


io  SOURCE  of  funding  numbers 


PROGRAM 

PROJECT 

task 

WORK  UN*T 

ELEMENT  NO. 

NO 

NO 

ACCESSION  NO 

65502N 

COMPUTER  ALGORITHMS  AND  ARCHITECTURES  FOR  THREE-DIMENSIONAL 
EDDY-CURRENT  NONDESTRUCTIVE  .EVALUATION. 


■  2  PERSONAL  AUTmOR(S) 


13a.  TYPE  Of  REPORT 

Final 


16  supplementary  notation 


SABBAGH  ASSOCIATES 


13b  TIME  COVERED 

from  12/2/86  to  1/2Q/8 


|U  DATE  OF  REPORT  (YW.  Montn.  Oty) 

9\  1/20/89 


’7  COSATi  CODES  18  SUBJECT  TERMS  (Continue  on  revert*  <f  necessity  end  identify  by  block  number) 

group  sub-group  Eddy-current  nondestructive  evaluation.  Three-dimensional 

02 _ image  reconstruction,  inverse  problems,  conjugate  gradient 

04  algorithm,  computer  architecture  and  DSP  integraged  circuit 


'9  ABSTRACT  (Contrnvif  on  reverse  if  necessity  end  identify  by  block  number) 


SEE  NEXT  PAGE 


20  O  S*R  8.  *.ON  '  AVAILABILITY  O?  ABSTRACT 

d  UNCLASS  ‘  EDUN.  V:TE0  Q  $aME  A5  RPT  C  3:  C  USE'S 

2’.  ABSTRACT  SECjRiTY  CLASSIFICATION 

UNCLASSIFIED 

::a  SAVE  0s  RES'ONS-SlE  iND'ViD'jAl 

Dr.  L.  E.  Sloter 

22o  TE.Eri’ONt  (inerjot  Arte  Coot! 

(202)  692-7445 

dff.ce  s'  s' jo. 

AIR-931A 

OO  FORM  1473,  it  var 


83  APR  fc  '.  or  r-a,  „iea  j-t  e»-a»siea 
AH  oi^f '  fo  t  o-'s  a't  obso.fi* 


C-ASS  ='CAf.QN  O? 


19.  ABSTRACT 


z  2a  this  report^we  develop  an  electromagnetic  model  for  three-dimensional  inversion  of  eddy- 
current  data,  an  inversion  algorithm  based  on  the  conjugate  gradient  technique,  and  a  special 
purpose  computer  that  we  estimate  can  execute  this  algorithm  in  times  comparable  to  high 
speed  main-frames.  This  computer  has  a  pipeline  architecture  and  is  designed  around  our  parallel 
implementation  of  the  inversion  algorithm  and  makes  use  of  high-speed  DSP  chips.  The  inversion 
process  achieves  a  higher  performance  measure  when  more  than  one  data  set  is  inverted.  The 


sequential  order  of  the  inversion  scheme  restricts  the  number  of  active  elements  in  the  pipe  for 
a  single  problem.  When  more  than  one  inversion  problem  enters  the  pipe,  then  more  than  one 
element  could  be  active  to  improve  the  overall  performance  of  the  system. 


The  basic  electromagnetic  model  starts  with  the  integral  equations  for  electromagnetic  scat¬ 
tering,  which  are  then  discretized  by  means  of  the  method  of  moments.  This  gives  us  the  funda¬ 
mental  inversion  model,  which  is  then  solved  using  the  conjugate  gradient  algorithm.  In  order  to 
accomplish  the  three-dimensional  inversion,  we  acquire  data  at  a  number  of  frequencies;  there¬ 
fore,  our  inversion  process  is  called  a  multifrequency  method.  The  choice  of  frequencies,  and  the 
number  of  frequencies  to  be  used,  depend  upon  the  conductivity  of  the  host  material,  and  the 
depth  resolution  sought.  r  a  \  ^ _ 


The  method  of  conjugate  gradients  has  a  number  of  attractive  features  for  our  purposes.  Chief 
among  them  is  that  it  allows  a  large  problem  to  be  solved  efficiently,  and,  because  it  is  an  iterative 
algorithm,  it  allows  us  to  take  advantage  of  the  special  Toeplitz  structure  of  the  discretized  model. 
We  also  derive  an  algorithm  that  allows  us  to  constrain  the  solution,  use  preconditioning  and  a 
Levenberg-Marquardt  parameter.  Preconditioning  is  often  useful  in  improving  the  convergence  of 
the  conjugate  gradient  algorithm,  and  the  Levenberg-Marquardt  parameter  is  needed  to  stabilize 
the  solution  against  the  effects  of  noise  and  modeling  inaccuracies. 


The  inversion  algorithms  may  require  a  priori  information  about  the  flaw  regions.  The  infor¬ 
mation  can  be  used  to  concentrate  the  inversion  efforts  on  regions  of  interest  rather  than  unflawed 
regions.  Statistical  pattern  recognition  and  computer  vision  techniques  have  been  examined  to 
achieve  this  goal.  The  purpose  of  applying  statistical  pattern  recognition  techniques,  is  to  detect 
the  flaw  regions  and  the  background  regions  in  the  spatial  domain.  In  addition,  a  graphical  tool 
can  be  used  to  analyze  the  raw  data  when  used  as  input  features,  and  evaluate  the  dassifiability 


of  the  measurement  (any  two  features). 
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CHAPTER  I 

DEVELOPMENT  OF  THE  MODEL 


1.  Introduction 

In  this  chapter  we  describe  mathematical  models  that  were  developed  to  be  the  basis 
for  our  inversion  algorithms.  These  models  make  use  of  a  bulk  conductivity  model  of  an 
anisotropic  material.  The  bulk  conductivity  approach  and  the  associated  Green’s  functions 
for  aflat  plate  are  described  in  Sabbagh  Associates’  report  SA/TR-3/88.  Here  we  develop 
direct  and  inverse  models  based  on  a  whip  source  probe  and  a  ring  source  probe.  In 
laboratory  tests  these  probes  have  been  used  successfully  to  detect  flaws  and  so  we  develop 
these  models  to  allow  for  the  reconstruction  of  three-dimensional  flaws. 


2.  Whip  Source  Direct  Model 

(a)  Computation  of  the  Incident  Fields  Due  to  a  Whip  Source 

The  infinitely  long  ‘whip’  is  oriented  parallel  to  the  y-axis,  with  x-coordinate  equal  to 
Xo  and  x-coordinate  equal  to  zq.  Hence,  we  have  for  the  current  density 

Jo(x,y,z)  =  I0S(x-zo)6(z-zo)ay,  -oo  <  y  <  oo.  (1) 

The  Fourier  transform  of  J  is  given  by 

J J_j(x  ~ 


We  have,  from  (15)  of  NSWC  TR  85-304: 

*,(*)  = 

2?T 


(3) 


where  the  tilde  denotes  a  function  of  (*„*,),  and  Gai(*|*o)  the  external  Green’s  func¬ 
tion.  From  here  on  we  suppress  the  subscripts  on  the  Green’s  function  and  replace  them 
with  the  superscript,  ‘(e)’. 

Hence,  the  electric  field  at  level  z ,  within  the  slab,  is  given  by 

E*.(x,  *.) 
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where  the  subscripts  refer  to  a  4  x  3  tensor  (because  there  are  four  transverse  field  com¬ 
ponents  in  e<,  and  the  applied  current  can  point  in  three  directions.) 

Upon  taking  the  inverse  Fourier  transform  of  (4),we  get 

E09(x,y,z)  J°°G^(kz,0-,z,zo)e-^t'Zo)dkt 
=F0x(x0  -x,0-,z,zo) 

T  t°°  ~  ^ 

Eoy(x,y,z)  =  ^  J  G<£(kx,0-,z,z0)e-’km(-t~to)dkx 

=F0y{x0  -  *>0;  z,  z0). 

These  are  the  functions  that  are  to  be  used  in  the  next  section. 

(b)  Computation  of  the  Scattered  Fields  Due  to  the  Whip  Source 
We  will  analyze  the  problem  shown  in  the  figure: 


Field  Point 
~^(x,y,z) 


We  have,  from  (15)  of  NSWC  TR  85-304: 

Ct(z)  =  J  Gia(z|z')  •  («) 

where  the  tilde  denotes  a  function  of  (Jfe,,  k9)\  i.e.,  (6)  is  in  the  Fourier  domain.  Gi2(*|z') 
is  the  “internal  Green’s  function” . 

We  assume  that  the  only  significant  currents  lie  in  the  transverse  plane;  i.e.,  J^(z')  = 
/£B)a,(z')  +  Jl*\(z'  ).  From  here  on  we  suppress  the  subscripts,  12,  on  the  Green’s 


function  and  understand  that  this  is  the  internal  Gro  ' 's  function.  We  also  suppress, 
for  the  moment,  the  independent  variable,  z and  understand  that  refers  tc  he 
anomalous  current  at  some  depth  z1.  Hence,  the  electric  fields  at  the  sensor  level  z  due  to 
the  anomalous  current  at  level  z'  are 


E^GnJl^  +  GuJM 
Ev=G2lJi*  +  GnJ!f\ 


where  the  subscripts  refer  to  a  4  x  3  tensor  (because  there  are  four  transverse  field  com¬ 
ponents  in  e<,  and  the  anomalous  current  can  point  in  three  directions). 

We  are  really  interested  in  Bz  at  the  sensor.  According  to  (10)(b)  of  NSWC  TR 
85-304,  we  have 

Bx  =  fi0Hz  =  ~—Ex  +  —Ey,  (8) 

U)  U) 

which,  according  to  (7),  implies  that 

Bz  =  -^[GnJla)  +  G12Jla)]- 

—kyGn  +  kxG2i  \ 

"  /  * 


■^[521^b)  +  g22^] 

( — kyG\2  +  kxG22  \  ~(a) 

{  -  n ' 


(9)(a) 


For  completeness  we  write  the  other  two  components  of  B: 


Bz  +  noG22Jla) 

By  =HoGA1J^  + 


(9)(*>) 


Let’s  call 


Hxx  =  MoGji,  Hxy  —  fioG22 

Byz  =  M0G4I,  Hyy  =  HqG\2 

Hzz  =  (~kyGn  +  ktG2  i)/u>,  Hzy  =  (—kyGi2  4-  kxG22)/ u>. 
Then,  from  (9)(a): 

Bz(x,y,z)  = 

or,  in  convolution  form: 

£«(*>»,*)  =4^2  JJ  HMX{x  -  -iy,z,z')j£Xt,ij)didii 

+  4 **  JJ  H**(x-t,y-ri',Z,Z')J*a)(t’v)dtdv 


(10) 


(11) 


(12) 
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In  order  to  compute  emf’s,  we  integrate  Bz  over  the  appropriate  sensing  coil.  This 
computation  should  be  done  in  (11),  because  it  is  very  easy  to  integrate  the  exponential  in 
(*,y).  The  result  is  to  introduce  new  H  functions  for  the  emf’s,  and  these  functions  can 
then  be  transformed  for  use  in  the  convolution  integrals  of  (12),  et  seq.  (The  derivation 
of  these  H  functions  is  described  more  fully  in  Appendix  A  of  Chapter  IV;  see  (A. 33),  et 
atq.) 

Let  the  conductivity  of  the  flawed  region  at  level  z'  be  the  scalar  function  a^\x,  y,  z1) 
(for  a  void  region  =  0).  We  are  assuming,  therefore,  that  the  flaws  are  isotropic. 
The  incident  field  at  z'  is  E0(*,y,  z')  =  Eox{z,y,  z')a*  +  Eoy(z,y,  z')a.y.  If  the  source 
of  the  incident  field  is  an  infinitely  long  ‘whip’  oriented  parallel  to  the  y-axis,  with  x 
coordinate  equal  to  q,  then  it  is  straightforward  to  show  that  Eoz(x,y,  z')  =  Fox(q  — 
*,0, z'),  Eoy(x,y,z')  =  Foy(q  —  *,0,2'),  where  Fox  and  Foy  are  given  in  (5).  The  zero  in 
the  second  argument  of  the  F’s  implies  that  the  incident  field  is  independent  of  y,  as  we 
know. 

The  anomalous  conductivity  tensor  is  the  difference  between  the  flaw  conductivity 
tensor  and  the  host  conductivity  tensor: 


V') 

0 

0 

Vn 

ffl2 

0  ' 

0 

i ) 

0 

— 

<J\2 

CT22 

0 

0 

0 

0U) 

m  0 

0 

^3S  . 

Hence, 

4a)(*> y)  =*»?(*.  v)*o.(«  -  *>0,z')  +  <r[l\x,y)F0y(q  -  *,0,z') 

•Jja)(*i V)  =0v,\z,y)FOx{q  -  *,0,z')  +  <r{£\z,y)F0y{q  -  *,Q,z'), 

and  we  note  that  =  012,  which  is  known.  Of  course,  if  we  work  in  the  principal- axis 
system,  then  0\ 2  =  0. 

When  this  is  substituted  into  (12)  we  get 


£*(*,y.*)=^2  JJ  HMx(x-Ciy-rj,ztz,)Fox(q-C>0,z,)a[l\(tri)d(dri 
+  4^2  JJ  H*Ax-t>V-'n\z,*')Foy(q-t,0,z,)<r{2l\t,ii)dtdTi 

+  4^2  JJ  lH**(x-t>y-rnz,z')F0y(q-t,o,z') 

+  Hxy(x  -CV-V,*,  z')Fo *(?  -  t ,  0,  *')]  rj)didrj. 


Because  we  sense  the  field  at  the  same  location  as  the  whip,  x  =  q,  we  can  replace  q 


by  x  in  (14)  and  get 


I 

I 

I 

I 

1 

I 

I 

I 

I 

I 


JJ  Tu(x  —  £,y  —  r}\  z,  z')a\*\t,r))d£dT] 

+  4^2  JJ  Tn(x  ~  t>y  ~  V>z,z')<Ti2l\t,Ti)dtdT)  (15) 

+  4^2  JJ  T12(x  -  (,y  -  r)-,z,z')a[l)(i,T})d£dri> 

where 

Tu(z-t,y-v;z,z')  =H,x(x  -  t,y  -  tj>z,z')Fox(x  -  t,0,z') 

Ti2(z-t,y-V,z,z')  =HZy{x  -t,y  -r)\Z>z')Foy{x  -  £,0,2*) 

Ti2(z  -  t,y-V,z,z')  =Hzx(x  -t,y-  rr,z,z')F0y(x  -  £,0,2') 

+  Hzy(x  -  £,y  -  i)\  z,  z')Fox(z  -  ^,0,2'). 

The  last  integral  in  (15)  is  known,  because  <7^  =  ffi2-  Indeed,  if  we  are  operating  in  the 
principal-axis  system  of  the  host  conductivity  tensor,  then  an  =  0. 

For  three-dimensional  flaws,  we  need  only  perform  a  spatial  integration  of  (15)  with 
respect  to  z'.  We  consider  the  slab  to  be  partitioned  into  Nz  discrete  ‘layers’  and  consider 
a^  to  be  constant  with  respect  to  z'  for  each  layer.  The  transfer  function  t\™\  for 
example,  for  a  layer  bounded  by  2»“^  and  2^  is 


Tn\z-t,V-ri  ;*,*»)=  f  Hzx(x  -  £,y  -  q;  z,z')F0x(x  -  {,0;  2,  z')dz'  (17) 

The  total  field  due  to  flaws  in  a  slab,  then,  is  the  sum  of  the  fields  due  to  flaws  in  each  of 
the  ‘layers’.  This  will  be  fully  explained  in  Section  4. 

Equation  (15)  is  the  same  equation  that  we  derived  for  reconstructing  <r^  using  an 
infinite  current  sheet  for  excitation.  The  only  difference  is  that  for  the  whip  the  transfer 
functions  T\\%  Tjj,  and  7x2  consist  of  the  product  of  two  arrays,  rather  than  the  product 
of  one  array  (Hzx, or  Hzy)  and  a  scalar  (Eqx,oi  Eoy),  as  was  the  case  with  a  current 
sheet  excitation.  Clearly,  if  the  excitation  source  is  bounded  in  the  y-direction,  then 
the  transfer  function  is  the  product  of  two  two-dimensional  arrays.  Hence,  it  appears 
that  we  can  apply  all  of  the  inversion  algorithms  that  were  developed  for  excitation  with 
infinite  current  sheets,  but  that  there  must  be  a  little  more  pre-processing  to  generate  the 
appropriate  transfer  functions. 


3.  Ring  Source  Direct  Model 

(a)  Computation  of  the  Incident  Fields  Due  to  a  Ring  Source 

The  source  is  a  single  filament  ring  parallel  to  the  slab  workpiece.  The  ring  has  radius 
r0  and  is  centered  at  (*o,  yo,  20  )•  The  current  density  is 

Jo(*»  y ,  *)  =  Io6(z  -  zo)(— a*  sin  6  +  a„  cos  0)6(r  -  r0)  (18) 

The  transform  of  J  is  given  by 

j.(*„*„,z)  =  <<*-*,)  Ay  £  (—a*  sin  0  +  ay  cos  6)6(t  —  r0)e’(km*+k*v^dxdy  (19) 
From  Appendix  C  of  NSWC  TR  85-304,  this  transform  is 

Jot*,,*,,*)  =  Jfr-oM*  - f°>  +  a,—)  (20) 

where  K  =  \jkx2  +  kv2. 

From  (15)  of  NSWC  TR  85-304,  the  fields  due  to  this  current  ring  are 
«<(*)  =  J  G{e)(z|zo)  •  J0(zo)«feo 

=  GW(*|*o)  •  (21) 

where  the  tilde  denotes  a  function  of  ( ks,kv )  and  G^'(z|zo)  is  the  external  Green’s 
function  which  is  a  4  x  3  tensor  (because  there  are  four  transverse  field  components  in 
e(  and  the  applied  current  can  point  in  three  directions). 

The  electric  field  at  level  z  within  the  slab,  then,  is 

£c(*,*o)  =  A^i£Aj,(t,ro)  [-i,<5<,?(*i*o)  +  t.eSj,(*i*o)] 

£«,(*,  *o)  =  ^  [-*,<5»  W*>)  +  *.<5g(*l*o)]  (22) 

Taking  the  inverse  Fourier  transform  of  (22)  we  obtain 
Eo*(*,V,z)  = 

^ / /”  K)  +  Jfe.GiVdfe., fcv)j  e-^(--*)+*»^-v»dfcxdfcv 

=  Fo.(*o  -  *,l/o  -  y;  *o) 

Eoy{x, y,z)  = 

J  Foo  [-jfc»<5«(jb-fcv)  +  *«G(2?(ibs,ibv)] 

=  aPov(*o  -  *,yo  -  y;  *o)  (23) 


(b)  Computation  of  the  Scattered  Fields  Due  to  the  Ring  Source 

The  computation  of  the  scattered  fields  due  to  the  ring  source  is  identical  to  that  of  the 
whip  source  (Section  2(b))  through  (12).  For  the  ring  source,  though,  the  incident  field  is 
not  indepc  _dent  of  y .  If  the  ring  source  has  x  coordinate  q  and  y  coordinate  r,  the  incident 
electric  fields  are  Fox(*,y,*')  =  Fo,(q-x,r-y,  z')  and  Eoy(x,y,  z')  =  Foy(q-x,r  -y,  z'), 
where  F o*  and  Foy  are  given  in  (23).  If  we  choose  x  and  y  to  correspond  to  the  principle 
directions  of  the  host  material,  then  oyi  =  a ji  =  0  and  the  anomolous  currents  are 

4a)(*.y)  =  <ruFoz{q-x>r  -y,z') 

4#)(X»V)  =  VnFoyiq  -  x,r  -y>z')  (24) 

When  equations  (24)  are  substituted  into  (12),  we  obtain 
Bt(x,y,z)  = 

“3  j  J ^  H„(x  -t,y-i V,z,  z')F0x(q  z>}t)((,  rjjd^dq 

+  4~a  J  ~rr,*, z')Fov(q  -  (,r  -  (25) 

The  sensor  is  a  single  filament  loop.  If  the  source  ring  remains  stationary  while  the 
sensor  loop  is  moved  around,  then  the  electric  fields  do  not  change  as  the  sensor  loop  is 
moved.  In  this  case,  Fo*  is  a  function  of  {  and  q  only  and  (25)  becomes 

B.{x,y,z)  =  ^  J 

+  (26) 

So,  Bm  is  computed  by 

B,  =  H„Um  +  HtvTLy  (27) 

where 

n,(*,y)  =  Fo«(x,y)aSt)(*,y) 
nv(*,y)  =  Foy(*,y)aS^(x,y) 

On  the  other  hand,  if  the  source  ring  and  the  sensor  loop  always  move  together  (and 
are  concentric),  then  q  =  x  and  r  =  y.  In  this  case,  (25)  becomes 

B,(*,y,z)  =  ^  y  y*  H„{x  -  (ty  -  ti)Fo,(x  -  t,y  -  ri)a[‘\ttri)dtdri 

+  47r*  /  /  **v*(x  ~  t'V~  v)Foy(x  -  Cy  ~  (28) 


and  Bt  is  computed  by 


Bt  =  fnaft  +  fnc#  (29) 

where 

Tn(*,y)  =  H„(x,y)F0x(x,y) 

r22(*,y)  =  Htv(x,y)F0v(x)y) 

Finally,  if  the  source  ring  and  the  sensor  loop  move  together  but  are  not  concentric, 
then  the  offset  must  be  considered  in  (25).  If  the  offset  in  the  x  direction  is  a  and  the 
offset  in  the  y  direction  is  b,  then  equation  (25)  becomes 

B*(*,y,z)  = 

irr5  J  Loo  H“^X  Z')F°*(X  ~t-a>y-V-b>  *  Vnfe  v)d£dT) 

+  4^2  /  Loo  Hys^X  ~^y~‘n'z'  z')F°v(x  -t-a,y-ri-b,  tf  >  v)d£dr]  (30) 
and  Bx  is  computed  by 

Bt  =  (31) 

where 

Ti  .(*,»)  =  «*.(*, v)f£m(*,y) 

T„(z,y)  =  Hn(x,y)  f£">(*,S) 

and 

fit,,)(z,y)  = 

f£">(x,y)  =  ^l-*‘-»F0v(zty) 

For  three  dimensional  flaws,  the  transfer  functions  must  be  integrated  with  respect 
to  z' . 
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4.  Whip  Source  Inverse  Model 

(a)  The  Discretized  Integral  Equation 


We  assume  that  we  are  operating  in  the  principal-axis  system,  so  that  the  last  integral 
in  (15)  vanishes.  If  the  whip  is  oriented  along  the  y-axis,  then  we  know  from  theoretical 
considerations  that  the  electric  field,  Fox(x  —  £,0,  z'),  vanishes.  Thus,  (15)  reduces  to  an 

integral  equation  for  a  single  unknown,  the  anomalous  conductivity,  zf)  at  layer 

z': 

Bz(*, y, z)  =  JJ  T22(x  -t>y-V>2> V, 2')d{d7).  (32) 


This  is  the  contribution  of  the  layer  that  is  located  at  z' .  The  net  contribution  of  all 
layers  is  given  by  the  integral  of  (32)  over  the  flawed  region: 


B(*,y) 


=  {i?  JLj{*  "  * 

=  J  {//”  dz'  ■ 


(33) 


From  here  on  we  will  suppress  the  sensor  2-coordinate,  because  it  is  fixed;  we  also  suppress 
the  subscripts  on  T  and  a^. 

Next,  expand  the  unknown  conductivity  in  a  series  of  pulse  functions  (defined  on  a 
regular  grid  of  spacing  (Sx,Sy,6z)): 

N.  A',  s, 

*)  =  £££  *1 mnP,(*/Sx)P~(y/6,)PM6z),  (34) 

1=0  m=0  n=0 


where  the  {<rjmn}  are  real-valued,  positive  constants.  Any  other  functions  that  lire  defined 
over  this  grid,  such  as  the  sine  functions,  could  work  as  well;  the  pulse  functions,  however, 
and  functions  derived  from  them  by  convolution,  are  particularly  nice. 

The  Fourier  transform  of  (34)  is 


s.  N,  N. 


47T2 


1=0  m= 0  n= 0 


fsin(ktSx/2)^fsin(ky6y/2)^  n 
*  kx6x/2  kySy/2  W*/**)» 


and  when  this  is  substituted  into  (33)  we  get 


SxSy 


N.  N,  N, 


(35) 


( 'v) (‘,'y(w)(T5r) 

e— >1*-  (*—(!+!  /2)«*)-«-  — (»»»+!  /2)*ir)]  Jjb*  ef  fe 


(36) 
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T^n\kZtkt,  u>)  is  the  Fourier  transform  of 


-Ly 


Hzy(x  ~^,y  -Tj-,w,z')Foy{x 


£>  0;a>,  z')dz 


I 

t 


(37) 


which  was  derived  in  (17).  We  are  explicitly  showing  the  dependence  of  the  transfer 
function  on  frequency,  in  anticipation  of  the  multifrequency  model  for  inversion. 

We  take  moments  of  (36)  by  multiplying  by  “testing  functions”,  and  then  integrating. 
For  testing  functions  we  will  use  the  same  pulse  functions,  Pi(x/6x),  P\i(y/6y),  that  were 
used  in  the  expansion  of  the  unknown  (thus,  this  is  Galerkin’s  variant  of  the  method  of 
moments).  The  integration  on  the  right-hand  side  introduces  another  Fourier  transform, 
so  that  the  result  is 


Blm  =  Jf  Pi(x/6x)PM(ky/6y)B(x,y)dxdy 


1—0  m—0  n=0 


sin(fc*5x/2)v  2  ,sin(fcy5y/2)..  2  (38) 
ktSx/2  '  '  ky6y/2  ' 


e-jlt,S*(L-l)+kt6y(M-m)]dkidk^ 


This  equation  can  be  written  as  the  sum  of  Toeplitz  operations 

N.  N.  N, 

=  T(n)(L  - 1> M  -  (39) 

n=0  1=0  m= 0 

where  the  two-dimensional  Toeplitz  matrix,  T^n\  is  given  by  the  integral  in  (38).  This  is 
the  discrete  version  of  the  integral  equation,  (33),  and  is  the  basis  of  our  inversion  method. 

(b)  The  Multifrequency  Model 

Equation  (39)  consists  of  (Nz  + 1)  x  (Nv  + 1)  equations  in  (Nz  + 1)  x  ( Ny  + 1)  x  (Nz  +  1) 
unknowns.  Hence,  we  need  more  equations.  The  easiest  way  to  generate  these  equations 
is  to  repeat  (39)  at  a  number  of  different  frequencies.  This  is  easy  to  do  in  the  lab  (where 
we  must  measure  B  at  these  frequencies),  and  on  the  computer  (where  we  generate  a  new 
transfer  function,  T,  at  the  same  frequencies).  This  approach  is  plausible  because  the 
anomalous  conductivity  is  assumed  to  be  independent  of  frequency.  Thus,  assuming  that 
we  use  Nf  frequencies,  we  have 


B\  —T\%\  ®  ffi  H - 4-  ®  <tn. 


Bnj  =Tnj,\  ®  o\  + - 1-  TN,,Nm  ®  <tn.- 


(40) 


MO 


Each  of  the  B's  in  (40)  is  a  (iVz+l)x(Arv+l)-dimensional  data  array,  the  subscript  denoting 
the  frequency  at  which  the  data  is  taken.  The  operator  ®  denotes  the  two-dimensional 
Toeplitz  operation 


N. 

T,,n  ®  =  S  E  ^"’(i 

1=0  m=0 


(41) 


We  write  (40)  in  a  block-matrix  form,  in  which  the  real  and  imaginary  parts  are 
separated 


■M.1 

ji,nm 

pW 

= 

rp(R) 

T(/) 

t(R) 

•••  I1,N. 

nil) 

BN,  J 

T(>) 

1 

TU) 

O 


<T\ 


0N.  J 


(42) 


Equations  (40)  or  (42)  are  our  basic  system  of  discrete  equations,  and  constitute  the 
‘multifrequency  model’  for  the  whip  source.  These  equations  may  contain  several  thousand 
unknowns,  so  we  must  use  efficient  methods  of  inverting  them.  In  Chapterll  we  will  apply 
the  method  of  conjugate  gradients,  together  with  the  Fast  Fourier  Transform  (FFT),  to 
accomplish  this  task. 

The  question  arises  as  to  the  best  number,  Nf,  of  frequencies  to  use,  and  what  is  the 
optimum  range  of  frequencies.  Generally,  this  can  be  answered  by  trial-and-error,  with 
the  following  ideas  as  a  guide.  Least-squares  methods,  such  as  the  conjugate  gradient 
algorithm,  often  work  better  with  overdetermined  systems,  because  the  variance  of  the 
error  is  reduced  as  the  number  of  equations  increases,  for  a  given  number  of  unknowns. 
Hence,  we  would  like  Nf  to  generally  be  much  larger  than  Nr.  It  is  time  consuming, 
however,  to  generate  too  much  data,  so  there  is  a  trade-off  that  can  only  be  determined 
by  conducting  numerical  experiments  with  typical  problems. 

The  same  can  be  said  in  determining  the  frequency  range,  but  we  know  intuitively,  if 
for  no  other  reason,  that  we  should  use  as  broad  a  frequency  range  as  possible.  We  can  be 
a  little  bit  more  definite  here,  and  rely  upon  the  phenomenon  of  skin  effect  to  guide  us.  If 
we  want  a  resolution  of  SM  in  depth,  then  our  upper  limit  of  frequency  should  produce  a 
skin  depth  that  is  smaller  than  6S,  though,  as  we  will  see  in  Chapter  III,  we  have  gotten 
good  results  with  simulated  data  at  lower  frequencies. 
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5.  Ring  Source  Inverse  Model 

The  inverse  algorithm  developed  in  the  previous  section  for  the  whip  probe  is  equally 
applicable  for  the  ring  probe.  For  the  whip  probe,  we  noted  that  if  the  grid  was  oriented 
parallel  to  the  principle  axes  of  the  workpiece,  =  0  and  the  third  integral  in  (15) 
vanishes.  Also,  since  the  whip  is  also  oriented  parallel  to  the  y  axis,  Fo*  =  0  and  the 
first  integral  in  (15)  vanishes.  Hence,  (32)  includes  only  the  second  integral.  For  the  ring 
probe  it  is  still  true  that  we  can  choose  to  orient  our  reference  grid  parallel  to  the  principle 
axes  of  the  workpieces  and  so  we  can  eliminate  the  third  integral  in  (15).  However,  the 
first  integral  must  be  retained.  The  basic  inversion  equation  for  the  ring  probe,  then, 
involves  7ri  and  <r[“^  as  well  as  Tjj  and  <r^.  For  isotropic  materials,  er$  =  and  the 
transfer  functions  can  be  combined.  For  anisotropic  materials,  we  can  either  treat  the  two 
principle  conductivities  as  independent  unknowns  and  double  the  size  of  the  problem  or 
we  can  assume  a  constant  ratio  for  the  conductivities  and  combine  the  transfer  functions. 
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CHAPTER  II 

APPLICATION  OF  CONJUGATE  GRADIENTS 


1.  Introduction 

The  discretized  system  of  equations,  (1-40),  which  resulted  from  the  application  of  the 
method  of  moments  to  the  operator  equations  in  Chapter  I,  will,  in  general,  have  a  large 
number  of  unknowns.  In  subsequent  chapters  we  demonstrate  some  problems  with  4000  to 
10,000  unknowns.  In  addition,  we  found  that  this  system  has  a  very  special  structure;  it 
was  Toeplitz  in  two  of  the  three  dimensions.  This  means  that  we  should  apply  a  solution 
technique  that  can  accomodate  a  large  number  of  variables,  while,  at  the  same  time,  taking 
advantage  of  the  special  structure.  This  suggests  the  use  of  iterative  techniques,  such  as 
the  conjugate  gradient  (CG)  method. 

We  have  successfully  applied  the  CG  method  to  a  number  of  problems  involving  volume 
integral  equations  in  nondestructive  evaluation  [1-3].  Much  of  the  rest  of  this  report  deals 
with  the  application  of  this  method  to  three-dimensional  inverse  problems. 

In  this  chapter,  we  will  simply  sketch  the  important  features  of  the  CG  algorithm; 
references  [4-6]  should  be  consulted  for  further  details  on  the  method. 

2.  The  Conjugate  Gradient  Method 

Let  us  write  the  complex  vector-matrix  equation  (1-40)  as  the  operator  equation 


where 


and 


Y  =  A  oX, 

‘  Bi  ' 

r  _ 

■  Bn,  . 

r  vi  1 


IVN.  J 


Ti,i  •••  Ti,n. 

’  V\  ' 

AoX  = 

1 

O 

* 

Tn,,i  •••  Tnj.n,. 

■  VN.  . 

(1) 

(2) 

(3) 

(4) 


Keep  in  mind  that  each  of  the  B's  and  <r’s  is  a  two-dimensional  array,  say  of  dimension 
32  x  32,  and  each  Ti.i  in  (4)  is  a  two-dimensional  Toeplitz  matrix,  Tij(l  —  L,m  —  M).  l,m 
index  the  'row’  in  each  of  the  two  dimensions,  while  L,M  index  the  ‘column’  of  each 
dimension. 
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We  will  need  the  adjoint  operator,  A*,  which  corresponds  to  the  conjugate  transpose 
of  the  block-matrix  in  (4): 


■3?i 

tH  1 

^Nj,  1 

■  B a  ■ 

A*  oY  = 

; 

0 

• 

rpH 

.Bn,  . 

T"  is  the  Hermitian  transpose  of  the  two-dimensional  Toeplitz  matrix  Ttj\  i.e.,  TXH}(1- 
L,m  —  M)  =  T*j(L  —  l,M  —  m),  where  *  denotes  the  complex-conjugate.  The  operator 
formalism  is  precisely  the  same  if  we  use  the  real  system,  (1-42),  except  that  the  Hermitian 
transpose  is  replaced  by  the  ordinary  transpose  in  defining  the  adjoint  operator  in  (5). 


We  remind  the  reader  that  the  ©  operation  that  appears  in  (4)  and  (5)  stands  for  the 
sum  of  a  number  of  two-dimensional  Toeplitz  operations,  as  in  (1-40). 


The  conjugate  gradient  algorithm  starts  with  an  initial  guess,  Xo,  from  which  we 
compute  Ro  =  Y  —  A  o  Xo,  P\  =  Qo  =  A*  °  Ro .  In  addition,  we  have  a  convergence 
parameter,  e.  Then  for  k  =  1,...,  if  Test  =  ||-Rt(|/||^||  <  e,  stop;  X*  is  the  optimal 
solution  of  (1).  Otherwise,  update  X*  by  the  following  steps: 

Sk  =  A  o  Pk 


\\Sk\\ 


Xk  =  Xk-\  +  tikPk 


Rk  =  Rk-i  -  akSk 


Qk=  A*  o  Rk 


bk  = 


HQfcll2 

||Qfc-i||2 


•Pfc+l  =  Qk  +  bkPk- 


(6) 


comment  (1):  The  algorithm  terminates  at  the  Mth  step  when  Qm+i  =  0,  so  that 
Xm+i  is  the  least-squares  solution  of  Y  =  A  o  X.  The  vectors  Q 0,  Q lt  are 

mutually  orthogonal,  as  are  the  vectors,  S\,  S2,  5$, _ In  addition 


if  j  <  k ; 
otherwise. 


comment  (2):  This  suggests  that  we  monitor  the  iterates  {(?*}  for  loss  of  orthogonality, 
and  restart  when  the  condition  >  tiQjf+iQk+i  is  satisfied,  where  ei  =  0.2 

(say).  When  this  occurs,  we  set  bk  =  0  in  the  last  line,  and  then  continue  (i.e.,  we  restart 
with  a  pure  gradient  step). 

comment  (3):  Allen  McIntosh,  Fitting  Linear  Models:  An  Application  of  Conjugate 
Gradient  Algorithms ,  Springer- Verlag,  1982,  gives  an  alternative  expression  for  6*: 

.  Q”(Qk  -Qk-i) 

* "  II Q*-i||*  ’ 
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which  seems  to  produce  Qt’s  that  are  more  orthogonal,  when  using  the  criterion  of  com¬ 
ment  (2).  This  definition  requires,  however,  that  an  extra  array  to  store  Qk-i  he  made 
available.  This  is  no  problem  if  we  monitor  for  orthogonality  for  the  purpose  of  restarting, 
because  that  array  is  required  anyway. 

The  convolution  and  correlation  operations  that  are  a  part  of  A  and  A*  are  evaluated 
by  using  the  FFT,  as  described  in  Appendix  A.  This,  together  with  the  fact  that  the 
storage  requirements  are  reasonably  modest,  are  the  reasons  why  the  conjugate  gradient 
algorithm  becomes  attractive  for  large  problems  in  our  model. 

3.  Conjugate  Gradients  with  Constraints 

The  conjugate  gradient  algorithm  that  was  just  described  does  not  constrain  the 
solution.  In  solving  an  inverse  problem,  we  often  need  to  constrain  the  solution  in  order 
to  get  meaningful  solutions.  Hestenes  [6,  Chapter  III]  presents  algorithms  that  involve 
general  linear  inequality  constraints  (such  as  bounds  on  the  solution).  We  show  one  such 
algorithm,  the  active  set  method,  using  Hestenes’  notation. 

The  problem  is  to  minimize  the  quadratic  function, 

F(x)  =  ]-ztAx  —  hTx  -I-  Cg, 

Ct 

on  the  set,  5,  of  all  points,  x,  satisfying  a  set  of  inequality  constraints 

gi(x)  =  w[z-  ki  <  0  (*  =  (7) 

We  assume  that  A  is  N  x  N  and  symmetric;  later,  we  will  consider  the  more  general  case 
that  can  be  solved  using  least-squares. 

A  special  case  of  (7),  which  will  be  of  interest  to  us,  are  the  bounds 

c«<z<<(T  (*  =  1,...,JV). 

This  can  be  put  into  the  form  (7): 

ffi(x)  =  -**  +  c*  <0 
9%+n(*)  =  <0  (*  = 

Hence,  M  =  2N\  tu*  in  (7)  is  the  unit  vector  pointing  in  the  negative  tth- coordinate 
direction,  for  *  =  and  in  the  positive  (»  -  !V)th-coordinate  direction,  for  i  = 

N  +  These  vectors,  of  course,  are  the  outer  normal  vectors  to  the  feasible 

region,  S,  which,  in  this  case,  is  a  cuboid  whose  edges  are  parallel  to  the  coordinate  axes. 

(a)  CG-Algorithm  for  Minimizing  F  on  S  [6,  p.  224] 

Step  1.  Select  a  point  x\  in  S.  Compute 

=  -F'(ii)  -  h~  Axi,  9n=9iM  (t  =  1,...,2W). 
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If  fa  =  0,  stop;  algorithm  terminates. 

Else,  let  I  be  the  set  of  indices,  i,  such  that  gn  =  0  (the  ‘active  6et’),  and  go  to  Step  2. 

Step  2.  If  I  is  empty,  i.e.,  no  constraints  are  active,  let  H  be  the  identity  matrix,  and  go 

to  Step  3. 

Else,  let  H  be  the  nonnegative  symmetric  matrix  that  annihilates  the  vectors,  tu»,  »  e  I. 

comment:Note  that  if  the  first  N  constraints  are  active,  or  the  last  N,  then  H  is  the 
sero  matrix,  because  that  is  the  only  matrix  that  jointly  annihilates  all  of  the  coordinate 
vectors.  This  is  the  case  when  x  is  at  a  vertex  of  the  feasible  region  (which  corresponds  to 
a  comer  of  the  cuboid).  In  general,  the  columns  of  H  are  orthogonal  to  the  normal  vectors 
of  the  surfaces  that  intersect  to  form  the  part  of  the  cuboid  on  which  *  lies.  (Consider  the 
case  in  which  x  lies  on  a  face  of  the  cube,  or  on  an  edge.)  It  will  be  quite  apparent  how 
we  apply  H ,  so  that  no  matrix  multiplies  will  be  involved. 

If  H  =  0  go  to  Step  5,  with  Xi  playing  the  role  of 

Else  go  to  Step  3. 

Step  3.  CG-subroutine.  Set 

Pi  =  fi  =  Ht! .  (a) 


comment  :If  H  is  not  the  identity  matrix,  then  when  it  operates  on  a  vector  it  merely  nulls 
certain  components.  Thus,  the  resulting  vector  lies  in  the  constraint  subspace  in  which  z 
is  located.  This  means  that  the  correction  vector,  p,  lies  in  the  constraint  boundary,  if  z 
starts  there. 

Starting  with  k  =  1  compute 


Cfc 


Else 


»k=Apk,  c*  =  ptrk,  *k  =  Pk»k}  a*  =  -T-, 

a* 

n  -I0’  *  G  1 

q  “  Wm,  iil, 

*fc+!  =  *k  +  ajtpfc,  9j,k+i  =  9jk  +  (j  =  1,...,2JV). 

If  for  some  j  ^  I,  gj,k+ a  >  0,  go  to  Step  4. 

»■*+!  =  »“*  —  f't+i  =  Hrk+i . 


If  ffc+i  =0,  or  k  =  N,  go  to  Step  5. 

Else 

l  *Ifk+ 1 

Pk+i  =  nt+i  +  f>kPk,  ok  = - j- — . 

Replace  Jfe  by  k  +  1  and  go  to  (b). 


(*) 

(«) 

(d) 

(«) 

(/) 

(17) 

(h) 
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Step  4.  Scale  back  to  the  boundary  of  the  feasible  region  and  update  the  active 
set.  Let  J  be  all  indices  j  I,  3  g},k+ 1  >  0.  Let  fi*  be  the  smallest  of  the  ratios 


Reset 


ajk  = 


i  €  J- 

9jk 


*i  =  *k  +  O’kPk,  ri=h-Ax!,  g}1  =  gj<k  -f  akqjk,  (j  =  1,  •  •  • ,  2N). 

Update  the  active  set  by  adjoining  to  I  all  indices  j  $  I,  3  g3\  =  0.  Then  go  to  Step  2. 

comment:  Instead  of  scaling  back  to  the  boundary  along  the  conjugate  direction,  we  can 
return  to  the  boundary  along  the  orthogonal  direction  from  the  current  estimate  point. 
This  seems  to  speed  up  the  algorithm. 

Step  5.  If  rfc+i  =  0,  stop;  algorithm  is  terminated. 

Else  select  shortest  v  of  the  form 

v  =  r*+ 1  —  Wiyi ,  summed  for  all  i  6  I  with  yi  >  0.  (*') 


comment: It  is  easy  to  do  this  because  the  wt  are  unit  coordinate  vectors.  The  yi  are 
celled  either  Lagrange  multipliers  or  components  of  the  dual  vector. 

If  -u  =  0,  stop;  Kuhn-Tucker  conditions  axe  satisfied,  and  algorithm  is  terminated  at 
minimum  point  of  F  on  5. 

Else,  choose  a  >  0,  such  that 

3j,fc+i  +  avTWj  <0  (i  =  1, . . . ,  2 N),  F(xk+i  +  av)  <  F(xk+i ).  (j) 

comment :The  last  inequality  holds  when  0  <  a  <  2vTrk+\/vTAv. 

Restart  the  algorithm  at  Step  1  with 


*i  =  *k+i  +  av 


as  the  initial  point. 

(b)  CG- Algorithm  for  Least-Squares  on  5 

comment:The  algorithm  follows  from  the  preceding  one,  after  observing  that  the  new 
functional  to  be  minimized  is 

=i XTA‘AX  -  Xt(A-Y)  +  i||yf , 
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where  we  are  returning  to  our  usual  notation.  The  constraints  are  the  same  as  before,  (8). 
Step  1.  Select  a  point  .Xi  in  5.  Compute 


Rx  =  Y-AXu  Q1=A'R1  =  -F'(X1 ),  9il  =  gi(X  i)  (*  =  1,...,2N). 

If  Qx  ss  0,  stop;  algorithm  terminates. 

Eke,  let  I  be  the  set  of  indices,  t,  such  that  gn  =  0  (the  ‘active  set’),  and  go  to  Step  2. 

Step  2.  If  /  is  empty,  i.e.,  no  constraints  are  active,  let  H  be  the  identity  matrix,  and  go 

to  Step  3. 

Else,  let  H  be  the  nonnegative  symmetric  matrix  that  annihilates  the  vectors,  Wx,  *  6  I. 
comment :This  is  the  same  matrix  as  above. 

If  H  =  0  go  to  Step  5,  with  Xi  playing  the  role  of  Xk+i- 
Else  go  to  Step  3. 

Step  3.  CG-subroutine.  Set 


P1=751  =  HQi. 


Starting  with  k  =  1  compute 


5*  =  APt,  at  = 


]gyi 

M2 


lik-\WTpk>  it  I, 


Xt+i  =  Xk  +  akPk,  9j,k+ 1  =  9jk  +  0  =  1,...,  2  N). 

If  for  some  j  1 9j,k+ i  >  0,  go  to  Step  4. 


Rk+i  =  Rk  —  &kSk,  Qk+l=A*Rk+l ,  (/) 

If  <Jfc+1  =0,  or  k  =  N,  go  to  Step  5.  ( g ) 

Else  _ 

h  =  ■  p*+>  =  $»+>  +  wv  (fc) 

Replace  Jb  by  k  +  1  and  go  to  (b). 

_ T — 

comment:  The  definitions  of  a*,  bk  force  the  orthogonality  conditions:  QkQk+ i  =: 
0,  STSk+i  =  0.  This  suggests  that  we  monitor  the  iterates  {Qk}  for  loss  of  orthogonality, 

*  —7’—  — j*  — 

and  restart  when  the  condition  >  eiQfc+i^fc+i  *8  satisfied,  where  ti  =  0.2 

(say).  When  this  occurs,  we  set  6*  =  0  in  (h),  and  then  continue. 
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Step  4.  Scale  back  to  the  boundary  of  the  feasible  region  and  update  the  active 
set.  Let  J  be  all  indices  j  I,  9  9j,k+i  >  0.  Let  a*  be  the  smallest  of  the  ratios 

“,k  =  -fii,  j  e  J. 

Reset 

-Xi  =  Xk  +  atPk,  Ri=Y  —  AXi,  Qi=A*Ri,  gji  =  gjtk  +5*//*,  (j  =  1,. . .  ,2N). 

Update  the  active  set  by  adjoining  to  7  all  indices  j  £  7,  3  =  0.  Then  go  to  Step  2. 

comment:  Instead  of  scaling  back  to  the  boundary  along  the  conjugate  direction,  we  can 
return  to  the  boundary  along  the  orthogonal  direction  from  the  current  estimate  point. 
This  seems  to  speed  up  the  algorithm. 

Step  5.  If  Qk+i  =  0,  stop;  algorithm  is  terminated. 

Else  select  shortest  V  of  the  form 

V  =  Qk+i  —  ^  WjAj,  summed  for  all  i  £  7  with  \i  >  0.  (i) 

If  V  =  0,  stop;  Kuhn-Tucker  conditions  are  satisfied,  and  algorithm  is  terminated  at 
minimum  point  of  F  on  S. 

Else,  choose  a  >  0,  such  that 

giMl+aVTWj<0  (j  =  1, . .  .,2N),  F(Xk+x  +  aV)  <  F(Xk+1).  (;) 

comment  :The  last  inequality  holds  when  0  <  a  <  2VTQk+i/VTA*AV. 

Restart  the  algorithm  at  Step  1  with 


X\  =  Xk+1  -f  a.V 


as  the  initial  point. 

(c)  Example  Calculation 


r°  ro  i  ii  pi 

4.  i  i  ,  ,  r.  o  . 

11  L  J  8 


Then  the  solution  of  the  normal  equation 


isX  = 


.  The  outer-normal  vectors  to  the  two-dimensional  constraint  region  (which  is 


a  square)  are 


Wi  = 


Start  with  X\  = 


-1 

0 

0 

0 


,  W2  = 


0 

-1 


,  WS  = 


,  = 


Hence,  the  active  index  se 


.  Then  9l(Xi)  =  0,  g^X,)  =  0,  5s(*i)  =  -1,  9i{X,)  =  -1. 
is  I  =  (1,2),  which  means  that  the  H  matrix  is  the  null-matrix. 


The  initial  residual  and  gradient  vectors  are,  respectively, 


'2' 

Ri  = 

0 

8 

,  Q1=A*R1  = 

'  8  ' 
10 

Enter  5:  Minimize  V  with  nonnegative  A2,  where, 


V  = 


8 

10 


Hence,  Aj  =  A2  =  0,  and  V  = 


-Aj 

8 

10 


-1 

0 


—  a2 


0 

-1 


8  +  Ai 

10  +  A2 


.  Next,  consider 


0  <  a  < 


51.1  +  a[8  10] 

52.1  +  «(8  10] 

5s,i  +  ®(8  10] 

5<,i  +  <*[8  10] 

328 


r  ~  1 


-1 

0 

0 

-1 

1 

0 

0 

1 


=  —  8a  <  0  =>  a  >  0 
=  —  10a  <  0  =>  a  >  0 
=  —  1  +  8a  <  0  =»  a  <  1/8 


=  -  1  +  10a  <  0  =»  a  <  1/10 


328 


102  -I-  182  +  182  748' 

Thus,  a  =  0.1,  and  we  leave  step  5  with 


*i  = 


0.8 

1.0 


0.8 

1.0 


Enter  step  1  with  this  value  of  X\  and  compute 


'2' 

'o  r 

_ 

'  1  * 

Ri  = 

0 

8 

— 

i  i 
i  i 

0.8 

1.0 

= 

-2.8 

5.2 
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Enter  step  2: 


H  = 


0  O' 

0  0 


Enter  step  5:  Find  smallest  V  with  nonnegative  As,  A4,  such  that 


V  = 


4  -  As 

5  -  A4 


Hence,  As  =  4,  A4  =  5,  and  V  = 


Hence,  stop  with  solution 


X  = 


A  = 


4.  Preconditioned  Conjugate  Gradients 

The  rate  of  convergence  of  the  conjugate  gradient  method  depends  upon  the  condition 
number  of  the  matrix  operator,  A.  Hence,  in  order  to  speed  up  the  convergence,  we 
precondition  A.  There  are  several  ways  to  do  this  (see  Allen  McIntosh,  Fitting  Linear 
Models:  An  Application  of  Conjugate  Gradient  Algorithms,  New  York:  Springer- Verlag, 
1982);  we  are  going  to  consider  only  one  method,  scaling  to  produce  columns  of  A  that 
are  unit  vectors. 

Return  to  the  basic  equation 


Y=AX 

—ABB~lX  (9) 

=ABU, 

where  B  is  some  invertible  operator  (matrix),  and  U  =  B~1X.  B  is  chosen  to  improve  the 
condition  number  of  A.  We  will  take  it  to  be  a  diagonal  matrix,  whose  nth  entry  is  the 
reciprocal  of  the  norm  of  the  nth  column  of  A. 

We  do  not  multiply  A  with  B,  because  that  would  destroy  the  special  convolutional 
structure  of  A.  We  introduce  a  new  operator,  A1,  which  is  the  composition  of  B  followed 
by  A.  The  adjoint  of  A'  is  then  A'*  =  BTA*.  Note  that  because  B  is  diagonal,  then 
BT  =  B.  In  multidimensional  problems,  in  which  we  don’t  actually  write  out  the  operator 
equations  as  matrix  equations,  we  intrepret  a  “diagonal  operator”  to  be  an  operator  that 
multiplies  the  nth  component  of  the  solution  vector  by  the  nth  component  of  B. 

The  conjugate  gradient  algorithm  remains  unchanged,  except  that  A  is  replaced  by  A'. 
The  solution,  U,  of  the  scaled  equation  is  then  unsealed  to  get  the  original  solution:  X  = 
BU.  Keep  in  mind  that  this  also  introduces  scaled  bounds  in  the  constrained  conjugate 
gradient  algorithm  for  U.  That  is,  if  c1  <  x%  <  <P,  then  ( c*/B% )  <  u’  <  (d*/5‘). 
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5.  Conjugate  Gradients  with  the  Levenberg-Marquardt  Parameter 

In  our  discussion  of  deconvolution  via  Fourier  transforms,  we  introduced  a  filter  pa¬ 
rameter,  a,  which  smoothed  the  solution.  Here,  we  introduce  a  similar  parameter,  A, 
called  a  Levenberg-Marquardt  parameter.  We  want  to  see  how  it  enters  into  the  conjugate 
gradient  algorithm. 

We  start  with  the  augmented  functional 
F(X)  =i||y  -  AX ||’  + 

=\xTA’AX  -  XtA-Y  +  i||VT  +  iA’pe l|’  (10) 

=\xta!‘ax  -  xTA-y  +  i|im2, 

where  the  augmented  operator,  A' ,  is  defined  by 


and  the  augmented  right-hand  6ide  is 


(12) 


I  in  (11)  is  the  N  x  N  identity  matrix  (where  N  is  the  number  of  unknowns),  and  the  zero 
vector  in  (12)  is  of  length  N. 

The  vectors  in  the  conjugate  gradient  algorithm  that  are  affected  by  these  definitions 
are  S  and  Q.  We  define  new  vectors,  S'  and  Q' ,  in  terms  of  the  old  ones  by 


Q'  =  -F'(X)  =A-Y-{A'A  +  \2)X 

=A*R-\2X  (14) 

=Q-X2X. 

It  is  easy  to  see  that  ||S'||2  =  ||S||2  +  A2||P||2. 

We  can  use  the  Levenberg-Marquardt  parameter  to  solve  a  constrained  least-squares 
problem.  First,  let  us  approximate  the  cuboidal  constraint  set,  that  was  defined  in  Section 
3,  by  the  inscribed  hypersphere,  whose  center  is  at  the  centroid,  Xe,  of  the  cuboid.  Then 
we  replace  the  functional  of  (10)  by 

F(X)  =  |||K- >U||’  +  jA’||*  -  *t||’,  (15) 
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which,  under  the  change  of  variables,  U  =  X  —  Xc,  becomes 

0(10  =  jl|y  -AX.-  AVf  +  (16) 

Hence,  the  previous  algorithm  is  unchanged,  except  that  the  inhomogeneous  term,  Y,  in 
(12)  is  replaced  by  Y  —  AXe. 

Now  the  question  arises,  how  do  we  determine  A?  In  a  statistical  approach  to  de- 
convolution,  A  can  be  given  in  terms  of  signal-to-noise  ratio,  as  is  done  in  Chapter  IX. 
Without  this  data,  however,  we  look  for  something  else.  In  our  discussion  of  decon¬ 
volution  via  Fourier  transforms,  we  showed  that  the  smoothing  parameter,  a,  could  be 
determined  by  solving  a  nonlinear  equation,  given  certain  prior  information.  We  do  some¬ 
thing  similar  with  the  Le venberg-  M arquardt  parameter,  but  first  we  quote  the  following 
theorem:(Charles  L.  Lawson  and  Richard  3.  Hanson,  Solving  Least  Squares  Problems,  En¬ 
glewood  Cliffs:  Prentice-Hall,  Inc.,  1974,  p.  193) 

Theorem.  For  a  fixed  nonnegative  value  of  A,  say,  \,  let  X  be  the  solution  vector  for  the 
problem  of  minimizing  (10),  and  let  fi  =  j||K  —  «4J£||2.  Then  Tl  is  the  minimum  value  of 
j|| Y  —  AX ||2  for  all  vectors  X  satisfying  |JXJ(  <  ||X||. 

The  proof  of  this  theorem  is  simple,  and  is  given  in  Lawson  and  Hanson.  The  meaning 
of  the  theorem  is  clear:  Given  the  radius  of  the  hypersphere  constraint  set  (which  follows 
from  the  original  hypercube  constraint  set),  j|X||,  we  minimize  (10)  by  means  of  conjugate 
gradients,  for  a  collection  of  A’s.  We  choose  that  solution  for  which  ||X||  =  ||X||.  This 
yields  the  optimum  A,  and  gives  us  the  optimum,  constrained,  least-squares  solution  of 
Y  =  AX. 

This  approach  is  equivalent  to  solving  for  A  using  a  trial-and-error  method.  This  is 
inevitable  when  using  conjugate  gradients.  If  we  were  solving  a  much  smaller  problem,  we 
could  determine  the  singular  value  expansion  of  A,  and  use  that  result  in  setting  up  an 
analytic  equation  for  A,  which  could  then  be  solved  using  Newton’s  method,  as  indicated 
in  Lawson  and  Hanson.  It  remains  to  be  seen  whether  this  approach  is  faster  than  that 
which  constrains  each  component  of  the  solution  vector  individually. 

There  is  a  serious  problem  with  this  approach,  however.  The  inscribed  hypersphere  in 
N-dimensional  space  has  a  volume  that  is  much  smaller  than  the  N-dimensional  constraint 
cuboid,  when  N  is  only  reasonably  large,  whereas  the  circumscribed  hypersphere  has  a  vol¬ 
ume  that  is  much  larger.  Hence,  using  the  inscribed  hypersphere  results  in  a  very  strong 
(and  undoubtedly  incorrect)  constraint,  whereas  using  the  circumscribed  hypersphere  re¬ 
sults  in  a  very  weak  constraint. 
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APPENDIX  A 


Efficient  Computation  of  Convolutions  and  Correlations 


If  we  attempt  to  solve  (1)  using  an  iterative  technique,  such  as  the  conjugate  gradient 
method,  it  will  be  necessary  to  evaluate  the  vector-matrix  product  many  times.  This  is  a 
PQ- step  operation,  where  P  is  the  number  of  rows  (equations),  and  Q  is  the  number  of 
columns  (unknowns)  of  the  matrix  T.  If  the  matrix  is  square,  with  dimension  N  x  N,  this 
process  can  be  reduced  to  Nlog2  N  operations  by  using  the  Fast  Fourier  Transform  (FFT) 
for  evaluating  discrete  convolutions. 

The  appropriate  theorems  (in  one  dimension)  that  relate  discrete  Fourier  transforms 
and  convolutions  and  correlations  are  (  <=>  denotes  a  discrete  transform-pair): 


If  g{j)  «=>G(n) 

Mi)  ++  H(n) 


.  N- 1  -  N- 1 

Then  X)  Mi  +  *)M*)  =  jj  53  M*)M*  -  i) 

t=0  k= 0 

<=»  G(n)H(—n) 
=G(n)H(N  -  n ) 

^  E  sWMj + *)  =4  E -  >>(*) 

k=0  k=0 

«=!>  G{-n)H{n) 

=G(N  -  n)H(n) 


(A.l  )(a) 


(A.l)(6) 
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(A.l)(c) 


N- 1 


N- 1 


ji  E  «<*)*(;  -  *)  =4  E  so  -  *w*) 

4=0  4=0 

■<=>  C7(n)/ir(n), 

where  j  =  0,  •  •  • ,  N  —  1,  n  =  0,  •  •  • ,  JV  —  1  in  all  of  these.  Several  points  should  be  made: 
first  note  that  correlation  summing  is  not  commutable,  and  that  one  must  use  negative 
frequencies  in  the  discrete  Fourier  transform  (which,  of  course,  introduces  the  term  N  —  n). 

Let’s  look  at  the  matrix  structure  of  convolution  and  correlation  sums  and  see  how 
we  can  use  FFT  techniques  to  compute  them.  We’ll  work  in  one  dimension.  Consider  the 
following  convolution  sum,  which  is  written  as  a  vector-matrix  equation: 


(A.2) 


Rewrite  this  in  the  expanded  form  (padding  with  zeros  to  get  a  power  of  two)  in  order  to 
achieve  a  circulant-matrix: 


'yo' 

‘nto  m_i  m_2  m_j " 

■*o‘ 

yi 

mi  m0  m_i  m_2 

*1 

y2 

m2  mi  mo  m_i 

*2 

•ys. 

.ms  m2  mi  mo  . 

■  *s  ■ 

'yo' 

‘  m0 

m_i 

m_  2 

m_s 

0 

m3 

m2 

mi 

'*0' 

yi 

mi 

m0 

m_i 

m_2 

m_s 

0 

m  3 

m2 

*1 

ya 

m2 

mi 

m0 

m_i 

m_2 

m_s 

0 

m3 

*2 

ys 

m3 

m2 

m  1 

mo 

m_i 

m_2 

m_s 

0 

*s 

* 

0 

m  3 

m2 

mi 

mo 

m_i 

m_2 

m_s 

0 

* 

m_s 

0 

ms 

m2 

mi 

mo 

m_i 

m_2 

0 

* 

m_2 

m_  3 

0 

ms 

m2 

mi 

mo 

m_i 

0 

.  * . 

m_i 

m_2 

m_s 

0 

m3 

m2 

mi 

mo  . 

.  0  . 

(A.3) 


where  the  *  denotes  a  discarded  entry.  Hence,  the  sequences  to  be  FFT’d  are:  (mo ,  irij ,  m2 , 
Tnj,0,m_j,Tn_2,m_i)  and  (*Oi*ij*2,*s,0,0,0,0),  and  the  output  sequence  is  (yo>yi,y2,ys, 
*,  *,  *,  *).  The  order  of  the  entries  in  the  sequences  is  very  important. 


Now  for  correlations: 


yo 
Vi 
Vi 
Vi 
* 
* 
* 
l  * 


'yo' 

‘mo 

mi 

m2 

mg' 

■*o' 

yi 

mi 

m2 

m3 

m4 

ya 

m2 

m3 

m4 

m® 

*2 

-ys. 

.m3 

m4 

mg 

m® . 

■  *s  - 

ided  form  (padding 
k: 

with 

zeros  to  get 

a  power  of 

m0 

mi 

m2 

ms 

m4 

ms 

m® 

0  ■ 

■*o' 

mi 

m2 

m3 

m4 

mg 

m® 

0 

mo 

*1 

m2 

m3 

m4 

m  g 

m® 

0 

mo 

mi 

*2 

m3 

m4 

mg 

m® 

0 

mo 

mi 

m3 

*3 

m4 

m5 

m® 

0 

mo 

mi 

m2 

ms 

0 

m5 

m® 

0 

mo 

mi 

m2 

ms 

m4 

0 

m® 

0 

m  0 

mi 

m2 

m3 

m4 

mg 

0 

0 

mo 

mi 

m2 

m3 

m4 

m® 

m®. 

.  0  . 

(A.4) 


(A.5) 
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where  the  *  denotes  a  discarded  entry.  Hence,  the  sequences  to  be  FFT’d  are:  (mo  >7711,1712, 
rnj, 7714,1715, m«,0)  and  (*o,*i,*j,*j,0,0,0,0),  and  the  output  sequence  is 

*»  *>  *)■  The  order  of  the  entries  in  the  sequences  is  very  important,  and  also  don’t  forget 
to  negate  the  frequencies  in  the  transform  of  the  z-sequence. 

To  summarize:  we  expand  the  original  data,  padding  with  zeros,  as  necessary,  to  get 
a  circulant  matrix,  and  then  take  FFT’s. 

APPENDIX  B 

Bi-Conjugate  Gradients 

Let  us  consider  the  following  equation, 

AX  =  Y,  (B.l) 

where  A  is  a  known  operator,  X  is  the  unknown,  and  Y  is  the  known  vector.  For  a 
non-Hermitian  operator,  the  conjugate  gradient  method  solves  the  normal  equation, 

A*  AY  =  A*Y,  (B.2) 

where  A *  is  the  adjoint  operator  of  A.  It  is  noted  that  the  condition  number  of  the  original 
equation  (B.l)  is  squared  in  the  solution  of  (B.2).  In  the  bi-conjugate  gradient  method,  one 
solves  the  non-Hermitian  operator  equation  (B.l)  directly.  For  this  algorithm  an  additional 
2N  storage  spaces  is  required,  where  N  is  the  dimension  of  the  unknown  vector  X. 

Method  of  Solution 


To  solve  (B.l),  one  starts  with  an  initial  guess,  Xo,  for  X,  and  then  defines  the 
residual  Rq  =  Y_—  AX o,  and  an  initial  search_direction  P0  =  Rq.  In  parallel,  we  define 
the  bi-residual  Rq  =  -ftjj  and  bi-directional  Po  =  P0*,  where  (*)  denotes  the  complex 
conjugate.  In  addition,  we  have  a  convergence  parameter,  e.  Then  for  A:  =  1,2, ...,  if 
Test=  ||iZjfc|j/||y||  <  e,  stop.  Otherwise,  update  X*  by  the  following  steps: 


Sk  =  APk 

ak  = 

<Sk,Pk> 

Xk+i 

Pk+l  =  Rk+l  +  bkPk 
•Rfc+i  =  Rk  —  akSk 


Sk  =  A'Pt 
bk  = 

<Rh,Rk> 

Xk  +  akPk  _ 

Pk+ 1  =  R*+i  + 

Rk+ 1  =  Rk  —  o,\. Sk- 


The  scalar  {at}  is  chosen  so  as  to  force  the  bi-orthogonality  conditions, 


(B.3) 


<  Rn,Rm>=<  Rn,Rm>=  0  I  0  <  T»  <  n  <  N, 


and  {6t}  i*  chosen  to  force  the  bi-conjugacy  conditions, 


<  Pn,Sm  >=<  Sn,Pm  >=  0  ;  0  <  m  <  n  <  JV, 
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provided  that  the  algorithm  does  not  break  down;  i.e.,  for  all  *  for  which 

<Si,Pi>^  0  ;  <Ri,Ri>^  0. 

The  algorithm  must  terminate  with  Rn  =  Rn  =  0  in  at  most  N  iterations.  It  is 
important  to  point  out  that  the  bi-conjugate  gradient  algorithm  does  not  minimize  the 
norm  of  the  residual,  ||/2||,  at  each  iteration.  Nevertheless,  if  the  algorithm  does  not  break 
down,  it  converges  at  a  much  faster  rate  than  the  normal  conjugate  gradient  algorithm. 
Prom  our  experience  with  large  size  and  ill-conditioned  operators,  however,  the  bi-conjugate 
gradient  algorithm  does  seem  to  break  down;  (i.e.,  |  <  Si,  Pi  >  |  <  6  ;  |  <  Ri,Ri  >  |  < 

6,  where  6  is  a  very  small  positve  number),  before  it  meets  the  convergence  criteria. 


CHAPTER  III 

RECONSTRUCTIONS  USING  SIMULATED  AND  LABORATORY  DATA 


1.  Introduction 

The  goal  of  our  work  is  to  reliably  reconstruct  three  dimensional  flaws  in  workpieces 
from  emf  measurements  that  are  inherently  noisy.  First,  though  we  tested  the  algorithm 
using  simulated  data,  that  is,  computer  generated  data.  Through  these  tests,  we  were 
able  to  learn  how  the  conjugate  gradient  (CG)  algorithm  performs  as  a  function  of  the 
frequency  range  of  the  emf  data  and  the  degree  of  overdeterminedness  of  the  system. 
Also,  the  equations  we  are  trying  to  solve  are  ill-conditioned  and  so  we  needed  to  know  if 
acceptable  solutions  could  be  expected  using  this  form  of  the  CG  algorithm.  The  results 
of  the  tests  using  simulated  data  are  very  encouraging.  We  also  include  some  preliminary 
results  based  on  actual  measured  data. 

2.  Reconstructions  Using  Simulated  Data 

(a)  Simulating  the  EMF  Data 

Two  types  of  synthetic  data  were  available.  Since,  in  the  fourier  domain,  we  are 
solving  an  equation  of  the  form  Ax  =  B,  we  can  provide  a  test  solution  xt  and  calculate 
Bt  —  Axt.  Then,  using  A  and  Bt,  we  can  try  to  reconstruct  xt  using  the  CG  algorithm. 
This  type  of  synthetic  data  will  be  referred  to  as  “exact”  data. 

Another  type  of  synthetic  data  resulted  from  our  model  of  the  direct  problem.  In 
the  direct  problem,  we  calculate  emf  data  from  a  flaw  definition  and  a  transfer  func¬ 
tion  derived  from  internal  and  external  Green’s  functions  (see  Chapter  I).  This  transfer 
function  is  discretized  using  the  method  of  moments  and  becomes  the  operator  A  in 
the  inverse  model.  The  discretization  introduces  error  so  that  the  emf  calculated  using 
the  direct  model  differs  slightly  from  the  emf  calculated  using  the  discrete  operator  A. 
The  difference  between  the  two  diminishes  as  more  terms  are  used  to  approximate  the 
operator.  This  type  of  synthetic  data  will  be  referred  to  as  “direct  model”  data. 

All  tests  using  simulated  data  used  the  “direct  model”  data. 

(b)  Material  and  Data  Collection  Parameters 

The  emf  that  is  measured  depends  on  the  material  parameters  and  and  the  parameters 
of  the  data  collection  system.  Figure  1  shows  the  typical  arrangement  of  workpiece, 
source  (whip)  and  sensor.  For  the  synthetic  data  tests  the  distances  shown  were  as 
follows. 


2*ource  =  0.013  in. 
z»mtcr  =  0.019  in. 
Zr.turn  =  0.273  in. 
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The  sensor  consisted  of  ten  square  turns  from  0.10  inch  to  0.55  inch.  The  workpiece 
was  0.11  inch  thick.  The  conductivity  tensor  was  (mho6/xn) 

'  104  0  0  ' 

0  104  0 
0  0  1. 

This  represents  a  graphite-epoxy  composite  material  that  is  isotropic  in-plane  and  essen¬ 
tially  nonconducting  through  its  thickness. 

Data  measurements  were  simulated  at  each  point  of  a  32  x  32  grid  in  the  sensor  plane 
at  a  resolution  of  0.1  inch.  The  thickness  was  discretized  into  four  layers,  each  also  a 
32  x  32  grid.  A  test  flaw,  then,  could  be  defined  as  a  collection  of  “voxels”  each  with 
dimension  0.1  x  0.1  x  0.0275. 

(c)  Convergence  Test 

A  convergence  test  variable  was  defined  in  terms  of  the  2,a  norms  of  the  simulated 
data  and  the  residual.  Specifically,  at  iteration  t,  the  convergence  test  variable  6,  is 

,  r-Ajqi, 

1 "  uni, 

where  Y  is  the  simulated  emf  data,  Xi  is  the  flaw  function  at  iteration  t  and  Y  —  AXi  is 
the  residual. 

(d)  Test  Set  #1 

A  flaw  in  the  shape  of  a  cross  was  chosen  as  the  standard  test  flaw.  This  test  flaw 
was  placed  in  the  center  of  each  of  the  four  layers  (one  at  a  time)  and  an  attempt  was 
made  to  reconstruct  the  flaw  from  emf  data.  The  conductivity  was  scaled  so  that  flaw 
locations  had  a  value  of  1  and  host  material  locations  had  a  value  of  0.  For  ease  of 
reference,  we  will  refer  to  a  particular  flaw  using  the  word  ‘flaw’  followed  by  digit(s) 
indicating  the  layer(s)  that  contain  the  test  flaw.  This  first  test  set,  then,  involves  fiawl, 
flaw2,  flaw3  and  flaw4.  Data  was  simulated  at  5  frequencies:  1,  3,  5,  7  and  9  MHz.  The 
unconstrained  version  of  the  CG  algorithm  was  used  and  the  initial  guess  was  no  flaw, 
that  is,  all  zeroes. 

The  number  of  unknowns  in  this  problem  is  32  x  32  x  4  =  4096.  Because  the  emf 
data  is  complex  and  our  solution  (conductivity)  is  real,  each  set  of  emf  data  provides 
2  x  32  x  32  =  2048  equations.  For  5  sets  of  emf  data  and  4  layers,  the  overdeterminedness 
of  the  system  is  2.5. 

The  purpose  of  this  set  of  tests  was  to  determine  the  ability  of  the  CG  algorithm  to 
isolate  the  flaw  to  the  proper  layer.  Also,  we  wanted  to  find  out  the  effect  of  depth  on 
the  convergence. 

Figure  2  shows  grayscale  plots  of  the  attempted  reconstruction  of  fiawl  at  four  stages 
during  the  iterative  process:  100,  500,  1000  and  2048  iterations.  Table  1  below  shows 


I 

|  the  value  of  the  convergence  measure,  as  well  as  the  minimum  and  maximum  values 

of  the  solution  for  each  layer. 


Table  1 


|  flawl 

Layer  1 

Layer  2 

Layer  3 

Layer  4 

Iterations 

* 

max 

max 

H2S 

max 

min 

max 

100 

0.00397 

1.343 

lljWjJ 

500 

0.00067 

R 

1.216 

R 

RBKexI 

? 

1000 

0.00027 

-0.293 

SBI 

-0.065 

JRfl 

R 

0.027 

2048 

0.00011 

1.364 

-0.156 

-0.037 

0.036 

Figure  3  shows  three-dimensional  plots  of  the  solution  after  2048  iterations.  The 
solution  exhibits  overshooting  at  the  transition  from  host  material  to  flaw.  That  is,  near 
the  base  of  the  flaw  “tower”  some  values  drop  to  -0.350  and  at  the  edges  of  the  top  of 
the  “tower”  some  values  reach  1.364. 

Figures  4  through  9  are  similar  plots  for  flaw2,  flaw3  and  fiaw4.  Likewise,  Tables  2 
through  4  show  the  convergence  measure  and  upper  and  lower  limits  of  these  solutions. 


Table  2 


flaw2 

Layer  1 

Layer  2 

Layer  3 

Layer  4 

Iterations 

Si 

ESS 

max 

U'MM 

max 

max 

RESS 

max 

100 

0.00972 

EjjKjMJ 

0.154 

iJM 

. 

|iK£| 

500 

0.00166 

ifwil 

f 

1000 

0.00079 

-0.170 

E El 

1.258 

2048 

0.00036 

-0.155 

0.133 

-0.136 

1.232 

BBS 

-0.182 

EB9 

Table  3 


flaw3 

Layer  1 

Layer  2 

Layer  3 

Layer  4 

Iterations 

Si 

ESS 

max 

ESS 

max 

min 

max 

max 

100 

0.01762 

Bawii 

fcjXjWjl 

mmM 

j*  jljJ 

RTfE* 

mmn 

500 

0.00431 

Bjiij;  3 

R 

R 

ESS 

)S*i£3 

1000 

0.00199 

R 

0.041 

HB' 

1.042 

-0.176 

2048 

0.00088 

EXEm 

0.046 

-0.184 

BiB!: 

1.049 

SBI 

Table  4 


|  flaw4 

Layer  1 

Layer  2 

Layer  3 

Layer  4 

Iterations 

sx 

ESS 

max 

min 

max 

max 

ESS 

max 

100 

0.03249 

BjIjWjg! 

IjMTJ 

-0.219 

mmm 

0.470 

-0.199 

500 

0.00600 

BSSl 

E 

sm 

322 

B|j|2 

0.512 

Km 

0.00268 

3 

0.031 

BBa 

J 

j 

0.452 

-0.242 

1.151 

K5M 

0.00123 

EEa 

0.043 

-0.184 

0M?1 

ESQ 

0.369 

1.181 
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The  results  of  this  first  set  of  tests  are  encouraging.  Flaw3  and  fiaw4  are  more  difficult 
to  isolate  and  this  is  to  be  expected.  Figure  10  is  a  plot  of  the  convergence  measure,  <$,, 
for  the  first  500  iterations.  Notice  that  the  deeper  the  flaw  is,  the  slower  the  convergence 
is.  The  "skin  effect”  phenomena  suggests  that  the  results  should  be  better  if  some  of  the 
data  were  taken  at  higher  frequencies.  With  this  in  mind,  more  tests  were  performed. 

(e)  Test  Set  #2 

The  error  in  our  solutions  reported  in  the  previous  section  was  greatest  in  the  layers 
immediately  adjacent  to  the  layer  that  contained  the  flaw.  The  worst  case  for  layer 
discrimination,  then,  would  seem  to  be  having  a  flaw  in  layers  2  and  4.  The  error  that 
would  appear  in  layer  3  may  be  large  enough  to  lead  us  to  believe  that  there  is  a  flaw  in 
that  layer  also.  We  will  call  this  flaw  arrangement  flaw24.  First  we  tried  to  reconstruct 
flaw24  using  data  taken  at  the  same  five  frequencies  as  above  (i.e,  at  1,  3,  5,  7  and  9 
MHz).  The  results  are  shown  in  Figures  11  and  12  and  in  Table  5. 


Table  5 


flaw24 

Layer  1 

Layer  2 

Layer  3 

Layer  4 

Iterations 

Si 

max 

min 

max 

max 

jlllJTfflJljll 

max 

0.00992 

EMU 

mvm 

1.096 

mm 

-0.166 

0.639 

1 

0.00251 

IPPpl 

KPI1 

R 

1.097 

gggl 

E' 

-0.245 

0.930 

0.00120 

iSfSI 

RSI 

1.110 

EE3 

ji 

-0.175 

0.998 

0.00053 

-0.128 

0.128 

-0.123 

1.118 

-0.402 

0.489 

-0.186 

1.045 

At  9  MHz,  the  skin  depth  of  the  workpiece  is  0.066  inches  or  2.4  times  the  layer 
thickness.  A  second  reconstruction  was  performed  for  fiaw24  using  data  simulated  over 
a  wider  frequency  range.  The  frequedes  used  were  1,  5,  10,  20  and  30  MHz.  At  30  MHz, 
the  skin  depth  is  0.036  inches  or  1.3  times  the  layer  thickness.  The  results  of  this  second 
reconstruction  are  displayed  in  Figures  13  and  14  and  in  Table  6  below.  As  expected, 
the  higher  frequency  data  improves  the  solution. 


Table  6 


flaw24 

Layer  1 

Layer  2 

Layer  3 

Layer  4 

Iterations 

Si 

ESI 

max 

W2EM 

max 

HE2EM 

max 

max  | 

100 

0.02541 

Br -‘ul 

1.381 

500 

0.00452 

R  1 

1.148 

1000 

0.00189 

-0.151 

1.239 

1.315 

2048 

0.00072 

BESl 

0.024 

-0.238 

1.242 

DESS 

1.398 

Earlier  we  noted  that  5  sets  of  emf  data  made  the  system  overdetermined  by  a  factor 
of  2.5.  To  see  if  additional  data  sets  from  the  same  frequency  range  would  improve  the 
solution,  we  tried  a  third  reconstruction  using  9  sets  of  emf  data  (1,  2,  3,  4,  5,  6,  7,  8  and 
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9  MHz).  This  increased  the  factor  of  overdetermination  to  4.5.  The  results  are  shown 
in  Figures  15  and  16  and  in  Table  7.  Note  that  the  results  in  Table  7  are  essentially 
identical  to  those  in  Table  5.  The  additional  data  sets  may  be  important,  though,  when 
data  with  noise  is  used. 
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flaw24 

Layer  1 

Layer  2 

Layer  3 

Layer  4 

Iterations 

Si 

El 

max 

min 

max 

ESI 

max 

M23M 

max 

Iggg^sjB 

0.00953 

toyl 

Dji] 

BjTCTjl 

mm 

1 

0.00244 

Bpij 

m 

11  n 

E ESI 

0.00117 

EXE! 

ECU 

H 

on  5 

1.000 

MEEm 

0.00052 

HE01 

0.133 

lERl 

1.120 

-0.399 

■851 

1.043 

In  Chapter  II  a  method  of  constraining  the  solution  was  considered.  A  final  recon¬ 
struction  was  performed  with  all  variables  constrained  to  the  interval  [-0.0001,1.0].  The 
results  are  shown  in  Figures  17  and  18  and  in  Table  8  below. 


Table  8 


flaw24 

Layer  1 

Layer  2 

Layer  3 

Layer  4 

Iterations 

Si 

max 

min 

max 

B29 

max 

WZSM 

max 

flRF 

1.000 

-0.001 

1.000 

9 

isa 

9 

1.000 

-0.001 

1.000 

■ft 

1.000 

R  9 

1.000 

-0.001 

1.000 

13fl89 

DESK 

1.000 

EESI 

1.000 

-0.001 

1.000 

These  results  are  not  very  good  compared  to  the  unconstrained  solutions  we  have  gen¬ 
erated.  More  work  is  needed  in  this  area  to  determine  whether  our  implementation 
of  constraints  should  be  modified.  The  results  are  interesting  in  that  the  convergence 
measure  Si  goes  from  1.0  to  0.104  in  about  100  iterations  and  then  changes  very  lit¬ 
tle  through  iteration  1460.  Then,  at  some  point  between  iteration  1460  and  iteration 
1470,  a  direction  for  significant  improvement  is  found  and  the  solution  improves  steadily 
through  iteration  1560  and  then  again  little  changes  until  some  point  between  iteration 
1960  and  iteration  1970  when  another  “good  direction”  is  found.  Figure  19  is  a  plot  of 
Si  illustrating  this  phenomenon.  More  work  is  needed  to  determine  if  convergence  can 
be  speeded  up  for  the  constrained  cue. 

To  make  sure  that  there  is  nothing  fundamentally  incorrect  with  the  constrained  CG 
algorithm,  we  used  the  actual  flaw  as  the  initial  guess.  This  guess  would  yield  $o=0 
and  terminate  the  algorithm  if  we  were  using  “exact”  data.  Since  we  are  using  “direct 
model”  data,  though,  our  initial  guess  is  not  the  best  fit.  Figures  20  and  21  and  Table 
9  show  the  results  of  this  last  run. 
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Table  0 


|  fiaw24 

Layer  1 

Layer  2 

Layer  3 

Layer  4 

Iterations 

* 

max 

min 

max 

max 

min 

max 

100 

0.01504 

KfjTj 

1.000 

1.000 

500 

0.01256 

lonTii 

1.000 

Ipfff « 

(339 

jfl 

1.000 

WEM 

0.01254 

0.110 

EBi 

1.000 

mllll  J 

0.704 

IftMi  ■ 

1.000 

■a 

0.01253 

BEa 

0.110 

-0.001 

1.000 

EE3 

0.706 

BBS 

1.000 

Notice  that  although  some  isolated  points  build  up  to  help  match  the  data,  the  flaw 
locations  from  the  initial  guess  are  unaltered.  Also,  the  convergence  measure  is  much 
smaller  in  Table  9  than  than  in  Table  8.  Hopefully  this  means  that  our  previous  test 
(Figures  17  and  18)  will  eventually  evolve  (albeit  slowly)  to  this  same  solution. 

Post-processing  Using  Classification  Theory 

Returning  to  the  unconstrained  results,  post-processing  based  on  some  classification 
theory  results  seems  to  “dean  up”  the  solutions.  Using  some  of  the  solutions  as  a 
“training  set”,  0.23  was  determined  to  be  the  optimum  value  for  partitioning  the  data 
into  two  classes,  host  material  and  flaw  (see  Chapter  VI).  If  we  filter  our  first  test  from 
Test  Set  #2  by  assigning  0  to  all  variables  less  than  0.23,  we  obtain  the  solution  shown 
in  Figure  22. 
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3.  Some  Preliminary  Results  Using  Measured  Data 

Several  tests  that  use  actual  measured  data  have  been  performed  to  validate  the 
model.  Results  of  a  few  of  those  tests  will  be  presented  in  the  following. 

A  0.11  inches-thick  work-piece  of  graphite  epoxy,  in  the  form  of  woven  fiber,  has 
been  used  for  measurement  in  the  laboratory.  Holes  at  various  depths  have  been  drilled 
into  this  workpiece,  as  shown  in  Figures  23-25.  A  multi-frequency  measurement  of  the 
work-piece  has  been  performed,  as  explained  in  Chapter  V.  For  the  test,  an  algorithm 
based  on  a  ten-layer,  five-frequency  model  has  been  chosen.  Figures  26-30  show  the  plots 
of  the  real  and  imaginary  parts  of  the  measured  data,  B,  for  the  frequencies  of  2,  4,  8, 
12  and  16  MHz. 

In  order  to  reconstruct  the  flaws  from  the  measured  data,  a  normalization  process 
is  performed.  First,  the  left  hand  side,  Y,  of  equation  (II  —  1)  is  calculated  by  using 
the  flaw  as  the  unknown  X\  these  results  are  shown  in  Figures  31-35.  Next,  the  read 
and  imaginary  parts  of  the  measured  data,  B,  are  normalized  with  respect  to  Y  at  each 
frequency.  This  normalization  procedure  attempts  to  compensate  for  amplitude  and 
phase  variations  during  the  data  acquisition.  The  sup-norm  of  the  real  and  imaginary 
parts  of  By  as  well  as  Y,  versus  frequency  sue  depicted  in  Figure  36.  It  is  noted  that  the 
real  and  imaginary  sup-norms  of  B  sue  greater  than  the  real  and  imaginary  sup-norms 
of  Y  at  all  frequencies.  The  normalized  value  of  B  is  then  used  to  solve  for  the  unknown 
X  of  (II  -  1). 

Figures  37-39  are  samples  of  the  real  and  imaginary  parts  of  the  64  x  64  elements 
of  the  operator,  A,  that  appears  in  (II  —  1),  for  different  layers  and  frequencies.  The 
elements  of  the  operator  are  highly  localized  around  the  origin,  as  indicated  in  Figures 
37-39.  This  simply  states  that  the  coupling  effect  is  localized  in  physical  space. 

Figures  40-42  show  the  results  of  the  reconstruction  of  the  flaw.  The  constraint 
of  positive  conductivity  is  used;  the  Levenberg-Marquardt  parameter  is  set  to  1.5,  the 
number  of  iterations  is  500  and  the  threshhold  is  0.3.  Figure  43  shows  the  flaw  and 
its  reconstruction  in  gray  scale.  The  flaws  of  layers  two  to  six  have  been  somewhat 
reconstructed.  Excessive  noise  is  accumulated  in  layers  one,  two,  seven  and  ten,  which 
makes  the  flaws  somewhat  unrealisable. 

Although  the  reconstruction  with  actual  data  is  not  perfect  at  this  preliminary  stage, 
it  is  very  promising  for  further  investigations.  This  shortcoming  is  in  part  due  to  the  fact 
that  the  theoretical  model  for  this  study  is  based  on  the  Born  approximation.  For  the 
future  study  a  more  accurate  and  realistic  model  for  reconstruction  is  proposed,  which 
is  discussed  in  the  following  chapter. 
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Figure  24:  Ten  layered  graphite  epoxy  with  some  drill  hole  flaws  in  various 
depths;  layer  5  -  layer  8. 
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Figure  25:  Ten  layered  graphite  epoxy  with  some  drill  hole  flaws  in  various 
depths;  layer  9  -  layer  10. 


Figure  27:  Real  and  Imaginary  parts  of  the  measured  data,  B ,  at  4  MHz. 
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Figure  32:  Real  and  Imaginary  parts  of  the  measured  data,  Y ,  at  4  MHz. 
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Figure  36:  The  sup-norms  of  the  read  and  Imaginary  parts  of  B  and  Y. 
-hi^are  real  and  imaginary  sup-norms  of  B ;*,o,  are  real  and  imaginary 
sup-norms  of  Y  respectively. 


Figure  43:  Gray  scale  of  the  flaw  and  its  reconstruction;  first  two  rows  is 
the  actual  flaw,  and  second  two  rows  are  the  reconstruction. 
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CHAPTER  IV 

RIGOROUS  FORMULATION  OF  THE  INVERSE  PROBLEM 


1.  Introduction 

Throughout  this  report  we  have  simplified  the  model  by  approximating  the  anomalous 
current  density,  •  E,  as  =  <r^  •  Eo,  where  Eo  is  the  incident  electric  field 

produced  by  the  whip  (or  any  other  source)  in  the  absence  of  the  flaw.  This  approximation 
linearizes  the  problem,  and  obviously  simplifies  the  task  of  inversion.  Nevertheless,  the  true 
problem  is  nonlinear,  because  it  involves  the  product  of  two  unknowns,  and  E. 

In  our  earlier  work  in  eddy-current  inversion  we  solved  the  nonlinear  problem  by  re¬ 
sorting  to  a  package  of  nonlinear  equation-solvers,  called  MINPACK.  In  addition,  we  used  a 
variable-metric,  nonlinear-programming  package,  called  VMCON,  in  an  attempt  to  satisfy 
constraints.  These  worked  well  for  rather  small  problems,  but  for  the  three-dimensional 
problems  that  we  are  currently  addressing,  it  becomes  necessary  to  use  storage-efficient 
algorithms  that  allow  us  to  take  advantage  of  special  structures  of  our  matrices,  such  as 
convolution.  Hence,  we  will  go  back  to  the  rigorous  problem  to  see  how  it  might  be  attacked 
in  the  fight  of  what  we  now  know  about  iterative  algorithms  for  solving  large  problems. 
There  should  be  no  doubt,  however,  that  the  rigorous  solution  of  a  large  nonlinear  problem 
will  be  time-consuming. 

2.  The  Nonlinear  Coupled  Integral  Equations 

Let  region  1,  where  the  excitation  source  and  sensors  are  located,  be  above  the  slab, 
and  region  2  be  the  6lab,  which  contains  the  anomaly.  Assume  that  we  measure  the 
magnetic  induction  field,  B,  from  which  we  compute  EMF’s  in  the  usual  way.  Then  the 
appropriate  pair  of  coupled  integral  equations  is 

E(r)  =E<°>(r)  +  f  Cg'Vlr')  ■ 

J  flam 

=E<°>(r)+  f  G^e)(r|r')-5{0)(r')E(r')dr'  (1) 

J  flam 

B(r)  -  B"»(r)  /  G(”'>(r|r') .  J«(r')  •  E(r')*'. 

J  flam 

The  superscripts  on  the  Green’s  functions  denote  their  type,  the  first  denoting  the  type  of 
field  (electric  or  magnetic),  and  the  second  the  type  of  current  source  (electric  or  magnetic), 
whereas  the  subscripts  denote  the  regions  which  are  coupled  by  the  Green’s  function;  the 
first  subscript  denotes  the  region  which  contains  the  field  point,  and  the  second  the  region 
that  contains  the  source  point.  The  superscript,  (0),  on  E  or  B  denotes  the  incident  field, 
which  is  the  field  that  exists  in  the  presence  of  the  slab  without  the  flaw.  The  left-hand 
side  of  the  second  equation  i6  integrated  over  the  sensor  roil  to  give  the  measured  EMF. 

This  system  is  nonlinear  (or,  more  precisely,  bilinear)  because  of  the  presence  of  the 
product  ff^(r')  •  E(r').  We  have  been  linearizing  the  problem  by  replacing  E(r')  in  the 
second  equation  by  E^°^(r'),  and  then  ignoring  the  first  equation. 


To  solve  this  system,  we  first  discretize  it  in  the  usual  way  by  means  of  the  method  of 
moments,  and  then  apply  an  iterative  technique  to  the  resulting  algebraic  equations.  The 
iterations  start  by  replacing  E  in  the  second  equation  by  and  then  solving  for 
by  using  one  of  the  iterative  methods  just  described.  This  step  uses  measured  data,  and 
is  the  ‘inverse’  phase  of  the  problem. 

Once  we  have  an  acceptable  approximation  to  d^a\  we  substitute  it  into  the  first 
of  (1)  and  solve  the  resulting  ‘direct  problem’  for  an  improved  version  of  E,  the  electric 
field  within  the  flawed  region.  This  step  will  usually  converge  more  rapidly  than  the 
inverse  phase  because  the  equations  are  better  conditioned.  We  can  use  the  same  iterative 
technique  here  that  was  used  during  the  inverse  phase. 

The  result  of  the  direct  phase  is  then  substituted  into  the  second  equation,  and  the 
second-level  inverse  problem  is  solved.  The  process  is  continued  until  the  error  in  the 
solution  is  of  the  order  of  the  error  in  the  measured  data. 

The  process  that  we  have  just  suggested  may  not  be  the  most  efficient  way  to  solve 
system  (1).  We  will  look  at  conjugate  gradients,  and  see  how  to  handle  the  nonlinearity 
there. 

Keep  in  mind  that  if  we  use  multifrequencies  on  this  problem,  a  number  of  distinct 
incident  fields  and  Green’s  functions  will  have  to  be  computed  and  stored,  one  for  each 
frequency.  Because  of  the  convolutional  nature  of  the  problem,  however,  we  will  be  able 
to  utilize  storage  efficiently. 

3.  Discretization  of  the  Coupled  System:  Method  of  Moments 

We  start  by  writing  the  Green’s  function  as  a  two-dimensional  Fourier  transform, 

G ■<“>(*-*', y  -  y*.  *,*')  ®  (2) 

with  a  similar  expression  for  G^me^.  From  here  on  we  will  drop  all  subscripts  that  are 
related  to  regions  1  and  2. 

When  (2)  is  substituted  into  (1)  the  integrals  over  x'  and  y1  become  Fourier  transforms. 
The  resulting  integral  equations  become: 

E(r)  =  E‘0)(r)  +  ^  jdz' G <ee>(z,  z')  •  J a{z')^k-^k^dkxdky  (3)(a) 
B(r)  -  B<°>(r)  =  Jdz'JJ' “  G <”">(*,  z')  ■  3a(z')e-^k'*+k^dktdky.  (3)(b) 

The  measured  data  are  usually  EMF’s  induced  into  flat  coils  that  are  oriented  parallel 
to  the  workpiece.  In  such  a  case,  therefore,  we  are  only  interested  in  j u/x  the  z-component 
of  (3)(b)  integrated  over  the  sensor-coil  position.  This  produces  another  integral  relation 
with  a  new  Green’s  function  to  replace  (3)(b): 

EMF( r)=^r  (3)(c) 
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where  is  the  2  row-vector  of  G^me^  multiplied  by  the  transfer  function  corresponding 

to  integration  over  the  sensor  coil.  The  Green’s  functions  that  appear  in  (3)  are  computed 
in  Appendix  A. 

The  discretization  of  the  coupled  system  of  integral  equations  is  done  by  subdividing 
the  region  of  space  occupied  by  the  anomaly  into  a  regular  grid  of  ( Nx  +  1)  x  (Ny  +  l)x 
(AT,  +  1)  cells,  each  of  size  6xx6yx  6Z,  and  then  expanding  the  electric  field  and  anomalous 
conductivity  tensor  using  pulse  functions  defined  over  the  grid: 


N.  Nt  N, 


and 


»wm  =  £  £  £ 

1=0  m=0  j=0  x  V  * 


The  pulse  function,  Py(s),  satisfies 


(5) 


=  Xj<*<j  + 1; 

3'  \0,  otherwise. 

Note  that  the  anomalous  current  density,  which  is  given  by  the  product  of  (4)  and  (5),  has 
exactly  the  same  expansion  in  pulse  functions  as  either  (4)  or  (5),  except,  of  course,  that 
the  expansion  coefficients  are  given  by  the  term-by-term  product  &imj  •  E/mj. 

We  are  going  to  use  Galerkin’s  variant  of  the  method  of  moments  to  complete  the 
discretization.  In  Galerkin’s  method,  we  ‘test’  the  integral  equations  (3)  with  the  same 
pulse  functions  that  we  used  to  expand  the  unknowns  in  (4).  For  example,  we  form 
moments  of  (3)(a)  by  multiplying  (3)(a)  by  Pi(x/6X )Pm {,ylSy)Pj{zj6z))l6xSy6t,  and  then 
integrating  over  each  cell.  This  yields  an  algebraic  system  for  E/mj : 


=  E, 


A'.  A'.  ,v, 

-  £  £  £  G$7’(!  fLUj  ■  Elmj, 

L=0  M= 0  7=0 


(7) 


where 


G<*e)(/  —  Z,m  —  M)  =  T-j{kXyky) 

/nin(kxSx/2)y  /sm(ky6v/2)y  dk^ 

'  kZ6X/2  '  '  kySy  f2  ' 
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Tjj  is  the  result  of  the  integrations  with  respect  to  z  and  z': 

=(**)  f*i+)  /*»(+)  __ 


r*V  r*n+)  x 

r„  (*„*,)  =  ^'g<">(2,2') 

r4+>  /*i(+)  _ 


=L/zL 

/*l+)  /■*>(+)  _ 

+L/!L->  dz'GD(z  ~  2° 


(9) 


The  subscripts,  I,  D,  R,  denote  the  incident,  depolarization,  and  reflected  Green’s  functions, 
respectively. 

We  will  show  how  to  compute  the  integrals  of  the  incident  and  reflected  Green’s 
functions  later.  The  depolarization  function,  however,  can  be  integrated  by  inspection; 
the  result  is  j(6kjSx/ivezt) a,az,  where  Skj  =  1,  if  t  =  j\  =  0,  otherwise. 

In  taking  moments  of  (3')  we  must  keep  in  mind  that  the  sensors  lie  in  a  fixed  z-plane 
above  the  workpiece.  Hence,  we  multiply  by  only  Pi(x/Sx)Pm(y/6y)  and  then  integrate  to 
get  the  algebraic  system: 


N.  N„  N, 


EMFim  =  Y,1l1L  g</°(*  -L,m-M)-  aLMJ  ■  ELMJ> 


L= o  M~  o  J= o 


(10) 


where 


G (/e)(I  -  L,m  -  M) 


juftoMy  JJ00e-j[k.6.il-L)+k,S,(m-M))  f(/e)(fcr,fc„) 


(11) 


ant^  =(Ee)  .*;(+)  _ 

Tj  (*.,*,)  =  J  dt'G<B‘\z,  z').  (12) 

-.(£*)  _ 

Keep  in  mind  that  T  ,  G*Be\  and  G*fie)  are  vectors,  not  dyads;  they  take  vectors  into 
scalars  (EMF). 

Equations  (7)  and  (10)  are  the  coupled  system  of  discrete  equations  that  define  the 
rigorous  inverse  problem.  The  system  is  nonlinear  (bilinear)  in  the  unknowns,  cimj,  Ejmj. 
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Note  that  in  two  of  the  three  dimensions,  the  indices  appear  in  a  Toeplitz  (or  convolution) 
form.  We  will  discuss  the  computation  of  the  matrices  shortly. 

4.  Multifrequency  or  Multiview  Reconstruction  Methods 

Eqnation  (10)  indicates  that  if  we  measure  the  EMF’s  at  a  single  z-level,  then  we 
don’t  have  enough  data  to  reconstruct  a  three-dimensional  flaw.  The  additional  data  to 
reconstruct  the  third  dimension  can  be  obtained  several  ways.  For  example,  we  can  excite 
the  system  from  a  single  coil  at  a  number,  N/t  of  frequencies,  where  Nf  >  Nz  +  1.  In  that 
case,  the  Green’s  functions  and  electric  field  in  (7)  and  (10),  as  well  as  the  measured  EMF’s, 
vary  with  frequency  (but  the  unknown  conductivities  are  assumed  to  be  independent  of 
frequency). 

If  we  let  the  integer,  n,  index  the  frequency  parameter,  then  (7)  and  (10)  become 

=*W») 

Nm  N,  Nt 

-  53  X/  ®iJ**(*  -  L,m-  M\n)  •  &lmj  •  EiMr(n)  (13)(a) 

L=0  M—0  7=0 

Nm  N,  Nm 

EMFtm(n)  =  ^2  X*  -  L,m  -  Af;n)  ■  <tlmj  ■  E lmj(h),  (13)(b) 

L= 0  M= 0  7=0 


where  G is  a  vector,  not  a  dyad. 

This  is  the  “multifrequency  reconstruction  algorithm”.  Note  that  here  the  Green’s 
functions  vary  with  frequency,  as  do  the  measured  EMF’s  and  the  computed  electric  fields. 

Another  method  is  to  excite  the  system  at  a  single  frequency,  but  with  a  number  of 
different  exciting  coils.  In  this  case,  only  the  electric  fields  and  measured  EMF’s  vary  with 
n,  which  now  indexes  the  location  of  the  exciting  coil.  Hence,  the  system  of  equations 
becomes 


N.  N,  AT, 

Eimj(n)  =E imj(n)  ~  5Z  5Z  ~  L,m  -  M)  •  clmj  •  EiAfj(n)(14)(a) 

L=  0  M= 0  7=0 
Nm  N,  Nt 

BMFlm(n)  =  '£'£.  Y.  G</e)(<  -  L,m  -  M)  ■  SLMJ  ■  E1M,(n),  (14)(b) 

L= 0  M= 0  7=0 

where  G  is  a  vector,  not  a  dyad. 

This  is  the  “multiview  reconstruction  algorithm”.  It  has  the  advantage  over  the 
multifrequency  algorithm  in  not  requiring  the  computation  and  storing  of  a  number  of 
different  Green’s  functions.  This  is  not  too  important  when  the  linearized  theory  is  used 
for  inversion,  but  could  be  quite  significant  in  the  present  context  of  the  rigorous  theory. 
Another  advantage  to  the  multiview  algorithm  is  that  it  can  be  applied  to  problems  where 
the  unknown  conductivity  is  frequency  dependent. 
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An  advantage  of  the  multifrequency  approach,  however,  i6  its  relative  flexibility  in 
the  laboratory,  especially  with  regard  to  hardware  implementation.  That  is,  it  is  probably 
easier  to  change  frequencies  than  to  change  excitation  sources. 

5.  The  Conjugate  Gradient  Algorithm  for  the  Coupled  System 

We  are  going  to  apply  the  conjugate  algorithm  to  the  coupled  system  of  either  (13) 
or  (14).  This  will  be  an  interesting  application  of  conjugate  gradients  to  a  nonlinear 
problem,  and  in  this  context  we  will  use  some  ideas  of  Stephen  J.  Norton,  “Iterative 
Seismic  Inversion”,  Geophys.  J.  Roy.  Astr.  Soc.  (submitted). 

We  will  consider  (14)  first,  and  work  in  the  coordinate  system  in  which  &  is  diagonal. 
Then  we  can  rewrite  (14)(b)  as 

N.  Nt  N, 

0  =  Rimn  =  - EMFim{n )  +  -  L,m  -  M)olm jxElm Jx{n)+ 

L= 0  M= 0  J=  0  ■ 

G%\l  -  L,m-  M)aLMjyELMjy(n)+ 

-  L,m-  M)<TLMJzELMJz{n ) 

N. 

=  ~EMFim{n)  +  ^2  0  (rjzEjx{n )  +  G^e)  0  <rjyEjy(n)  +  0  <rjzEjz(n)  , 

J=o 

(15) 

where  Rimn  is  the  Zmnth  component  of  the  residual  vector,  and  0  denotes  a  two-dimensional 
discrete  convolution. 

We  form  the  squared-norm  of  the  residuals 

N.  N,  Nv 

&(<TLMJx,<rLMJy,0LMJz)  =  yi  I!]  Rimn Rlmn,  (16) 

1=0  m=0  n=0 

where  Nv  is  the  number  of  “views”  (i.e.,  the  number  of  source  locations),  and  then  differ¬ 
entiate  with  respect  to  &lmj  =  {°LMJx><rLMjy,CLMJZ}- 


d * 

dlTLMJ 


Nm  Nt  AT, 


— 2Re  ^2  ^2  ^tmn  S  ^  (/  —  L,m  —  j(n)  \  -f 

1=0  m=0  n=0  '  ' 


LMj{n ) 
d&LMJ 


—IGradiMJ- 


In  (16)  and  (17)  the  asterisk  denotes  complex  conjugation.  We  are  using  dyadic  notation 
in  (17).  For  example,  the  first  term  within  the  square  brackets  is  the  diagonal  part  of  the 
dyadic  product  of  the  vectors  G;Ee'‘  and  E£M  j(n).  We  can  write  (17)  in  the  form  of  a 
two-dimensional  convolution,  also,  but  we  will  not  need  it  here. 
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E-lm  j(n)  satisfies  the  conjugated  version  of  (14)(a)  (but  keep  in  mind  that  d  is  real); 
what  can  we  say  about  the  dyad  dE*LMJ{n)fd&iMjl  This  can  probably  most  easily  be 
computed  by  using  a  finite-difference  approximation,  in  which  (14)(a)  is  solved  for  two 
values  of  a,  and  then  the  difference-quotient  is  formed.  It  is  possible,  however,  that  some 
insight  into  the  computation  of  the  derivatives  can  be  obtained  by  directly  differentiating 
(14)(a)  with  respect  to  keeping  in  mind  that  a  is  diagonal.  The  result,  after  a 

straightforward  computation,  is 


-  l',m  -  m1)  ■  Diag(E<<m«j/(n)) 


flEjmj(n) 

d&t'm'ji 


N.  Nt  Nm 

-EEE  GS"’(i  iLUJ  ■ 

L= 0  M=0  J=  0 


dE  LMjjn) 

d&l'm'j' 


(18) 


This  is  the  tensor  version  of  (14)(a).  After  solving  for  5E we  keep  the 
“diagonal”  entries  in  which  =  (/',  m',  jf').  These  correspond  to  local  gradients,  in 

which  the  variations  of  conductivity  and  electric  field  are  within  the  same  cell. 

Actually  solving  (18)  is  out  of  the  question,  because  of  the  computations  involved. 
A  possibility  is  to  keep  only  the  local  gradients  in  (18),  and  solve  the  resulting  equation. 
We  might  wish  to  simplify  the  problem  further  and  approximate  the  solution  for  the  local 
gradients  by  the  term  on  the  left-hand  side  of  (18)  when  (l,rn^j)  =  ( V,m',j '): 


U/mj(n)  = 


Oft’,  (0,0)  £,„„(») 

.G(">,(0,0  )£,„,. (») 


G<;;>„(0,0)£,m„(n) 

G&>(0,0  )E,mi,(n) 
<0'  0,0)£,„„(n) 


^>>**(0, 0)-E)  m,'*(n) 

G&Uo.O  )£,„„(») 

G<J">(0,0)£,„>I(n). 


(19) 


This  would  be  a  form  of  Born  approximation  for  the  local  gradients. 

The  conjugate  gradient  algorithm  starts  with  the  iterative  step  + 

dkf^LMJ-  We  choose  a*  to  minimize  +  at/^j)  with  respect  to  at,  for  a  given 

/Imj-  We  will  shortly  determine  the  optimum 

Let  us  first  derive  an  expression  for  the  Jfeth  residual,  (15),  using  a  linearized  expression, 
in  which  the  electric  field  is  replaced  by  its  ( k  —  l)th  approximation: 


•^imn  =  “  EMF[m(n)  +  ^2  ^  ^  G^fe\l  -  Ltm-  M)  •  (5^ M lj  + 

L= o  M- o  J=o 

+«»EEE  G</°('  -*.—*)■  Am;  •  4‘m‘A«) 

i=0W=0J=0 

(20) 
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Upon  substituting  this  result  into  (16),  we  get 


N.  N,  Nv 


Nmn  + 


)  (*%£’  + 


)' 


(21) 


nm  tv*  , 

*“’=EEE(* 

1=0  m=0  n=0  ^ 

-(-  ofc2Re[F^  •  #*-*>]  +  aJ||F<*>||*, 

where  we  are  using  vector-matrix  inner-product  notation.  Hence,  is  minimized  when 


a* 


Re[fl(*-i).F<*>] 


(22) 


and  when  this  is  substituted  into  (21)  we  find  the  minimum  value  to  be 


*(»>  _  *€.-«  _  1MF<>1  •  «(t-11)]8 

||fW||i 


(23) 


The  greatest  decrease  occurs  when  F ^  =  R^k~^  in  (23).  What  this  means  in  terms 

/(*) 

lmj>  the  change  in  &lmj  can  be  determined  by  returning  to  (20): 


Re 


.  F(k) 
Nm  N, 


-»•£  E  XX»’[E  E 

IsO  m=0n=D  ^1=0  M=0  J=0 

= E  E  E tlL,  [r'EEE -  s.«  -  uvti&w) 

L=0  M=0  J=0  |=0m=0n=0  '  / 


(24) 


Upon  comparing  the  summation  term  within  the  large  square  brackets  of  (24)  with  the 
expression  for  the  gradient,  (17),  we  see  that  the  summation  term  is  the  linearized  gradient, 
which  is  (17)  with  j(n)/ dB lm  j  =  0.  Hence,  (24)  shows  that  the  maximum  decrease 

in  the  norm  of  the  residuals  occurs  when 


This  is  the  steepest-descent  direction. 

There  is  an  important  orthogonality  relation  that  holds: 

N.  N.  ...  _ 

5^  *  Gra<^LMJ  ~ 

L= 0  M=0  7=0 
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In  order  to  derive  this  result,  substitute  the  final  version  of  (20)  into  the  expression,  (24), 
for  the  linearized  gradient  and  get  the  following  recursion  relation: 


N.  N,  Nv 


=R'E  E  E  k~n  ’  +  Diag  (o</*>*(l  -  L,m  -  M) *&*/(»)) 

1=0  m=0  n=0  '  / 

=5r -  £.m  -  !/)£&”» 

/=  0  m=0  n=0  ' 


When  the  definition  of  F$n  from  (20)  is  substituted  into  (27),  and  the  resulting 
expression  is  multiplied  by  and  then  summed,  we  get 

EEE  f£,j  ■  VrJ&j  =EEEC-  SMfcy  +  afcRe||F(fc)||2 

L=0  M=0  J=0  I<=0  M=0  J=0  (2g) 

=Re[R(fc"1)  •  F<fc>]  -  Re^*'1)  •  F(*>] 


We  have  used  (22)  and  (24)  in  arriving  at  the  final  result. 

The  conjugate  gradient  algorithm  starts  with  a  pure  gradient  step,  (25),  for  Jfc  =  1, 
and  then  continues  with 

flMJ  =  Wi1;  +  hfiMJ  •  (29) 

will  be  chosen  to  minimize  the  denominator  in  (23),  thereby  guaranteeing  an  improved 
convergence  rate. 

Before  deriving  an  expression  for  6* ,  we  will  derive  another  expression  for  the  numer¬ 
ator  of  at,  which  appears  in  (22). We  start  with  (24),  and  use  (29): 

Numerator  of  a*  =Re[J^*_1^  •  F^] 

=EEE  ■  JMiii’  (so) 

1=0  M= 0  J-0 

=  0ra<tk-1)\\2, 

where  the  final  result  is  due  to  the  orthogonality  relation,  (26). 

The  derivation  of  bt  starts  with 

iif^f 

l=0m=0n=0 

N.  N ,  N.  N.  N,  N.  V  } 

fk)  fk) 

L= 0  M=0  J-0  L'=0  M'—O  J'= 0 


where 


r>(*) 

1  LMJ.L'M'J' 


=  E  E  El>iH(G(/’>(l  —  i, m  —  M) lEgji'W)- 

i=0  m=0  n=0 

Diag(G,/')'(i  M’ )e£*&(")) 


(32) 


is  Hermitian  in  (LMJ;  L'M'J1).  Equations  (31)  and  (32)  are  derived  from  the  definition  of 
■^imn  i®  (20),  together  with  an  interchange  of  summations,  and  a  rearrangement  of  terms. 

Upon  substituting  (29)  into  (31),  and  then  minimizing  the  result  with  respect  to  fc*, 
we  get  for  the  optimum  value  of  bk 


_ile  2^i= 0  ’  •  ‘  2^J'=0  1  LMJ\L‘M[ 


Nm 


,(fc) 


bk  = 


tGradf  ^ 


lLMJ 


J  L'M'J' 


]Cl=0  ‘  ‘ '  £j'*=0  ^LMJ-.L'M'J'f  LMJ  ’  f  L'M'J 


yN.  r(fc) 


7(k-l)  -(k-1) 


(33) 


It  is  straightforward  to  rewrite  (27)  in  terms  of  T 


(*> 

LMJ.L'M'J" 


Gra^LMJ  =  Gra^LMJ  +  a*Re  E  E  E  fvM'J'? 


(*)• 

LMJ.L'M'J" 


L'=0  M'=0  J'=0 


(34) 


There  are  two  fundamental  conjugate  gradient  orthogonality  relations: 

E  "  E  ^IMJ  •  ?L'M' J'Rei^LM J-, L'M'J'}  ~  (35)(a) 

t=o  J' =o 

E  E  E  J  =  0.  (35)(b) 

i=o  M= o  J=o 


The  first  relation  is  easily  proved  by  substituting  bk  from  (33)  into  (34),  multiplying 
by  Re{ri‘i, J;L, w }, }fv m-j, ,  and  then  summing  over  (LMJ)  and  (L'M'J1). 

In  order  to  prove  the  second  relation,  we  multiply  (34)  by  Gra<fLMj  and  sum 


E E E^m1;  =EEE <M5i# •  g-4mV 

L=sO  M= 0  J=0  L=0  M=0  ■7=0 

-f  a*Re  2^  •  •  *  2^  Gra<rLMJ  ■  J l'M'J'1  lmj-.l' 
L=0  J'~  0 


(36)(a) 
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Next  substitute  (Srad^y  =  f'lhj  —  bkj^LMJy  ^hich  follows  from  (29),  into  the  second 
term  on  the  right-hand  side  of  (36)(a),  and  rewrite  the  right-hand  side  as 


IJC^rad^  1}||2  +  a*  Re  ^  •••  ’Y]  J'fhMJ  '  fuM'. 

L- o  J'=o 


(36)(b) 


!<*) 


The  second  term  under  the  summation  sign  vanishes  because  of  the  first  orthogonality 
relation,  (35)(a),  and  the  first  term  under  the  summation  sign  is  equal  to  ||^fc^||2,  because 

of  (31).  Therefore,  when  we  recall  that  at  =  —  ||(?rad(k~D |j2/||^fc^||2,  which  follows  from 
(22)  and  (30),  we  see  that  (35)(b)  vanishes,  which  proves  the  theorem. 

Our  last  task  is  to  derive  a  simpler  expression  for  fcf  In  order  to  do  this  we  must 
assume  that  the  change  in  the  model  (i.e.,  conductivity)  is  sufficiently  small  to  permit  us 
to  say  that  the  electric  field  within  the  anomalous  region  does  not  not  change  too  much 
from  iteration  to  iteration.  If  this  is  the  case,  we  can  write  ss  r^-1),  so  that 


— . . .  Vw*  r(k-1)  sO-1) 

,  _  m2-jL= o  Z^J'gQ1  LMJiL'M' J'uraaLMJ  '  Il'M'J ' 

k  ||F(fc-1>||2 


(37) 


where  we  have  used  (31)  and  (33). 

Next,  take  the  complex  conjugate  of  (34),  replace  a*  in  the  resulting  expression  by 
(22)  and  (30),  replace  k  by  k  —  1,  and  get,  after  rearranging: 


Nm 

«*£  E 


Ny  N. 


E# 


-*>  r(fc-D 
M'J'l 


LMJ-.l'M'J' 


L'=0  M'=0  J'= 0 


Upon  taking  the  dot  product  of  (38)  with  ^radf^^J ,  and  then  summing  over  LMJ,  we 
get 


RcV'  f  r(M)  II***  1}ll2  u^r ,  .Kfc- 

K  2-/  *  ‘  ‘  2-/  1  LMJ-,L'M'J'JL'M'J‘  •  &ro“LM  J  ~  7=  — — — ||G^ rad> 

||Graat*-2lj|2 


Dll2 


(39) 


L=  0  J'=  0 


where  we  have  used  (35)(b).  Upon  comparing  (39)  with  (37),  the  result  follows  immedi¬ 
ately: 

*  ||5ro<J(‘-J)||J ' 

Note  that  a*  and  fe*  are  negative  here,  whereas  they  were  positive  in  our  previous 
notes  on  the  conjugate  gradient  algorithm.  This  is  due  to  the  fact  that  we  have  defined 
our  residual  vector  to  be  the  negative  of  the  previous  definition. 
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We  can  summarize  these  results  by  saying  that  the  nonlinear  conjugate  algorithm 
is  similar  to  the  linear  one,  except  that  the  operators  must  be  updated  at  each  iteration, 
because  the  electric  field  is  updated.  In  addition,  we  may  have  to  ensure  that  the  increments 
in  the  solution  vector  are  not  too  large. 

Let  V  be  the  array  of  complex  scalars,  {EM Fim(n)},  where  Imn  index  the  array 
elements,  and  define  the  operator 

N.  N,  N, 

»W(/)=EEEG</' 

1=0  M=0  7=0 

Nm  N,  N. 

Dia* 

L=0M=0  7=0 

together  with  its  adjoint 

A'(l>(X)  =  R«E  E  E  ( G '/•>•(/  -  L,m  -  M) ) .  (42) 

1=0  m=0  n=0  '  ' 

Note  that  A  produces  an  array  of  complex  scalars  from  an  array  of  real  vectors, 
whereas  A *  produces  an  array  of  real  vectors  from  an  array  of  complex  scalars.  The 
electric  field  at  the  Jfeth  step  is  computed  from  (I4)(a),  using  the  fcth  approximation  of 
&LMJ- 

We  can  derive  a  recursion  relation  for  the  electric  field  when  j  =  iL  j 

in  the  following  way.  Substitute  this  expression  into  (14)(a),  in  which  the  electric  field  is 
replaced  by  its  A:th  approximation.  Subtract  from  this  result  (14)(a),  in  which  the  electric 
field  and  conductivity  tensor  are  replaced  by  their  ( k  —  l)th  approximations.  Then,  if  we 
say  that  E^f  j  K  >  we  6e*  the  final  result 

Eimi  =  Elmj  +  a*  5Z  G77 e)(l  —  L,m  —  M)  •  fLMJ  ■  (43) 

1=0  M=0  7=0 


The  second  term  on  the  right  hand  side  of  (43)  is  the  correction  vector  induced  in 
the  electric  field,  due  to  the  correction  term  in  the  conductivity.  Clearly,  (43)  gives 
explicitly  in  terms  of  Hence,  this  is  an  explicit  recursion  relation  for  the  electric 

field  vector.  We  can  get  an  implicit  relation  by  simply  replacing  the  electric  field  vector 
in  the  correction  term  of  (43)  by  E^,  which  means  that  we  would  have  a  large  system  to 
solve  in  (43).  Therefore,  the  explicit  recursion  relation  has  an  obvious  advantage  over  the 
implicit.  We  should  also  remember  that  (43)  was  derived  under  the  assumption  that  the 
changes  in  conductivity  are  not  large,  from  iteration  to  iteration.  This  means  that  a *  is 
small. 

The  conjugate  gradient  algorithm  starts  with  an  initial  guess,  Xo,  from  which  we 
compute  Rq  =  Y  —  A^Xo,  P\  =  Qo  =  A*^Ro-  In  addition,  we  have  a  convergence 
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parameter,  e.  Then  for  k  =  1,...,  if  Test  =  ||.Rk||/||y||  <  e,  6top;  Xk  is  the  optimal 
solution  of  (14)(b).  Otherwise,  update  Xk  by  the  following  steps: 


Sk  =  A^Pk 


IIQfc-1  II 
llftll 


Xk  =  Xk- 1  +  a.kPk 


Rk  =  Rk- 1  —  akSk 

Qk  =  A'wRk 

,  IIQtll2 
*  IIQ*-i  ||2 
Pk+ 1  =  Qk  +  bkPk- 


(44) 


comment  (1):  For  a  linear  problem,  solved  in  infinite  precision  arithmetic,  the  algorithm 
terminates  at  the  Afth  step  when  Qm+i  =  0,  80  that  Xm+ i  is  the  least-squares  solution 
of  Y  =  AX.  We  have  already  shown  that  the  vectors  Qo,  Qi,  Q 2 ,  -  -  - ,  are  mutually 
orthogonal,  as  are  the  vectors,  Si,  S2,  S3,. . ..  In  addition 


if  j  <  k  , 
otherwise. 


comment  (2):  This  suggests  that  we  monitor  the  iterates  {Qt}  for  loss  of  orthogonality, 
and  restart  when  the  condition  IQ^Qk+i]  >  ci Q^+iQk+i  is  satisfied,  where  ei  =  0.2 
(say).  When  this  occurs,  we  set  6*  =:  0  in  the  last  line,  and  then  continue  (i.e.,  we  restart 
with  a  pure  gradient  step). 

comment  (3):  Allen  McIntosh,  Fitting  Linear  Models:  An  Application  of  Conjugate 
Gradient  Algorithms,  Springer- Verlag,  1982,  gives  an  alternative  expression  for  bk’. 

.  QUQi-Qk-,) 

IIQ1-1IP  ’ 

which  seems  to  produce  Qk’ s  that  are  more  orthogonal,  when  using  the  criterion  of  com¬ 
ment  (2).  This  definition  requires,  however,  that  an  extra  array  to  store  Qk-i  be  made 
available.  This  is  no  problem  if  we  monitor  for  orthogonality  for  the  purpose  of  restarting, 
because  that  array  is  required  anyway. 

The  Toeplitz  operations  that  are  a  part  of  A  and  A *  are  evaluated  by  using  the  FFT, 
as  described  in  Appendix  A  of  Chapter  II.  This,  together  with  the  fact  that  the  storage 
requirements  are  reasonably  modest,  are  the  reasons  why  the  conjugate  gradient  algorithm 
becomes  attractive  for  large  problems  in  our  model. 

APPENDIX  A 

Internal  Green's  Functions,  Gu,  G22,  in  the  Bulk  Model 
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We  are  going  to  sketch  a  procedure  that  organizes  the  computations  in  such  a  way 
that  the  dependence  i6  displayed  explicitly,  so  that  integrations  with  respect  to  z'  are 
computed  more  easily.  There  is  a  bit  of  preliminary  work  that  must  be  done.  We  will  use 
our  previous  equations  and  notation  as  much  as  possible. 


C~  ® 


-A,* 


e~  V 


v  4-  y\  , 

0 

1  i  L'Ko. 

We  rewrite  the  equation  of  continuity  at  the  bottom  interface  in  the  following  form, 
which  is  suggestive  of  a  scattering  operator: 


U 

XI  £/  ' 

[t>2|v4| -Viol  -  Vso]  =-[vi|wa]  , 


(A.l) 


where 


d'  =d'e*l<i0+*'> 
f1 

c‘  =c'e_A,(,0+,) 
s1  _ — 


(A. 2) 


Equation  (A.l)  defines  the  scattering  operator  at  the  bottom  interface: 

[/'Urn 


’Me'  * 


(A-3) 


=[t/2|v4| -Uioj -Vso]  1  [-Vi  |  -  tJs] 


B12 

•#21 

B22 

B%x 

Bi2 

Bax 

Ba2 

(A.4) 


It  is  easy  to  express  these  equations  in  terms  of  the  expansion  coefficients  defined  at 


=r>  k 


(A.5) 


g~  A»(z0  +  *') 


(A.6) 


•e-M*o+*')  o  0  0] 

=(-)  0  e-As(*0+*')  0  0  c“M*o+*')  0 

S  -  o  0  105fc0  e-*«(*o+*') 

10  0  o  lJ 

’Biie"w^,<l+I ’)  #12C-(*i+*»)(*o+*') -i  '  ' 

Jg21g-(>i+>»)(*0+*')  522g-2A,(*o+*') 

=  Bsie”Al(xo+*,)  J932e_A»<*0+*') 

.1 341e-x'(‘0+z')  j?42e-A»(*o+*'> 

A  similar  analysis,  starting  at  the  top  interface  and  working  down  to  z'(+),  yields 

■ a ' 

b  ^(+)  [dl  f.  7v 

c  =5  f  ,  (A.7) 


^+)  = 


■T„eA‘*  T„ex** 

T2ieA‘*'  Tne^' 

rsleJAl*'  rS2e(A>+As>*' 

.r4  ie<A‘+A»>''  Ta  2e2XtI' 


Tn  T\  2' 
Tai  7m 
7si  Tj2 
.  Ta\  Ta 2  . 


=  [V2o|®4o|  -  ®1 1  -  Vs]"1  [V2IV4] 


(A.8) 


(A.9) 


is  the  scattering  operator  at  the  top  interface.  The  discontinuity  equation  at  z  =  z' 
becomes 

c  -  c' 
d-d * 


e  —  e 

i  f-r 


=  [®i  l®<]  V-J, 


(A.10) 


or 


(All) 


c  —  c  =hi 
d-d'  =h2 
e  -  e'  =hs 

/-/'  =  *>4, 

and  the  h’ s  are  independent  of  zo  and  z' .  Don’t  forget  that  we  have  a  different  set  of  h's  for 
each  type  of  current  (electric  or  magnetic),  and  for  each  orientation  of  the  current  source. 

Now,  on  to  the  algebra!  When  written  out,  (A. 7)  and  (A. 5)  become: 

«=s<^+s<+>/ 

e-Sj+'i  +  Sj+V 
e=SiVd+si?f 


<C=sSrV  +  5<rV 


(A- 12) 


>12 


/• -SfrV  +  SfcV 
a  -sfrV  +  sfrV 


h  =sirV  + 


>41 


42 


e  , 


and  when  (A. 11)  is  substituted  into  the  middle  four  equations  in  (A. 12),  we  get,  after  some 
steps 

and-a12f  =ri 
—  021  d  +  O22/  =r2  > 


(A- 13) 


where 


«>.  =1  -  sir’#’  -  s<r>s<+> 

=5Sr^« 1 + sir’siJ 1 
0ll  =s<r)s<+> + ;+> 

»»  =1  -  ’  -  ^r'4+> 

<•.  =Aj  -  ^r’Ai  -  ■sir’** 

r2  —  Sj, 


(A-14) 


*~*h,  -  Stfh,. 


When  we  expand  an  and  022,  using  (A. 6)  and  (A. 8),  we  find  that  they  are  independent 
of  z',  but  we  find  that  an  is  proportional  to  exp(Ai  —  Aj)z'  and  a2i  is  proportional  to 
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exp(— (Aj  —  As)*').  Hence,  when  we  solve  (A. 13)  using  Crainer’6  rule,  the  denominator 
determinant  is  independent  of  s'  (as  we  knew  it  had  to  be).  This  means  that  z'  appears 
only  in  various  exponential  terms  in  the  numerator  of  the  solution.  This  makes  it  easy  to 
integrate. 

Let’s  continue.  After  solving  (A. 13)  we  substitute  (d,  /)  into  the  first  two  equations 
in  (A. 12),  and  get  what  we  are  really  looking  for:  the  field  above  the  slab.  The  result  is 
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where 

Nx  =ex'z'h2  -  fciJ?ne-2X“°e-x‘*'  -  fisB12e-{A‘+A»>*°e-^z' 

N2  =ex**' hi  -  hiB2ie~(Xl+x*)*° e~Xl*'  -  hsB22e-2Xa*0e-x»*' 

D  =(1  -  BnTue-™'*0  -  B12Tiie~(Xl+Xz>'°)(  1  -  B2iTi2e^Xl+x^Zo  -  B22TA2e-2Xatl>) 

-  (B21Tsle”(A‘+x»>*«’  +  H22T4ie_2A**0)(HnT82e“2Al,°  +  J?i2T42e~^,+A*^°). 

(A.16) 

Equation  (A. 15)  is  needed  for  computing  Gi2,  where  the  subscripts  refer  to  regions  of 
space  (1  being  above  the  slab,  and  2  being  the  interior  of  the  slab).  In  order  to  compute 
G22>  we  need  the  coefficients  c,  e,  d! ,  /'.  These  are  gotten,  after  solving  (A. 13),  by  using 
(A.ll)  and  (A. 12): 

c  =$<+></ +$£>/ 

e-s&’rf+SSJV 

d,=d-h2 


(A.17) 


There  is  another  way  of  getting  the  results  of  this  note  that  reflects  the  physics  of 
the  problem.  It  is  less  complicated  algebraically,  and  will  produce  expressions  that  are 
easily  integrated  with  respect  to  z'.  The  technique  »  based  upon  the  notion  of  incident 
and  reflected  Green’s  functions,  and  uses  the  scattering  operators  that  have  already  been 
defined. 

We  start  with  (A. 11)  and  decompose  the  fields  within  the  slab  into  incident  and 
reflected  parts.  Because  the  h’s  in  (A. 11)  are  independent  of  z0  and  z',  being  properties 
of  only  the  exciting  point  source  and  the  material  within  the  slab,  they  are  used  to  define 
the  incident  field  within  the  slab: 


f  h2V2e  M*-*' ')  -f.  fi4v4e  if  z'  <  z  <  0 

\  —hiVieXl(*~x')  —  h3v3ex*(x~x'\  if  —  zq  <  z  <  z'. 


(A. 18) 


Note  that  the  incident  field  is  decomposed  into  an  upward-traveling  wave  above  z1 ,  and  a 
downward-traveling  wave  below  z1. 


The  incident  field  contains  all  those  singularities  of  the  Green’s  function  that  are 
associated  with  the  source  and  6ink  points.  The  reflected  Green’s  function,  therefore,  is 
regular,  and  is  the  response  due  to  the  presence  of  the  boundaries  at  z  =  0  and  z  =  —zq. 
The  reflected  field  is  given  by  the  remainder  of  the  total  proposed  field  and  the  indicent 
field: 


Fr  =cvtex'<M-x')  +  ev3ex^x-M,)  +  dW^'-0  +  fv^e 
=c"vieXlX  +  e"v3eXxZ  +  d"v2e~x'x  +  /"C4e~A*z 


where 


d"  =  d'eXlX> ,  f"  =  feXxX' 
c"  =  ce~XlX' ,  e"  =  ee-A**\ 


(A.19) 


(A.20) 


Note  that  Fr  is  continuous  across  z\  and  that  the  (c",  e”)  terms  are  downward¬ 
traveling  (which  means  that  they  originate  at  the  upper  boundary,  z  =  0),  whereas  the 
(d",  f")  terms  are  upward-traveling  (which  means  that  they  originate  at  the  lower  bound¬ 
ary,  Z  =  -20). 

Eventually  we  will  have  to  integrate  Fr  with  respect  to  z  and  z'.  The  integral  with 
respect  to  z  is  trivial,  being  simply  the  integral  of  the  exponential  terms.  The  integral 
with  respect  to  z',  on  the  other  hand  involves  the  coefficients,  c",  d",  e",  and  these 
are  only  implicitly  defined  in  terms  of  z'.  This  is  the  same  problem  that  motivated  these 
notes  in  the  first  place;  now  we  are  going  to  develop  another  formalism  for  expressing  these 
coefficients  in  terms  of  z\  in  order  to  evaluate  the  integrals. 

In  terms  of  our  new  notation,  (A.3)  and  (A. 4)  are  equivalent  to 
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and  (A. 7)  and  (A. 8)  produce 
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(A. 22) 


The  physics  implied  by  (A. 21)  and  (A. 22)  is  interesting.  The  first  term  in  (A. 21) 
corresponds  to  the  reflection  from  the  lower  boundary  of  the  downward-going  part  of  the 
incident  field,  and  the  second  term  represents  the  reflection  from  the  lower  boundary  of 
the  downward-going  wave  that  has  been  reflected  from  the  top  boundary.  In  (A. 22)  the 
terms  correspond  to  reflections  from  the  top  surface. 

Equations  (A.21)  and  (A.22)  can  be  solved  for  c",  d" ,  e",  f"  in  several  ways.  One 
way  is  to  use  an  iterative  scheme,  which  starts  out  by  assuming  that  c"  =  0,  e"  =  0  and 
then  computing  [d",  /"]  from  (A.21).  This  result  is  then  substituted  into  (A.22),  from 
which  an  updated  value  of  [c",  e"]  is  obtained.  This  completes  the  first  cycle;  the  second 
cycle  starts  with  the  updated  [e",  e"j  being  substituted  back  into  (A.21).  This  iterative 
method  mimics  the  multiple  reflections  that  sire  produced  in  this  system. 

An  alternative  method  of  solution  is  simply  to  create  a  4  x  4  linear  system  from  (A.21) 
and  (A.22),  which  is  then  solved  conventionally.  No  matter  how  (A.21)  and  (A.22)  are 
solved,  these  results  are  all  that  are  needed  to  compute  G22  (don’t  forget  the  depolarizing 
term).  In  order  to  compute  G12,  we  need  [a,  fc],  as  before.  Rather  than  use  (A. 15),  we 
can  derive  a  second  equation  from  (A. 7)  and  (A. 8): 
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(A. 23) 


This  equation  can  be  solved  once  (d",  /"]  is  known. 

Because  the  coefficient  matrices  in  (A.21),  (A.22)  are  independent  of  z',  we  can  easily 
derive  an  equation  for  the  integrals  of  [c”,  d",  e",  f")  with  respect  to  z\  In  fact,  the 
equation  is  (A.21),  (A.22),  with  the  inhomogeneous  terms  replaced  by  the  appropriate 
integrals  with  respect  to  z'.  Such  integrals  are  easily  computed  because  of  the  presence  of 
the  exponentials. 

Let’s  charge  on  and  solve  (A.21)  and  (A.22)  the  old-fashioned  way,  by  elimination  of 
pairs  of  unknowns.  For  example,  if  we  substitute  (A.22)  into  (A.21),  thereby  eliminating 
[c",eM]  in  favor  of  [d",f"]t  we  get 


where  B  is  the  matrix  that  appears  in  (A.21),  and  T  is  the  matrix  that  appears  in  (A.22). 
The  solution  of  (A. 24)  is 
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and  when  this  is  substituted  into  (A.22),  we  get 
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Because  we  are  only  solving  2x2  systems,  the  inverses  can  be  easily  carried  out  by  using 
Cramer’s  rule,  as  was  done  in  solving  (A. 23). 

The  solution  vector  can  be  written  as 


(A.27) 


where  the  matrix  elements  are  independent  of  z',  and  are  called  reflection  coefficients.  Note 
that  there  is  a  matrix  of  reflection  coefficients  for  each  of  the  three  components  of  electric 
or  magnetic  current  density,  because  h  depends  upon  the  current  source. 

In  terms  of  the  2x2  matrices  defined  earlier,  we  have 
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=  -  [T -AT]"1!. 

The  Dyadic  Green’s  Function  G [‘2\z,z') 

The  “electric-electric”  dyadic  Green’s  function  has  the  following  interpretation:  its 
first  column  is  the  electric  field  vector  produced  by  a  point  source  of  electric  current  that 
is  oriented  in  the  ^-direction,  its  second  column  is  the  electric  field  vector  produced  by 
a  point  source  of  electric  current  that  is  oriented  in  the  y-direction,  and  similarly  for  the 
third  column.  Of  course,  we  are  still  working  in  the  Fourier  domain,  so  that  the  only 
spatial  variables  are  ( z,z '),  where  z  is  the  field  point,  and  z’  the  source  point.  Hence,  we 
can  write 

G(ee)(z,z')  = 
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The  superscripts  on  the  field  components  in  each  column  denote  the  direction  of  the  applied 
electric  point-current  source.  The  delta  function  term  that  appears  in  the  zz  component 
of  the  dyad  is  due  to  the  fact  that  Ex  contains  a  term  that  is  directly  propo  v  ir  1  to  the 
s-component  of  the  applied  electric  current  density.  Hence,  when  this  current  density  is  a 
delta  function,  as  it  is  when  computing  the  dyadic  Green’s  function,  the  same  delta  function 
appears  in  the  Green’s  function.  This  term  is  the  “depolarization”  Green’s  function. 

The  tangential  components  of  the  electromagnetic  field  that  appear  in  (A. 29)  are 
computed  from  (A. 18)  and  (A. 19),  where  the  coefficients  in  (A. 19)  are  given  in  (A. 27). 
We  conclude,  therefore,  that  the  matrix  in  (A. 29)  consists  of  two  parts:  an  incident  part, 
due  to  (A. 18),  and  a  reflected  part,  due  to  (A. 19).  Hence,  we  write  G^ee^(z,z')  =  G/(z  — 
z')  +  Gd(z  —  z')  -|-  G r(z,z').  Note  that  the  incident  and  depolarization  Green’s  dyadics 
depend  upon  the  difference  between  the  source  and  field  points. 

Let’s  go  into  the  computation  of  the  dyadic  Green’s  function  in  a  little  more  detail. 
We  start  with  (A. 18)  and  (A. 19),  keeping  in  mind  that  there  are  three  different  h's  (corre¬ 
sponding  to  each  possible  orientation  of  the  point  source  of  current)  for  the  electric  Green’s 
function,  and  three  different  h’s  for  the  magnetic  Green’s  function.  We  are  interested  in 
only  the  electric  dyadic  Green’6  function  now.  Hence,  each  h  will  give  us  a  column  of  the 
Green’s  dyad;  how  do  we  get  the  rows?  Remember  that  the  eigenvectors  stand  for 
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Hence,  the  entries  of  (A.29)  can  be  easily  picked  out  of  this  vector  structure.  For 
example,  the  first  column  of  the  incident  part  of  (A.29)  is  obtained  upon  referring  to 
(A.18): 
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In  order  to  compute  the  reflected  part,  we  start  with  (A. 19),  and  get 

E  <*)  =c"<*>  vJV1*  +  e'^uJV**  +  dw(*>^1)e-A‘*  +  (A.32)(a) 

£$x)  =cM^*^v[2)eAl*  +  e"(x)vJ2)eA**  -f  rf"(l)t42)e"Al*  +  /"(*)^2)e_A»*  (A.32)(b) 
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where  [c"<*\  e"<*),d"(*>,  /"(*>]  is  given  by  (A.27),  with  the  h’s  given  by  [h(2x),  h(4x),  h(,x)]. 

It  should  be  clear,  now,  how  to  compute  the  other  two  columns  of  the  incident  dyadic 
Green’s  function. 

The  Vector  Green’s  Function  Gife^(z,z') 

We  start  with  the  dual  of  the  bottom  row  of  the  dyad  in  (A. 29),  which  produces  the 
z-component  of  the  magnetic  field: 

H  ,(*,*')  = 

[(-*,Fir)  +  ksEyX^)/ufio  {-kyE{y)  +  kxEiy))/wfi0  (-kyEix)  +  kmWf))lup*  ]  • 

(A. 33) 

If  S(kx,ky)  is  the  transfer  function  corresponding  to  integration  over  the  sensor  coil,  then 

G(,?*>(z,  z')  =  S{k„  t').  (A.34) 

This  is  the  expression  that  is  used  in  (12)  (see  the  comment  below  (3)(c),  also).  The 
electric  field  is  evaluated  in  region  1  of  the  figure  (z  >  0),  so  that  we  need  [a,  6]  from  either 
(A.15)  or  (A.37). 
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Consider  (10);  with  a  whip  source  that  produces  1  y- directed  electric  field,  Eq,  in  the 
principal- axis  system,  we  have 

EMF,m  =  -L,m-  M)E0LUJt^j.  (A. 35) 

L=0  M=0  J=  0 


Now  jg  pven  jn  (n)  and  (12),  which  uses  the  y-component  of  (A. 33)  and 

(A.34) 

G<e'M(z,z')  =  S(kt)ky)(-kyE^  +  kzEM)/u>no,  (A. 36) 

where,  from  the  figure, 
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E^  =  -  a^ajoe"^  +  bMe~x,>z. 


(A-37) 


^sr—  a(»)  and  fc(»)  are  given  by  either  (A. 15)  or  (A.23)  with  [h]  =  For  example,  in 

(A. 16)  we  have 
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(A. 38) 


whereas  (A.23)  and  (A.27)  produce 
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Clearly,  we  are  using  separation  of  variables,  in  which  the  coefficients,  [a,  6],  carry  the 
z'  dependence.  Hence,  in  evaluating  the  integrals  over  z'  in  (12),  we  need  only  integrate 
these  coefficients  with  respect  to  z'.  The  results  are  readily  apparent  from  (A. 38)  and 
(A.39),  because  the  z'  dependence  is  explicitly  stated  in  the  exponentials;  the  coefficient 
matrices  are  independent  of  z'. 
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CHAPTER  V 


LABORATORY  DATA  COLLECTION 
AND  MODEL  INTERFACE 
FOR  FLAW  INVERSION 


Jeff  C.  Treece 


1.  INTRODUCTION 


Collecting  and  storing  laboratory  data  became  a  fairly  complicated  task  as  a  large 
amount  of  data  were  taken  under  a  wide  range  of  physical  circumstances.  This  chapter 
gives  an  overview  of  the  software  and  hardware  involved  in  the  complete  data  acquisi¬ 
tion  process.  We  first  discuss  the  data  collection  software,  which  touches  on  the 
hardware  design  for  a  complete  understanding  of  the  motivation  behind  the  software. 
We  also  review  the  sensors  and  excitations  used  for  the  project;  these  sensors  are  men¬ 
tioned  here  to  complete  the  documentation  of  handling  and  keeping  track  of  the  data. 
Finally,  we  present  some  sample  laboratory  data. 


2.  DATA  COLLECTION  SOFTWARE 

The  data  acquisition  itself  was  just  one  part  of  the  overall  process,  and  began  with  a 
“C”  language  data  acquisition  program  designed  to  run  on  MS-DOS  computers.  We 
used  two  different  MS-DOS  computers  during  the  project,  and  the  software  was 
developed  to  operate  well  with  almost  any  type  of  MS-DOS  computer,  even  one  without 
a  hard  drive.  During  the  project,  we  upgraded  from  a  PC-XT  compatible  computer  to 
12MHz  PC-AT  compatible  laboratory  computer,  allowing  us  to  acquire  data  much  more 
efficiently.  Some  changes  to  the  data  acquisition  program  were  required  to  get  the  pro¬ 
gram  to  run  on  the  new  computer,  these  changes  are  documented  in  the  Sixth  and 
Seventh  Quarterly  Reports. 

Data  collected  in  the  laboratory  went  through  many  handling  steps  before  its  usage 
in  the  computer  modeling  program.  This  section  outlines  the  steps  that  we  used  in 
acquiring  and  storing  the  data  collected  in  our  laboratory  at  Indiana  University  (IU).  The 
laboratory  was  set  up  and  operated  through  IU’s  “Partners  in  Applied  Research”  (PAR) 
program.  Hardware  used  in  the  project  was  developed  prior  to  this  work,  and  is  dis¬ 
cussed  in  the  reports  for  our  contract  with  Naval  Surface  Warfare  Center,  contract 
number  N60921-86-C-0172,  which  expired  in  December  1988  [SI].  It  will  be  necessary 
to  refer  to  features  of  the  hardware  in  this  report,  but  for  detailed  information,  the  reader 
is  encouraged  to  study  the  reports  for  the  NSWC  project  and  Reference  [T2].  It  is 
sufficient  for  this  project  to  model  the  electronics  system  as  in  the  block  diagram  of  Fig¬ 
ure  V-l.  It  is  very  important  to  note  that  the  signals  produced  by  the  electronics  are  not 
the  true  in-phase  and  quadrature  EMF  readings  required  by  the  model  [S2],  but  are 
derived  from  the  EMF  values  after  some  amount  of  amplitude  scaling  and  phase-angle 
rotation.  The  actual  EMF  values  must  be  determined  from  the  laboratory  measurements 
by  forming  a  linear  combination  of  the  measured  “zero”  and  “ninety”  signals.  Refer¬ 
ence  to  Figure  V-l  is  made  later  in  this  chapter  when  discussing  the  algorithms  that  we 
used  for  computing  the  actual  EMF  values. 


V-l 
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Figure  V-l.  Model  of  the  mixer  circuit. 


2  a.  Steps  in  the  Acquisition  Process 

To  obtain  a  sufficiently  large  base  of  data  from  different  flaws  requires  special  tools 
and  methods  for  compactly  storing  and  retreiving  the  data.  All  data  collected  were  even¬ 
tually  stored  on  our  Alliant  UNIX  machine,  and  the  UNIX  environment  was  used  to 
access  the  growing  database,  process  and  plot  the  data,  and  maintain  information  about 
the  laboratory  experiments.  The  following  list  summarizes  the  steps  acquiring  data  and 
making  it  ready  for  the  model. 

1.  Collect  measurements  in  remote  laboratory.  We  will  use  the  term  “test”  to  denote  a 
given  set  of  measurements  resulting  from  one  scan  over  a  sample  of  material.  A  test 
typically  consists  of  measurements  at  a  multitude  of  frequencies  on  a  grid  of  points  in 
the  X-Y  space,  and  is  stored  on  a  floppy  disk  or  a  hard  drive  in  a  number  of  different 
files. 

2.  Transfer  the  data  from  laboratory  computer  to  Alliant  computer.  This  task  was  best 
performed  by  first  using  an  archive  program  to  consolidate  the  data  into  one  file  (we 
used  the  public-domain  program  zoo),  then  copying  the  archive  file  to  an  IU  UNIX 
machine  that  could  be  used  to  send  data  to  our  Alliant  by  UUCP  (UNIX-to-UNIX 
copy).  Alternatively  we  sometimes  hand-carried  floppy  diskettes  to  an  IBM-AT 
networked  to  our  Alliant. 

3.  Store  the  data  in  an  orderly  manner.  Our  method  was  to  use  a  unique  directory  in  the 
labdata  account,  named  by  the  year  and  date  that  the  test  was  performed,  that  con¬ 
tained  a  few  standard  files:  the  compressed  archive  in  one  of  several  formats,  tar,  zoo, 
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or  ar.  Readme,  overview  information  about  the  measurement;  db.entry,  database 
information  in  a  specific  format;  and  header,  specific  information  about  the  measure¬ 
ment.  Readme  and  header  were  also  added  to  the  archive. 

4.  Add  db.entry  to  the  database  so  that  page  number,  data  measurement  statistics,  sam¬ 
ple,  excitation,  and  other  information  could  be  accessed  through  the  database  pro¬ 
gram. 

5.  Run  a  program  called  profile  developed  by  Bishara  Shamee  and  Jeff  Treece  to  make  a 
profile  page  of  grayscale  images  to  add  to  our  booklet  of  data. 


2  b.  Phase  Shifts  and  Gain  Errors 


Note  that  there  are  two  DC  outputs  of  the  circuit  represented  by  Figure  V-l,  whose 
phase  and  amplitude  depend  in  part  on  several  phase  errors  and  the  amplifier  gain.  In 
short: 


V, 


=  sin  ( to  t ) 

=  A  sin  (cor  +  &  ) 


=  ALPF< 


sin  ( cor  +  )  sin  ( cor  +  <j>,  +  <t>A  )  \ 


-  A '  cos  ( 4>w  -  <j>4  -  $4  )  +  V^o 


(V-O) 

(V-l) 


(V-2) 


v2=alpf 

-A'  cos 

Where  Vj  and  V2  are  the  outputs  of  the  in-phase  and  quadrature  mixers,  respectively, 
and  (j>,  is  the  phase  of  the  signal  applied  to  the  mixer  input.  Equations  (V-2)  and  (V-3) 
can  be  used  to  calculate  system  parameters  when  some  of  the  variables  are  known.  Also, 
it  is  possible  to  calculate  the  mixer  phase  shift  using  the  SPICE  circuit  analysis  program. 
One  can  also  extract  the  amplifier  phase  shifts  from  the  data  sheets.  Finally,  another  set 
of  two  equations  can  be  obtained  after  modifying  the  mixer  circuit  to  bypass  the  amplifier 
(thus  bypass  the  amplifier  phase  shift).  These  equations  and  calculations  yield  an  over¬ 
determined  system,  however,  the  equations  are  non-linear,  and  several  ideal  assumptions 
have  been  made.  Moreover,  die  laboratory  measurements  are  taken  with  some  degree  of 
experimental  error.  The  solution  to  the  set  of  equations  was  attempted  using  M1NPACK. 
A  more  complete  discussion  of  this  calculation  was  presented  in  the  Seventh  Quarterly 
Report. 


|  sin  ( cor  +  +  <J>G  )  sin  ( cor  +  <k  +  <|>A  )  j 

(  +  $<?  ~  ~  )  +  Vos90  (V-3) 
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2  c.  New  Method  of  Finding  Phase  Shifts,  Based  on  Field  Model 


The  MINPACK  method  of  determining  phase  shifts,  reviewed  in  this  chapter  and 
described  in  the  Seventh  Quarterly  Report,  has  certain  drawbacks,  one  of  which  is  its 
complexity.  The  complexity  of  the  process  makes  the  error  margin  uncertain.  Because 
of  this  fact  and  that  we  wanted  to  test  results  using  an  independent  method,  we  briefly 
looked  at  a  second  method  of  calculating  phase  shift  and  gain  multipliers  for  the  lab  data. 
The  second  method  was  based  simply  on  calculating  the  expected  EMF  values  1  and 
determining  from  the  calculations  the  proper  scaling  factor  and  phase  angle  for  the  actual 
frequencies  measured. 

To  calculate  the  expected  EMF,  we  must  know,  apriori,  the  properties  of  the 
material.  For  this,  we  selected  various  metals,  whose  conductivities  were  documented  in 
[Wl].  We  selected  thin  planar  samples  of  the  metals  and  measured  the  transmitted 
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Figure  V-2.  Sensor  and  excitation  arrangement  for  calibration  experiments. 


field  using  the  setup  of  Figure  V-2.  We  were  able  to  comput  the  expected  field  using  the 
experiment  geometry  and  sample  conductivity  as  inputs  to  the  model.  The  outputs  of  the 
model  were  assumed  the  correct  EMF  values,  and  the  outputs  of  the  lab  system  were 
assumed  to  be  Vj  and  V2  as  shown  in  Figure  V-l.  A  simple  linear  transformation  can  be 
used  to  transform  Vi  and  V2  into  the  calculated  EMF  values.  Then,  assuming  the  same 
hardware  setup  for  other  tests,  other  laboratory  data  can  be  converted  using  the  same 
linear  transformation. 

Plots  of  Figure  V-3  show  laboratory  data  and  model  calculations  for  the  field 
transmitted  through  a  sheet  of  aluminum  foil.  Similar  experiments  were  repeated  for 
other  samples:  copper,  stainless  steel,  iron,  and  thick  aluminum.  Thick  samples  were 
desirable  because  the  thickness  was  convenient  to  measure,  but  were  undesirable  because 
they  attenuated  the  signal  at  high  frequencies  (demonstrating  their  EM  shielding  proper¬ 
ties).  Consistent  results  were  obtained,  indicating  that  a  good  agreement  is  seen  after 
scaling  the  lab  data  by  a  factor  of  about  10  and  correcting  the  phase  by  almost  180 

1  The  calculation  itself  is  not  simple,  but  the  concept  is  simple.  The  calculation  is  based  on 
work  reported  in  Chapter  I. 
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Figure  V-3.  Laboratory  data  (top)  and  model  calculations  (bottom)  for 
transmitted  field  through  thin  sheet  of  aluminum  at  5MHz.  On  the  left  is  the 
“zero”  signal;  on  the  right  is  the  “ninety”  signal.  Scaling  and  rotating 
(phase  angle)  the  lab  data  gives  a  close  approximation  of  the  model  calcula¬ 
tions.  The  required  rotation  and  scaling  are  not  alarming.  See  text  for  a  full 
explanation. 


degrees.  The  180  degree  phase  error  is  not  distrubing:  such  an  “error”  could  easily  be 
the  result  of  connecting  the  sensor  or  excitation  with  the  wires  reversed  or  losing  a  minus 
sign  in  a  calculation.  The  consistency  of  these  experiments  indicates  that  this  method 
might  be  an  excellent  way  to  calibrate  the  laboratory  data.  The  two  biggest  problems  in 
obtaining  accurate  results  seem  to  be  getting  an  accurate  conductivity  estimate  and 


accurately  measuring  the  thickness  of  the  material. 


2  d.  Using  High-Frequency  Data 

At  high  frequencies,  phase  shifts  and  gain  errors  in  the  laboratory  hardware  can  not 
be  ignored.  In  competition  with  the  hardware  limitations  are  the  needs  of  the  computer 
model,  which  uses  a  multi-frequency  algorithm  that  benefits  from  a  wide  range  of  fre¬ 
quencies.  By  compensating  in  software  for  phase  shifts  and  gain  errors  and  redesigning 
some  of  the  sensors  and  exciting  coils,  we  were  able  to  collect  data  over  a  wider  range  of 
frequencies  than  expected  at  the  onset  of  the  project.  The  inductance  of  our  sleeve  exci¬ 
tation  source  was  about  80jiH,  and  the  sleeve  had  a  self-resonance  frequency  in  the  1- 
2MHz  range.  In  order  to  expand  our  usable  frequency  range,  we  first  looked  for  an  exci¬ 
tation  source  with  less  inductance  (though  data  collected  with  the  sleeve  are  not  useless: 
it  is,  however,  more  difficult  to  characterize  the  current  that  is  driven  through  the  coil, 
and  thus  determine  the  phase  and  magnitude  of  the  excitation  current).  We  next  used  a 
whip  excitation,  which  had  nearly  the  same  cross-sectional  geometry  as  the  sleeve,  but 
was  formed  of  only  one  turn  of  wire  instead  of  about  forty.  We  were  able  to  collect  data 
using  the  whip  at  frequencies  up  to  25MHz.  Smaller  sensors  also  allow  us  to  make  meas¬ 
urements  at  higher  frequencies  -  sometimes  up  to  50MHz  2 .  We  were  not  able  to  collect 
data  above  50MHz  because  that  is  the  highest  frequency  we  could  generate  with  our 
HP81 16A  signal  generator. 

The  small  excitations  do  not  produce  uniform  fields,  and  the  first  version  of  the 
computer  model  assumed  uniform  fields.  Later,  changes  were  made  to  the  model  (see 
Chapters  I  and  III)  that  allowed  us  to  use  data  from  whip,  ring  and  other  sensors  and 
excitations.  We  discovered  several  trade-offs  when  using  the  smaller  sensors:  they  typi¬ 
cally  gave  higher  resolution,  measured  only  the  near-surface  region  of  the  sample,  tended 
to  be  more  sensitive  to  lift-off,  and  gave  a  smaller  flaw  signal.  Since  the  small  sensors 
gave  higher  resolution  measurements,  we  were  often  able  to  detect  micro-features  of  the 
material,  such  as  fiber  weave  structure  and  fiber  orientation. 

The  software  phase  and  gain  compensation  took  place  in  three  stages:  we  first 
characterized  the  laboratory  harware  by  performing  experiments,  then  calculated  parame¬ 
ters  required  for  rotating  and  scaling,  and  then  we  converted  the  data  by  calculating 
actual  EMF  values,  taking  the  excitation  current  into  account.  The  first  step  was 
hardware-dependent:  for  each  new  sensor/excitation  arrangement,  new  experiments  were 
required  to  characterize  the  phase  shifts.  For  this  reason,  EMF  values  could  be  accu¬ 
rately  calculated  for  only  certain  sensors  and  excitation  (the  ones  for  which  we  per¬ 
formed  the  required  experiments).  The  second  step  was  accomplished  by  running  a  pro¬ 
gram  using  the  data  of  the  first  step  as  input.  The  program  calculated  best-fit  values  for  a 
number  of  unknowns.  The  entire  method  is  described  in  detail  in  the  Fifth  Quarterly 
Report  for  contract  number  N60921-86-C-0172,  between  Sabbagh  Associates  and  Naval 
Surface  Warfare  Center  [SI].  An  alternative  method  for  the  second  step  was 


2  The  data  at  high  frequencies  had  much  more  error,  that  is  the  phase  and  amplitude  were  much 
less  accurate  than  at  the  low  frequencies. 
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contemplated  near  the  end  of  this  project,  and  some  preliminary  results  are  presented  in 
Section  0.  At  certain  frequencies,  the  final  step  can  possibly  fail  (see  Seventh  Quarterly 
Report,  Appendix  D,  Section  e).  If  the  final  calculation  fails  due  to  large  specific  values 
of  phase  shifts,  the  data  am  not  useful  at  that  frequency  for  use  in  the  model. 


3.  SENSORS,  EXCITATIONS,  SAMPLES,  AND  DESCRIPTIONS 

We  tested  new  sensors  to  optimize  the  design  for  flaw  and  feature  detection.  Cer¬ 
tain  sensor  configurations  were  particularly  good  at  detecting  certain  features.  For  exam¬ 
ple,  a  directed  current,  such  as  our  whip  excitation,  was  found  to  be  good  at  detecting 
fiber  “tows.”  3  Certain  innovative  sensors,  such  as  the  “figure-eight”  sensor/excitation 
described  in  the  sensors  listing,  were  useful  for  detecting  flaws  and  other  small  features. 
We  have  noted  that  certain  sensor  designs  are  optimum  for  particular  applications;  we 
have  already  discussed  an  excitation  that  is  good  at  detecting  signals  from  tows,  which 
indicates  micro-structure  of  the  material.  However,  carrying  the  example  further,  large 
tow  signals  hinder  detection  of  flaws,  so  those  sensors  that  are  good  at  detecting  tows 
might  not  be  good  at  detecting  impact  damage. 


3  a.  List  of  Sensors  Used 


SENSORS 


1.  Sensor  "HW1" 

Pancake  coil,  13  turns,  28  AMG.  Diameter  0.16'' 

2.  Sensor  "HWla" 

HWl's  replacement.  Pancake  coil,  13  turns,  28  AWG. 
Outer  diameter  0.125";  Inner  diameter  0.06" 

3.  Sensor  "HW2" 

Pancake  coil,  4  turns,  34  AWG.  Diameter  0.09" 

4.  Sensor  "HW3" 

Pancake  coil,  3  turns,  34  AWG.  Diameter  0.06" 

5.  Sensor  "HW4" 

Pancake  coil  3  turns,  34  AWG.  Diameter  0.06" 

6.  Sensor  "HW5" 

Pancake  coil,  3  turns,  34  AWG.  Diameter  0.03" 

7.  Sensor  "HW6" 


3  Tows  are  the  result  of  the  bunches  of  fibers  that  are  used  when  the  materials  are  made.  Con¬ 
ductivity  is  better  in  the  tow  regions  than  between  tows,  resulting  in  a  striped  image.  Some  such 
images  are  presented  in  the  Fifth,  Sixth,  and  Seventh  Quarterly  Reports.  Tow  is  spinners’  cant, 
now  used  by  makers  of  graphite-epoxy.  Satin  weave,  which  has  warp  and  woof  without  a  twill,  is 
another  borrowed  textile  phrase,  taken  from  weavers  at  the  loom. 
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Pancake  coil,  3  turns,  34  AWG.  Diameter  0.04" 

8.  Sensor  "HW7" 

Pancake  coil,  2  turns,  34  AWG.  Diameter  ?" 

9.  Sensor  "macOl" 

Circuit  board  sensor,  0.6"  x  0.565".  10  turns,  each  having  a  width  of  20 
mils . 


10.  Sensor  Mmac08" 

A  circuit  board  sensor  made  up  of  8  macOl's  with  a  spacing  of  30  mils 
between  each  sensor.  Center-to-center  distance  between  sensors  is 
0.6". 


11.  Sensor  "jrb2" 

0.375”  diameter  barrel,  probe  tip  0.25”  diameter.  This  probe  has  a 
black  housing  with  three  gold  stripes.  There  is  a  female  connector  for 
coaxial  cable  in  back,  and  a  small  (less  than  0.06"  diameter)  coil  at 
front.  The  number  of  turns  in  the  probe  tip  is  unknown. 

Inductance  ?  henries.  Resistance  ?  ohms. 

12.  Sensor  "jrb2ns” 

This  is  simply  the  jrb2  probe  without  its  ferrite  shield. 

13.  Sensor  "HW8" 

Pancake  coil,  5  turns,  28  AWG. 

Outer  diameter  0.375";  Inner  diameter  0.03" 

14.  Sensor  "HW9" 


Pancake  coil,  5  turns. 

36 

AWG. 

Diameter 

0.03" 

15.  Sensor  "HW10" 
Pancake  coil,  5  turns. 

36 

AWG. 

Diameter 

0.035 

16.  Sensor  "HW11" 
Pancake  coil,  5  turns. 

40 

AWG. 

Diameter 

0.03" 

17.  Sensor  "HW12" 
Pancake  coil,  5  turns. 

40 

AWG. 

Diameter 

0.025 

18.  Sensor  "capl" 
Circuit-board  sensor. 

0.745"  x 

0.385" 

Fingers  separated  by  20  mils.  2  groups  of  fingers.  The  inner  group  has 
four  fingers,  each  outer  group  has  two.  The  outer  groups  are  connected. 

19.  Sensor  "cap2" 

Circuit-board  sensor.  0.65"  x  0.65” 

Fingers  separated  by  roughly  30  mils.  Two  separate  groups  of  three 
fingers  each.  The  fingers  zig-zag  twice. 

20.  Sensor  "cap3" 

Circuit-board  sensor.  0.5”  x  0.75" 

Fingers  separated  by  roughly  20  mils.  Two  separate  groups  of  five 
fingers  each. 

21.  Sensor  "cap4” 

Circuit-board  sensor.  0.6"  x  0.575" 

Six  groups  of  two  fingers  roughly  20  mils  apart.  The  groups  are  separate. 

22.  Sensor  "cap5" 

Circuit-board  sensor.  0.6"  x  0.55" 

Two  separate  groups  of  three  fingers,  each  separated  by  40  mils. 

23.  Sensor  "cap6" 

Circuit -board  sensor.  0.75"  x  0.55” 
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Both  sides  art  different.  One  aide  has  two  groups  of  three  fingers,  the 
other  side  has  three  groups  of  two  fingers,  where  the  outer  two  groups 
are  connected.  The  fingers  are  roughly  50  mils  apart. 

24.  Sensor  "cap7" 

Circuit-board  sensor.  0.275"  x  0.575" 

Six  separate  groups  of  two  fingers  separated  by  20  mils. 

25.  Sensor  "cap8" 

Circuit-board  sensor.  0.39”  x  0.37" 

Both  sides  are  different.  One  side  has  four  groups  of  two,  the  other 
side  has  one  group  of  four  flanked  by  two  groups  of  two.  The  fingers 
are  20  mils  apart.  One  finger  in  the  four  finger  group  is  broken. 

26.  Sensor  "HW13" 


Pancake  coil,  5  turns. 

40 

u> 

< 

Outer 

diameter 

0.15" 

27.  Sensor  "HW14" 
Pancake  coil,  5  turns, 

32 

AWG. 

Outer 

diameter 

0.1" 

28.  Sensor  "HW15" 
Pancake  coil,  5  turns. 

32 

AWG. 

Outer 

diameter 

O 

o 

<0 

+/- 

0.01 

29.  Sensor  "HW16" 
Pancake  coil,  5  turns, 

32 

AWG. 

8 

r» 

• 

N 

diameter 

0.09" 

+/- 

0.01 

30.  Sensor  "HW17" 
Pancake  coil,  5  turns. 

32 

AWG. 

Outer 

diameter 

o 

o 

VO 

3 

+/- 

0.01 

31.  Sensor  "280" 

Sue  Vernon's  86.280  probe.  The  barrel  is  a  highlighter  pen.  A  coax 
female  connector  is  on  one  end,  and  an  iron  cup,  835  mils  in  diameter, 
is  on  the  other.  Epoxy  seals  the  cup.  Inner  diameter  -  350  mils. 
Barrel  «•  2,875  mils  long.  Inductance  •*?  Henries.  Resistance  -?  ohms. 
A  label  says: 

86.280  291 
1  c  6  03 

32.  Sensor  "335" 

Sue  Vernon's  86.335  probe.  The  barrel  is  a  highlighter  pen.  A  coax 
female  connector  is  on  one  end,  and  an  iron  cup,  1"  in  diameter, 
is  on  the  other.  Epoxy  seals  the  cup.  Inner  diameter  -  450  mils. 
Barrel  «  2,140  mils  long.  Inductance  •  ?  Henries.  Resistance  -?  ohms. 
A  label  says: 

86.335  26  16  pa 
400  3b9 

33.  Sensor  "disk" 

A  disk  drive  head.  See  the  lab  book  for  details.  This  is  basically 
a  split-core  transformer,  as  far  as  we  can  tell. 

34.  Sensor  "HW18" 

One  loop,  34  AWG,  0.2"  diameter. 

35.  Sensor  "HW19" 

One  loop,  34  AWG,  0.1"  diameter. 

36.  Sensor  "eyel” 

Two  loops  with  one  turn  in  each  loop.  The  windings  oppose  so  that  the 
background  signal  is  eliminated,  outer  diameter  .13"  +/-  0.01",  34  AWG 

37.  Sensor  "XI" 

Four  handwounds  in  series.  Each  Hw  has  two  turns  of  38  AWG  wire 
around  the  peg  of  a  dual-row  header.  The  four  sensors  were  then  bent 
into  an  x  shape  such  that  opposing  legs  canceled  each  other  in  a 
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uniform  field.  This  isnsor  mss  then  place  in  a  ring  excitation. 


38.  Sensor  "rings" 

Two  handwounds  with  two  turns  of 
a  peg  of  a  dual-row  header.  The 
background  signal  is  eliminated. 

39.  Sensor  "Dl" 

Two  handwounds  with  two  turns  of 
a  peg  of  a  dual-row  header.  The 
background  signal  is  eliminated. 


34  AWG  wire  that  were  wrapped  around 
windings  oppose  so  that  the 


34  AWG  wire  that  were  wrapped  around 
windings  oppose  so  that  the 


40.  Sensor  "microl" 

Small  micro  sensor  made  at  Purdue. 


40.  Sensor  "micros" 

Small  micro  sensors  connected  in  series,  placed  side  by  side. 

50.  Sensor  "HW20" 

Two  loops,  34  AWG,  .70”  +/-  0.01”  diameter. 

51.  Sensor  "HW21" 

Two  loops,  34  AWG,  .59"  +/-  0.01"  diameter. 


3  b.  A  List  of  the  Excitations 


EXCITATIONS 


1.  Excitation  "linch  loop” 
diameter  1",  34  AWG. 

2.  Excitation  "big  aleeve" 

A  11"  x  3.88"  x  0.25”  piece  of  Plexi-glass,  wound  by  28  AWG  wire  16 
times.  Each  winding  is  0.125”  apart. 

3.  Excitation  "junk  sleeve" 

4.  Excitation  "large  aleeve" 

9”  in  length  when  flat.  39  conductors  (standard  computer  ribbon  cable) 
connected  together  to  form  a  solenoid.  Width  1.94”. 

5.  Excitation  "small  sleeve" 

9.5"  when  laid  flat.  14  conductors  (standard  computer  ribbon  cable) 
connected  together  to  form  another  solenoid.  0.7”  width. 

6.  Excitation  "bridge” 

A  circuit  board  sensor,  mac08,  is  attached  to  a  Plexi-glass  backing, 
which  measures  9.06"  x  5.06"  x  0.24".  The  mac08  leads  run  out  of  the 
bottom  edge  of  the  Plexi-glass,  back  to  another  board  with  one  200  ohm 
potentiometer  and  2  other  resistors  for  each  lead. 

7.  Excitation  "HW2" 

See  above  entry  in  SENSORS.  Osed  as  an  excitation,  HW2  was 
arranged  in  a  "ring"  configuration. 

8.  Excitation  "HW6" 

See  above  entry  in  SENSORS  and  "HW2". 
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9.  Excitation  "harp" 

An  11"  x  2"  x  0.25"  rectangular  aolanoid  wrapped  around  Plexi -glass.  It 
has  17  turns  of  28  AWG  wire,  each  spaced  roughly  0.125”  apart. 

10.  Excitation  "loop4" 

26  AWG,  diameter  0.875”.  Loop  excitations  are  used  in  aniaotrop^  tests. 

11.  Excitation  "loop21” 

12.  Excitation  "large  loop” 

28  AWG,  diameter  2.22” 

13.  Excitation  ”small  loop” 

28  AWG,  diameter  1.11" 

14.  Excitation  "macOl” 

See  above  entry  in  SENSORS. 

15.  Excitation  "monster” 

28  AWG  wire.  10”  x  10.6”.  A  30  turn  solenoid  whose  windings  are  in  the 
direction  of  the  larger  axis.  The  windings  are  spaced  five  to  six 
sixteenths  of  an  inch  apart. 

16.  Excitation  "ringl” 

28  AWG,  diameter  0.625" 

17.  Excitation  "ring2" 

18.  Excitation  ”ring3” 

Slightly  bent.  28  AWG,  diameter  1" 

19.  Excitation  "ring4” 

28  AWG,  diameter  0.1” 

20.  Excitation  "ring5" 

34  AWG,  diareter  1.5"  +-  0.06" 

21.  Excitation  "ring6" 

22.  Excitation  "ring7" 

36  AWG.  diameter  0.625"  4~  0.06" 

23.  Excitation  "ringS" 

36  AWG,  diameter  0.6"  4-  0.06" 

24.  Excitation  "ring9" 

36  AWG,  diameter  0.625"  +-  0.06" 

25.  Excitation  "ringlO" 

36  AWG,  diameter  0.64”  4-  0.06" 

26.  Excitation  "ring20" 

40  AWG,  diameter  0.125” 

Has  broken  lead  wire. 

27.  Excitation  "ring21" 

A  double  ring,  one  superimposed  upon  the  other.  34  AWG,  diameters  1" 

4-  0.06". 

28.  Excitation  "whip” 

A  single  loop  of  28  AWG  wire.  The  loop  is  8”  when  stretched  flat. 

29.  Excitation  "window" 

A  circuit  board,  7.06"  x  7"  x  0.06",  with  copper  on  top.  The  left  and 
right  edges  are  solid  copper  columns,  each  one  0.375"  wide.  Joining  the 
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two  are  regularly-spaced  copper  rows,  each  0.06"  wide,  that  completely 
fill  the  board  from  top  to  bottom.  There  are  100  rows  in  all. 

30.  Excitation  "loop5" 

28  AWG,  diameter  0.25”. 

31.  Excitation  "ringll" 

36  AWG,  diameter  0.625"  +-  0.06" 

32.  Excitation  "ringl2" 

36  AWG,  diameter  0.625”  +-  0.06” 

33.  Excitation  "fig8" 

A  sideways  figure  8.  The  leads  are  attached  to  loop  B.  Diameter  of  loop 
A  is  0.625"  •!—  0.06".  Diameter  of  loop  B  is  0.625"  +—  0.09".  The 
fig8  is  made  out  of  36  AWG  wire. 

34.  Excitation  "280" 

See  the  above  entry  in  SENSORS. 

35.  Excitation  "335" 

See  the  above  entry  in  SENSORS. 

36.  Excitation  "whip2" 

A  single  loop  of  28  AWG  wire.  The  wire  is  0.75"  long  when  stretched  to 
a  separation  of  roughly  60  mils. 

37.  Excitation  "ring23" 

34  AWG,  diameter  0.2" 

38.  Excitation  "ring24" 

34  AWG,  diameter  0.1" 

39.  Excitation  "ring25" 

34  AWG,  diameter  0.12" 

40.  Excitation  "ring26" 

34  AWG,  diameter  0.06" 

41.  Excitation  *ring27" 

42.  Excitation  "ring28" 

34  AWG,  diameter  0.33" 

43.  Excitation  "Rings" 

Two  handwounds  of  34  AWG  wire  that  were  wrapped  around  a  peg  of 
a  dual-row  header.  The  windings  oppose  each  other. 

44.  Excitation  "micro2" 

Small  micro  sensor  made  at  Purdue. 
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3  c.  Samples 


DESCRIPTIONS  OF  SAMPLES 


1.  Sample  "l” 

Sample  1  la  the  well-known  satin  weave  square  sample  that  has  twelve 
flat -bottom  drilled  holes.  The  dimensions  are  6.125"  wide  by  6.125” 
long  by  0.11"  thick.  The  top  side  has  twelve  visible  holes  (6  0.5" 
diameter  and  6  0.25”  diameter).  The  top  side  is  smoother  than  the 
bottom  side.  Weaves  of  dimensions  approximately  0.1"  are  visible  from 
both  sides.  The  pattern  appears  to  be  "under-4,  over-1."  See  the  lab 
notes  for  more  details. 

2.  Sample  "panl2" 

6.125”  x  6.125”  x  0.11"  satin  weave.  Same  weave  pattern  as  Sample  1. 
Yellow  alignment  lines  run  horizontally  as  shown  in  the  lab  notes.  The 
top  side  is  smoother  than  the  bottom  side. 

3.  Sample  "3" 

6.1"  x  6.1”  x  0.1"  satin  weave.  There  is  a  hole  almost  precisely  in  the 
center  of  the  top  side  (top  side  Is  the  aide  that  is  "roughest") .  The 
hole  is  0.2"  in  diameter,  appears  to  be  flat-bottomed,  and  is 
approximately  55mils  deep. 

4.  Sample  "pil" 

6"  x  6”  x  0.065"  checkered  cloth  pattern,  but  not  satin  weave.  Top  side 
is  rougher  than  the  bottom  side.  The  panel  flexes  more  easily  along  the 
y  axis  than  along  the  x  axis,  as  in  the  lab  notes.  A  tag  attached  to 
the  sample  says: 

10  layers 

11  mil  ohm 
T  -  67  - 

wt  -  57.1140  gm 

v  -  39.837  p  -  1.434 

high  resistivity  reasons  unknown 

5.  Sample  "pa4" 

5.9”  x  5.9"  x  0.08".  Top  side  has  eleven  holes  of  various  sizes;  see 
the  lab  notes.  The  panel  seems  rigid  in  both  directions. 

6.  Sample  ”4" 

6"  x  6"  x  0.085".  The  top  side  is  smooth,  the  bottom  side  is 
cloth-patterned.  The  bottom  has  "5245C"  and  "#4"  written  in  yellow 
grease  pencil.  There  is  a  yellow  X  connecting  the  four  corners  of  the 
sample.  We  will  call  the  "5245C”  edge  the  "0"  direction.  From  the  top, 
one  can  see  predominate  (fiber?)  lines,  diagonal  from  upper-left  to 
lower-right.  A  blue  grease  pencil  arrow  points  to  the  left  on  the  top 
side.  The  board  bends  easier  about  the  "0”  direction  axis  than  about 
the  axis  of  the  blue  arrow.  The  board  is  warped  in  the  "easy” 
direction  (about  50  mils  in  the  center) . 

7.  Sample  "2" 

11.375"x  11.5"  x  0.085"  cloth  pattern.  The  top  says  "top/front."  The 
panel  has  blue  grease  pencil  Xs  near  the  center  on  both  sides. 

Horizontal  grooves  of  various  widths  (from  0.125”  to  1”,  two  sets)  have 
been  stamped  on  the  top  in  four  rows.  The  grooves  average  less  than  5 
mils  deep.  7”  from  the  left  edge,  and  1"  and  2"  from  the  bottom  edge, 
are  2  depressions.  The  lower  one  is  roughly  10  mils  deep,  the  upper  one 
is  around  5  mils  deep.  2"  above  the  second  depression  is  a  third,  0.5” 
long  and  roughly  5  mils  deep.  On  the  back  face,  a  crack  extends  from 
top  to  bottom,  3.875"  inward  from  the  right  edge.  The  sample  is  warped 
near  this  crack,  by  (on  the  average)  5  mils.  The  whole  sample  bends 


along  tha  axes  parallal  with  the  crack  much  easier  than  along  other 


8.  Sample  "RAEl" 

13"  x  12.87S"  x  0.215"  cloth  pattern.  The  front  has  "RAEl"  and  "0-90” 
on  a  sticker  in  the  top  right  comer.  Various  irregular  depressions 
are  circled  with  yellow  grease  pencil,  see  the  lab  notes  for  more 
details. 

9.  Sample  "RAE2" 

4.94"  x  4.38"  x  0.125"  rigid  cloth  weave.  The  panel  has  blue  tape 
saying  "RAE2"  in  the  upper  right  corner.  The  back  shows  fibers  in  x 
patterns.  The  front  has  two  circled  regions.  See  the  lab  notes  for 
more  details. 

10.  Sample  "RAE3" 

6.125"  x  4.31"  x  0.195"  rigid  cloth  weave.  The  front  and  back  are  plain 
except  for  a  circled  depression  on  the  front.  The  depression  is  located 
2.5"  up  and  3.5"  over  from  the  bottom  left  corner. 

11.  Sample  "RAE4" 

7.125"  x  6.125"  x  0.395"  rigid  cloth  weave.  The  panel  has  no  diagonal 
fibers.  No  obvious  flaws  are  on  it.  Blue  tape  saying  "RAE4"  is  in  the 
upper  right  comer. 

12.  Sample  "RAE5" 

12.25"  x  11”  x  0.32"  rigid  cloth  weave.  There  is  a  tape  saying  "RAE5" 
in  the  upper  right  comer.  Also  in  the  top  right  corner  is  written 
"0+/-45. "  The  panel  has  5  layers  of  roughly  0.06"  each.  The  front  has 
two  circled  depressions.  The  upper  left  circle  says  "28.3  J.” 

13.  Sanple  "5" 

12"  x  11.88"  x  0.15".  Top  is  smooth,  bottom  is  a  canvas-type  of 
surface.  A  blue  X  is  drawn  on  the  front.  Lines  run  from  left  to  right 
on  the  top.  The  top  is  scraped  in  the  bottom  left  quadrant.  The  edges 
are  rough  and  irregular.  Fibers  run  vertically.  The  panel  is  rigid, 
except  for  vertical  axes.  Two  cracks  are  in  the  back.  One  is  0.75" 
from  the  right  edge,  the  other  is  6.75"  from  the  right  edge.  This  is 
the  unidirectional  sample  made  by  Olaf  Rask. 

14.  Sample  "bar" 

2"  x  12”  x  1.44".  The  top  is  cloth  weave,  the  bottom  is  satin  weave. 

The  sample  is  obviously  layered,  the  Middle  of  the  bottom  is  warped 
upward  by  roughly  15mils.  The  top  says  "4-10"  and  has  two  circular 
depressions:  one  squarely  in  the  center,  and  one  4  inches  beneath  the 
first.  The  top  has  shallow  grooves  spaced  around  1"  apart;  these 
grooves  are  apparently  joints. 

15.  Sample  "A" 

9"  x  9"  x  0.11".  The  top  is  smooth,  the  bottom  is  a  tight  canvas  weave. 
Fibers  run  diagonally,  +  and  -  22.5  degrees.  A  red  "A"  and  an  arrow 
indicating  fiber  direction  are  drawn  in  the  upper  left  quadrant.  The 
panel  has  no  obvious  flaws. 

16.  Sample  "B" 

9"  x  9"  x  0.11.  This  panel  has  the  same  physical  attributes  as  Sample 
A.  A  red  "B”  and  and  arrow  are  drawn  in  the  upper  left  quadrant,  as  with 
Sample  A.  Again,  no  visible  flaws  ere  present. 

17.  Sample  "ceramic" 

Two  ceramic  disks,  one  white,  one  black.  The  white  one  has  a  circular 
depression  in  it,  85  mils  in  diameter  and  30  mils  deep.  The  white  disk 
has  a  465  mil  diameter  and  is  242  mils  deep;  the  black  disk  is  465  mils 
in  diameter  and  212  mils  deep.  The  black  disk  is  abraded  near  the  edges. 


18.  Sample  "xas" 

12"  x  12"  x  0.17”.  This  sample  is  described  fully  in  the  notes  included 
in  the  package  from  England.  It  is  a  32-ply  carbon-epoxy  alab  with  nine 
0.2-ply-wide  delaminations  ranging  in  diameter  from  0.8"  in  the  upper  left 
to  2"  in  the  lower  right.  A  label  in  the  upper  right  corner  of  the  top 
face  aays: 

070  €88  /  B.C. 

XAS  914 

32  ply  0/90/+or-45 
with  delaminations 

19.  Sample  "carba" 

20.  Sample  "carbb” 

21.  Sample  "C302" 

10"  x  13.25"  x  .14"  carbon  carbon  weave  sample.  This  sample  has  it's  lower 
right  corner  cut  out.  The  sample  is  warped  down  the  middle,  however  it  is 
not  obvious.  The  front  of  the  sample  is  labeled  with  C302  written  on  it,  also 
the  direction  of  the  warp  is  indicated  here.  The  sample  also  has  a  purple 
square  painted  on  the  front . 

22 .  Sample  "copper" 

3.995"  x  3.930"  x  .024"  piece  of  copper. 

23.  Sanple  "lead” 

4.5"  x  4.6"  x  .106"  piece  of  lead. 

24 .  Sample  "stainless" 

16.88"  x  7.9"  x  .062"  piece  of  stainless  steel. 

25.  Sample  "slide” 

1.074"  x  3.135"  x  .274  carbon-carbon.  The  top  aide  has  two  large 
flaws,  also  the  top  right  corner  is  ground  down.  Backside  has  a  glossy 
appearance . 

SOLAR  CELLS 


1.  2L65-66-503 

Solarcell  marked  with  "O"  on  the  bottom  in  black.  Has  obvious  defect 
going  across  the  cell,  80  mils  up  from  the  centerline  (silver  strip  at 
the  bottom) . 

2.  5J79-63-501 

No  defects  are  obvious,  though  there  is  a  slight  hump  80  mils  up  from 
the  centerline,  with  the  silver  strip  at  the  bottom.  Also,  there  may  be 
a  defect  along  the  upper  right  side.  See  lab  notes  for  more  details. 

3.  5J69-20-06-501 

There  seems  to  be  a  defect  175  mils  below  the  centerline,  and  the  cell's 
appearance  is  "patchy,"  more  reflective  than  cells  tl,  42,  or  <4. 

4.  4L96-44-501 

No  obvious  flaws,  though  there  is  a  slightly  raised  area  about  200  mils 
above  the  centerline. 

5.  4M59-30-501 

6.  5K19-41-501 

7.  2M16-49-503 

8.  5K15-43-501 
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9.  2L87-08-501 

There  is  an  obvious  diagonal  crack  from  the  lower  left  edge  going  up  to 
the  right. 

10.  3M61— 41-503 

There  is  a  lengthwise  raised  area  300  mils  from  the  right  edge. 

11.  4M11-59-501 

This  cell  has  a  crack  on  the  lower  right  edge.  There  is  also  a  smashed 
spot  in  the  lower  right. 

'.2.  4M59-35-501 

There  is  a  crack  at  the  upper  left  edge. 

13.  3M62-65-501 

No  glaring  defects. 

14.  4M78-48-501 

No  obvious  defects. 

15.  3M72-22-501 
Nothing  obvious . 

16.  4M59-45-501 
Nothing  obvious . 

17.  4048-37 
Exposure. 

18.  4048-60 

19.  4048-14 

20.  4048-54 

21.  4048-07 


CIRCUIT  BOARDS 


1. 

7.94"  x  3.625"  x  0.06"  Copper-colored  rectangle.  See  the  lab  notes  for 
a  picture  of  it. 


2. 

9"  x  3"  x  0.06".  Top  is  silver-green,  bottom  is  dark  green.  There  are 
many  large  holes,  but  few  pinholes.  There  are  2  solder  globs  on  the  top 
and  the  bottom.  See  the  lab  notes . 


3. 

4.5"  x  5.44"  x  0.06" 

This  board  is  green.  There  is  a  2.63"  x  0.38"  projection  out  of  the 
bottom,  located  0.63"  from  the  left  edge.  The  top  has  a  5  x  7  grid  of 
contacts.  Six  chips  are  still  on  the  board.  The  back  "looks  like  a 
parking  lot."  See  the  lab  notes  for  a  diagram. 


4. 

9"  x  3"  x  0.06" 

This  board  is  a  clone  of  *2,  except  it  has  no  solder  globs,  and  its 
pinholes  are  smaller. 


5. 

5.44"  x  3.38"  x  0.06"  green  board.  The  top  and  bottom  both  have  rows  of 
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silver  near  tha  top  and  bottom,  roughly  0.125"  wide,  but  the  top  has  a 
silver  column  at  the  left  edge,  0.125"  wide.  In  the  middle  of  both 
sides  are  pinholes,  arranged  in  6  groups  of  16.  These  pinholes  cover  a 
rectangular  area  totaling  4"  x  1".  The  bottom  left  corner  of  the 
pinhole  rectangle  is  0.75"  from  the  left  and  0.75"  from  the  bottom. 

6. 

8.5"  x  5"  x  0.06"  green  board.  The  top  face  say  "top"  in  the  lower  left 
corner.  There  is  a  gold  strip  on  the  left  edge  of  the  top  and  bottom. 
The  top  has  sparse  wiring,  the  bottom  somewhat  more  wiring.  The  only 
holes  are  close  to  the  left  and  right  edges,  in  the  center.  See  the  lab 
notes  for  more  details. 

7. 

4.06"  x  2.44"  x  0.05"  brown  board.  The  left  has  a  contact  projection 
with  dimensions  1.5"  x  0.25",  located  0.5"  up  from  the  bottom  edge.  The 
front  says  "s2193626”  near  the  lower  left.  The  upper  right  and  lower 
left  corners  have  0.625"  long,  0.063"  high  cuts  extending  into  the  board 
from  the  right  along  the  top  and  bottom  edges.  Again,  see  the  lab  notes 
for  a  diagram. 


8. 

6.69"  x  6"  x  0.06"  board,  dark  green  front,  silvered  back.  There  is  a 
contact  projection  on  the  left  edge  of  the  top.  It  is  0.25"  wide,  2.5" 
long,  and  starts  1.75"  from  the  bottom. 


9. 

8.25"  x  2.88"  x  0.07"  board.  The  front  says  "44112"  near  a  resistor  in 
the  upper  left.  A  4.12"  x  0.375"  contact  projection  starts  2"  from  the 
right  edge,  near  a  capacitor.  Two  rows  of  22  vertical,  soldered, 
pinhole  columns  cover  the  front  and  back  faces.  Two  0.125"  diameter 
holes  are  positioned  at  2"  and  6.25"  from  the  left  side,  0.5"  from  the 
top. 


10.  copper 3 

16.25"  x  10"  x  0.06"  brown  board.  The  front  face  of  this  board  has 
"10424"  in  the  upper  right  quadrant.  The  back  face  has  a  mailing  label 
along  the  top  edge.  There  is  a  rectangular  hole  in  the  upper  left 
corner  of  the  front.  Gold  contacts  are  along  the  bottom  edge.  The 
wiring  and  pinhole  rows  are  predominately  vertical  on  both  sides.  The 
board  is  made  up  of  6  layers.  The  top  and  bottom  edges  of  the  front  are 
warped  upwards  by  roughly  0.3”,  and  the  left  and  right  edges  are  warped 
downwards  by  around  0.1".  The  board  is  flexible  in  all  directions. 


3d.  Database  Descriptions 


FEATURES 


1 .  tows 

Linear  disturbance  in  data  due  to  the  sample's  fiber  direction. 

When  two  or  more  fiber  directions  are  present,  a  weave-like 
pattern  appears  in  the  data. 

2.  flaws 

Any  man-made  or  natural  cut,  crack  or  blemish  located  on  the  sample 
being  analyzed. 

3.  anisotropy 

Used  to  obtain  directional  properties  of  a  sample.  A  pattern  is 
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obtained  by  using  a  single  current  loop  placed  either  on  top 
(l)Reflected  or  on  the  bottom  of  the  sample  (2)  transmitted. 

4.  foil 

A  piece  of  foil  used  to  represent  a  "target".  This  piece  of  foil  can 
be  any  size,  shape  and  at  any  location.  (See  Read. me  files  for  detailed 
information  on  location.) 

5.  aluminum 

A  piece  of  aluminum  used  as  a  sample. 

6.  orientation 

Aluminum  "targets"  cut  in  various  shapes  and  sizes.  These  "targets" 
are  then  placed  on  the  sample  being  analyzed.  Useful  for  giving  true 
directions  of  sample  with  respect  to  the  positioner. 

7.  monosensor 

A  single  sensor  being  used  for  eddy  current  detection,  as  opposed  to 
using  several  sensors. 

8.  solarcell 

A  semiconductor  solarcell  was  used  as  the  sample. 

9.  subtraction 

The  subtraction  of  one  data  set  from  another.  This  is  generally 
performed  on  data  coming  from  two  samples  of  identical  makeup  and 
orientation,  except  one  of  the  two  samples  may  have  a  man-made  flaw. 

10.  saturated 

Measured  signal  at  the  A/D  converter  exceeded  the  maximum  readable 
value.  The  data  obtained  at  these  points  are  generally  useless.  Asure 
indication  is  areas  in  the  grayscale  that  are  all-black  or  all-white. 

11.  traces 

Copper  traces  on  a  pc  board. 

12.  weaves 

Heave  patterns  can  be  seen  on  plots  associated  with  satin  weave 
graphite  epoxy  composite  materials. 

13.  SBIR88 

These  data  were  used  in  a  proposal  for  the  SBXR  1988  solicitation. 

14.  O.lufcap 

A  . luf  capacitor  was  placed  at  the  output  of  the  lowpass  filter  on 
amplifier.  The  resistor  value  used  in  the  LPF  was  10k. 

15.  2.7ufcap 

The  lowpass  filter  at  the  output  was  using  a  2.7uf  capacitor.  The 
resistor  used  was  10k 

16.  O.Olufcap 

The  lowpass  filter  at  the  output  was  using  a  . Oluf  capacitor.  The 
resistor  used  was  10k 

17.  transmitted 

The  anisotropy  test  was  done  in  transmission  (the  sensor  and  excitation 
were  on  opposite  sides  of  the  material) . 

18.  reflected 

The  anisotropy  measurement  was  reflected  (the  sensor  and  excitation 
were  on  the  same  side  of  the  material)  . 

19.  tee 

The  sensor  arrangement  was  "tee."  For  a  tee  arrangement,  the  sensor 
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and  •xcitation  art  pancake  coils,  and  the  sensor  is  placed  tangentially 
to  the  material.  The  excitation  is  placed  tangentially  to  the  sensor 
(in  a  plane  parallel  to  the  sangtle) . 


20.  impact 

Damage  is  visible  in  the  measurement  that  was  the  result  of  inpact. 

21.  cracks 

This  primarily  applies  to  solar  cells.  A  crack  in  the  material  shows 
up  in  the  measurements. 

22.  incase 

This  feature  primarily  applies  to  solar  cells.  The  measurements  were 
made  while  the  sanples  were  still  in  their  cases. 

23.  blast 

The  data  was  taken,  but  the  real  reason  for  the  measurement  was  to 
expose  the  sample  to  a  "blast"  of  excitation. 

24 .  striping 

Discontinuities  in  the  Y  direction  in  the  plots  caused  by  multiple 
sensors  covering  adjacent  areas  of  the  sample. 

25 .  lemon 

An  anisotropy  test  feature.  The  reflected  or  transmitted  field  has  an 
"oval”  appearance  (stretched  in  one  direction) . 

26.  fourfold 

An  anisotropy  test  feature.  The  reflected  or  transmitted  field  has 
four  or  more  "lobes". 

27.  balazs 

The  data  collected  is  bad  for  some  reason.  For  example,  a  sensor  may 
have  come  loose  during  the  test  and  only  half  of  the  data  are  valid. 

26.  metal 

data  were  collected  from  transmitted  field  through  a  metal  sheet  so 
that  calibration  data  could  be  collected  to  determine  gain  and  phase 
shift  of  the  amplifiers 


4.  DATA-HANDLING  SOFTWARE 

The  phase  angle  and  gain  depends  on  many  things,  including  the  particular  sensors 
and  excitations  used.  The  numerical  accuracy  of  the  data,  and  hence  the  success  of  the 
inversion  algorithms,  hinges  on  the  proper  characterization  of  these  sensors  and 
hardware.  An  improved  method  of  storing  data  was  developed  during  this  project  that 
involved  using  a  different  coding  of  data  and  a  different  data  compression  algorithm. 

As  mentioned  above,  proper  handling  an  appropriate  usage  of  the  data  can  be  as 
complicated  a  task  as  actually  acquiring  the  data.  One  difficult  problem  is  documenting 
the  experiments  well  enough  that  the  phase  shifts,  gain,  and  physical  parameters  can  be 
accessed  when  required  by  the  conversion  program.  In  addition  to  having  the  informa¬ 
tion  available,  proper  usage  of  the  data  demands  that  the  user  have  knowledge  of  how  the 
data  conversion  tools  operate.  Another  difficulty  is  locating  the  data  required  after  it  has 
been  stored  away  in  its  archive. 
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Documentation  consisted  of  more  than  notes  about  sensor  size,  distances,  and  physi¬ 
cal  locations  of  samples  and  sensors;  we  also  performed  experiments  to  characterize  the 
hardware.  Results  from  the  experiments  were  used  data-handling  software.  Descriptions 
of  sensors,  samples,  excitations,  and  data  features  were  updated  for  the  laboratory  data¬ 
base,  and  a  menu-driven  program  was  developed  to  locate  particular  data  based  on 
features  or  other  keywords.  Documentation  was  improved  for  the  software  tools  required 
for  extracting  and  managing  data  to  help  the  user  understand  the  operation  of  the  rela¬ 
tively  complicated  data  conversion  utilities. 


4  a.  Keeping  Track  of  the  Data 

Our  method  insured  a  standard  way  of  storing  the  data,  and  reduced  the  amount  of 
disk  space  required,  making  it  easy  to  recover  raw  data  for  a  numerical  experiment.  By 
placing  vital  information  in  the  database,  one  later  has  the  abiltiy  to  search  for  specific 
samples  or  other  keywords  and  immediately  locate  the  data.  The  ability  to  select  certain 
features  of  the  data  helped  us  to  gain  insight  about  the  data  collection  process.  We  can, 
for  example,  locate  all  tests  that  demonstrate  the  “tow”  feature,  and  then  see  what  sen¬ 
sors  and  excitations  were  used  to  gather  those  data.  A  listing  produced  by  the  database 
program  is  included  in  this  report.  Sample  entries  from  the  laboratory  data  have 
appeared  in  previous  quarterly  reports,  and  do  not  appear  here.  The  complete  summary 
of  all  data  appears  in  a  collection  of  outputs  from  profile ,  that  we  appropriately  call  the 
Tons-o-Data  data  book.  One  such  entry  was  made  for  every  test  made  for  this  project 
and  others.  Also  included  in  this  report  is  a  description  of  the  operation  of  the  database¬ 
searching  utilities.  Standard  UNIX  and  public  software  utilities  were  used  to  locate  and 
manipulate  the  data:  awk,  sed,  tar ,  grep,  or,  zoo ,  compress,  and  troff.  The  following 
listing  is  an  introductory  notice  that  gives  the  user  an  idea  of  how  the  database  operates. 

INSTRUCTIONS  FOR  "DATA"  USAGE 


You  have  accessed  a  simple  (really)  database  that  will  give  you 
information  about  data  acquired  in  the  lab.  Some  on-line  definition 
files  can  be  viewed  from  the  main  menu.  These  definition  files  explain 
characteristics  of  sensors,  samples,  excitations,  and  features  from  the 
entire  collection  of  lab  data.  A  few  simple  database  operators  allow 
you  to  pick  out  particular  sets  of  data  based  on  your  selection  of 
features  or  other  stored  information.  A  report  will  be  generated  that 
can  be  sent  to  the  printer,  a  computer  file,  or  the  terminal  screen. 

In  order  to  effectively  use  the  database,  you  should  understand  how  the 
information  is  stored... 

1)  Raw  information  about  lab  measurements  is  edited  into  a  form-like  file 

that  is  human  readable.  That  file  is  called  "db". 

2)  The  raw  file  is  condensed  into  a  more  compact  file  that  is  not  very 

readable.  Information  is  stored  one-data-set -per-line.  Separate 
bits  of  information  within  the  line  are  separated  with  the  "6"  sign. 
The  compact  file  is  called  "db. condensed". 

To  retrieve  information,  the  compact  file  is  "grep"ed  for  patterns  that 
match  the  keywords  specified  by  the  database  user.  The  keywords  can  be  a 
list  of  words,  and  can  involve  logical  operators:  OR,  AND,  and  NOT.  For 
example,  one  can  search  for  all  data  listing  features  "flaws"  but  NOT 
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"balms."  Not*  that  ainea  tha  flla  is  grapped  on  a  line-by-lina  basis, 
there  is  no  way  of  tailing  what  fiald  tha  pattarn  matches.  Tha  pattarn  "1" 
will  match  "1"  in  tha  SAMPLE  fiald  (saa  fialds  below) ,  but  will  also  match 
"1988"  in  tha  DATE  fiald.  If  you  wish  to  match  fialds  pracisaly,  you  can 
surround  tha  pattarn  with  "8";  thus  tha  pattarn  "818”  matchas  the  "1"  in 
tha  SAMPLE  fiald  (unfortunataly  "818"  matchas  "1"  in  other  fialds  such  as 
"AVERAGES"  —  saa  balow) . 

Tha  fialds  appaaring  in  tha  db  file  are  listed  below: 

1)  (empty)  Empty;  nothing  in  this  fiald 

2)  PAGE  Pag*  number  in  tha  Tons-o-Data  manual 

3)  SAMPLE  The  official  sample  name 

4)  DIRECTORY  The  pathname  to  tha  data  on  the  Alliant 

5)  DATE  Data  recorded  by  data  acquisition  program 

6)  SENSOR  The  official  sensor  name 

7)  EXCITATION  The  official  excitation  name 

8)  ORIENTATION  Orientation  of  sample,  if  applicable 

9)  TOPORBOTTOM  Top  — >  sansor  scanned  over  top  of  sample 

10)  RANGE  Two  numbers  representing  lower  and  upper  frequency 

11)  TEST  Name  of  tha  test  performed  (e.g.  "flaw”) 

12)  QUALITY  On*  word  describing  data  quality  (good,  fair,  poor) 

13)  DIMENSIONS  Size  of  scan  in  inches  —  X  and  Y 

14)  RESOLUTION  Step  size  in  inches  —  X  and  Y 

15)  AVERAGE  Number  of  points  averages  together  per  measurement 

16)  VERSION  Laboratory  hardware  ID  number 

17)  FEATURES  Features  demonstrated  by  the  data 

For  more  information,  refer  to  documentation  for  the  data-taking  program 
and  documentation  for  tha  conversion  program. 


4  b.  Data  Conversion 

To  correct  in  software  for  phase  and  gain  errors,  a  program  was  run  in  the  lab  to 
record  the  outputs  of  the  A/D  converters  (all  channels)  under  simplified  conditions. 
These  data  were  used  determine  the  phase-shift  and  gain  effects  present  in  the  circuitry. 
A  program  called  adcounts .  c  reads  the  A/D  converters  on  the  lab  computer  at  a 
number  of  frequencies  to  provide  required  phase  and  gain  information.  A  second  pro¬ 
gram  is  a  c-shell  program  on  the  Alliant  UNIX  computer  that  takes  adcounts  output 
and  formats  a  “parameter”  file  to  be  used  by  the  parameter-calculating  program.  A  third 
program,  calcparms.f,  uses  the  MINPACK  equations  and  the  parameter  file, 
params.dat  to  determine  amounts  of  scaling  and  rotation.  The  output  of  the 
calcparms  program  is  called  calcparms.dat.  A  fourth  program,  makesystem 
takes  the  calcparms.dat  file  and  produces  mixer,  parms.  Finally, 
mixer  .parms  is  used  by  convert  when  the  phase  and  gain  calculations  are  made. 
This  process  is  further  described  above,  in  the  DATA  COLLECTION  section. 

convert  was  used  to  translate  raw  laboratory  data,  stored  in  a  compact  coded  for¬ 
mat,  into  user-ready  EMF  data,  processed  as  needed  to  compensate  for  the  phase  and 
gain  errors  described  above.  A  great  deal  of  input  is  required  by  convert  before  the 
EMF  values  can  be  extracted.  The  values  of  Nx  and  Ny,  along  with  the  number  of  sen¬ 
sors  used  and  the  ordering  of  the  data  files,  are  always  required  when  attempting  to 
reorder  the  data,  convert  must  also  know  the  storage  format  of  the  data:  several  dif¬ 
ferent  storage  formats  were  used  throughout  this  project  (see  Image  Compression  and 


Archival).  Also,  the  laboratory  setup  and  parameters  are  needed  to  use  the  data  in  the 
model  (e.g.  the  list  of  frequencies  and  excitation  voltages  used).  If  the  user  attempts  to 
compensate  phase  and  gain  or  normalize  to  the  exciting  coil  current,  then  the  values  of 
amplifier  gains  and  exciting  coil  voltages  are  required.  Other  information  is  useful  for 
keeping  track  of  how  the  data  is  collected:  date  of  the  measurement,  version  of  data- 
taking  program  used,  number  of  measurements  averaged  to  get  one  data  value,  sensor 
used,  excitation  coil  used,  the  sample  material  measured,  and  the  physical  dimensions  (in 
inches,  for  example)  of  the  area  covered  by  the  test. 

In  early  versions  of  the  data-taker  program,  the  data  files  were  named  with  digits 
indicating  the  sensor  number  and  Y  value  of  the  pass.  This  storage  is  very  different  than 
the  storage  method  of  later  versions.  Whenever  a  such  a  revision  is  made  to  the  lab 
software,  the  required  method  of  conversion  changes.  If  a  new  there  is  a  new  hardware 
setup,  the  sensor  order  is  possibly  reversed,  and  in  some  cases,  sensors  can  be  hooked  up 
with  reverse  polarity.  Thus  for  every  hardware  or  software  revision,  a  new  hardware  ID 
number  was  created  that  was  used  as  an  input  to  convert.  The  program  can  decide 
from  the  unique  ID  how  to  reorder  the  data  and  invert  certain  channels  if  required.  All  of 
the  above  parameters  are  read  by  convert  in  the  form  of  a  header,  which  is  always 
created  by  the  data-taker  program.  The  following  list  is  a  summary  of  the  information 
found  in  header: 


Hardware  ID  Number 
Date  of  Experiment 

Description  (Sample,  Sensor,  Excitation) 
Beginning  X  Value  in  Inches 
Ending  X 

Number  of  Data  Values  in  X 
Beginning  Y  Value  in  Inches 
Ending  Y 

Number  of  Data  Values  in  Y 
Number  of  Frequencies  Used 
Number  of  Sensors  Used 
Number  of  Data  Points  Averaged 
hist  of  Frequencies  in  KHz 
List  of  Exciting  Coil  Voltages 
List  of  Gain  Settings  (0-H;  1-L) 
Data-taker  Version  ID 


Included  in  the  list  of  inputs  to  convert  is  the  revision  level  of  the  data-taking 
program.  The  revision  of  the  data-taker  is  important  because  a  number  of  factors  within 
the  data-taker  can  affect  the  accuracy  of  the  measurements.  If  a  measurement  is  repeated 
for  verification,  it  is  important  to  use  the  original  revision  of  the  program.  Also,  when 
the  laboratory  numbers  are  compared  to  model  calculations,  the  accuracy  of  the  measure¬ 
ments  is  an  important  consideration.  Under  some  conditions,  the  data-taker  might  give 
erroneous  results,  and  the  revision  level  helps  debug  problems  with  the  data. 

A  number  of  hardware  system  parameters  must  be  determined  before  the  phase  and 
amplitude  of  the  signal  can  be  calculated  (see  DATA  COLLECTION).  These  parameters 
are  assumed  to  be  stored  in  a  file  called  mixer .  parms  in  the  lab  data  home  directory. 
When  convert  runs,  it  searches  for  the  file  and  allows  phase  and  gain  calculations  if 
the  file  is  present.  If  the  file  is  not  present,  then  the  program  will  warn  the  user  that  no 
phase/gain  calculation  can  be  attempted.  Sometimes,  convert  can  not  perform  the 
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requested  calculation  of  EMF  values.  In  this  case,  the  files  are  still  created  but  are  not 
normalized:  the  result  is  “raw”  values  from  the  A/D  converter.  Thus,  the  user  must  pay 
attention  to  any  warnings  given  by  the  program. 

The  contents  of  mixer. parms  are  calculated  by  calcparms:  mixer  phase 
shift,  exciting  coil  current,  exciting  coil  magnitude,  amplifier  phase  shift  ( H  and  L  gains), 
individual  channel  DC  offsets  (H  and  L  gains),  and  individual  channel  gains  (H  and  L 
gains).  The  values  of  the  parameters  change  with  different  exciting  coils,  sensors, 
cables,  and  mixer  circuits;  thus,  a  new  mixer.parms  file  is  required  each  time  a  new 
experiment  is  set  up  (whenever  the  actual  phase  and  magnitude  is  required).  Since  it  is 
not  always  feasible  to  produce  a  new  parameter  file  every  time  an  experiment  is  run,  we 
often  "live  with”  the  phase  and  gain  errors  that  result.  Any  time  that  the  “standard” 
parameter  file  is  not  accurate  enough,  additional  data  can  be  taken  to  improve  the 
experiment’s  accuracy.  It  is  impossible  for  convert  to  estimate  the  accuracy  of  the 
conversion,  except  when  the  frequency  is  known  to  be  out  of  range. 

Documentation  for  convert  exists  for  the  FORTRAN  version  of  the  program. 
When  the  format  of  data  storage  was  changed  (data-taker  Revision  1.17),  convert  was 
re-written  in  C  language  to  assure  better  portability  to  other  computers,  such  as  MS-DOS 
machines.  The  operation  of  the  C  version  of  convert  closely  follows  the  documentation 
in  two  internal  documents:  A  Guide  to  Using  Convert,  and  Convert  Usage  Instructions, 
both  written  by  Jeff  Treece.  Some  amount  of  usage  instruction  is  given  when  invoking 
convert  -h. 


4  c.  Image  Compression  and  Archival 

Counts  measured  from  the  A/D  converters  were  originally  stored  in  jeffcode.  4  A 
slightly  more  compact  storage  would  result  from  storing  a  pair  of  readings  in  three  bytes, 
but  it  would  be  impossible  to  edit  or  view  the  file  with  standard  text  tools.  The  disk- 
space  savings  would  be  25  percent,  ignoring  the  newline  character  at  the  end  of  every 
line  in  the  coded  file.  The  standard  file  format  for  all  hardware  Revision  Numbers  4.1 
and  4.4,  prior  to  software  Revision  1.17,  has  two  lines  of  jeffcode  numbers  for  each  X 
value  in  the  file:  one  line  contains  the  in-phase  measurements  and  the  other  line  contains 
the  quadrature  measurements.  Each  of  the  two  lines  has  a  coded  measurement,  occupy¬ 
ing  two  characters,  for  each  frequency  recorded.  Each  line  is  terminated  with  the  UNIX 
newline  character  (10  decimal).  Jeffcode  must  be  first  decoded  since  it  is  not  a  native 
machine  data  type.  Coding  and  decoding  jeffcode  can  be  done  using  bit  operations  to 
speed  up  the  process,  but  the  process  always  takes  some  amount  of  computation.  A 
decoded  number  is  always  in  the  range  of  0-4095  decimal,  inclusive,  and  represents  the 
number  of  counts  measured  by  the  A/D  converter.  The  analog  input  range  of  the  A/D 
converter  is  -10V  to  10V;  thus  the  counts  are  convened  into  a  voltage  by  the  equation 
V  =  -10  +  C  x  20/4096,  where  C  is  the  number  of  counts  measured  by  the  A/D 

4  Jeffcode  is  described  in  previous  quarterly  reports.  It  was  created  for  convenience,  not  for 
efficiency.  Every  byte  of  the  coded  file  is  a  readable  ASCII  character,  making  it  possible  to  edit 
and  view  die  file. 
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converter. 

Data  for  each  measurement  are  stored  in  a  compressed  archive.  Three  archive  pro¬ 
grams  have  been  used:  UNIX’s  ar,  UNIX’s  tar,  and  the  public  domain  program 
zoo.  Zoo  has  its  own  internal  data  compression;  when  the  other  two  are  used,  the 
resulting  archive  file  compressed  with  UNIX’s  compress  utility.  Data  from  the  lab 
stored  in  jeffcode  does  not  typically  have  a  good  compression  rate  (not  much  disk  space 
is  saved)  because  The  like-features  of  die  data  are  dissociated:  the  spatial  structure  of  the 
image  is  not  maintained.  One  way  to  improve  the  compression  is  to  store  the  gradient 
(along  columns)  of  the  file.  This  simple  operation  improves  compression  by  saving  only 
variations  and  not  the  absolute  numbers. 

Another  way  to  improve  compression  is  to  store  the  data  in  a  reordered  space.  The 
normal  ordering  (adjacent  X  values  stored  in  successive  memory  locations)  makes  the 
most  sense;  if  stored  this  way,  regions  of  small  variations  compress  well.  There  are  other 
advantages  to  storing  data  in  the  sensible  ordering;  for  example,  the  data  can  be  read 
directly  into  arrays  using  low-level  read  statements,  resulting  in  much  more  efficient  I/O 
in  processing  programs.  Compression  can  often  be  improved  by  separating  least- 
significant  bits  and  most-significant  bits  into  separate  files.  The  most  significant  bits 
have  little  variation,  and  compress  very  well.  This  discussion  applies  only  to  general- 
purpose  data  compression  utilities,  such  compress.  There  may  exist  other  compres¬ 
sion  utilities  that  take  advantage  of  the  geometry  of  the  data  file  that  would  perform  ade¬ 
quately  without  reordering  the  data.  Prior  to  Revision  1.17  of  the  data- taker,  storage  was 
in  jeffcode-,  after  that  time,  the  data  were  written  with  the  four  most-significant  bits  in  one 
file  and  the  eight  least-significant  bits  in  another  file.  The  naming  convention  for  the  data 
files  became: 

A  number  of  experiments  were  performed  using  the  techniques  discussed  above. 
Timing  results  are  not  very  accurate  because  they  were  performed  on  a  UNIX  machine 
while  it  was  tending  to  other  tasks.  The  timing  can  be  used  to  get  a  rough  idea  of  the 
computer  resources  required  to  store  and  extract  the  data.  Experiments  using  several  dif¬ 
ferent  methods  yielded  the  following  results: 
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Lab  Data  Files: 


Output  of  Convert: 


namelO.Oh 


namelO.OO 


Base  Name  | 

Frequency 

Phase 

Bit  Significance 


Base  Name  | 
Frequency 
Phase 


Figure  V-4.  File  naming  convention  for  the  data  files  and  convert  out¬ 
put  files.  The  “Base  Name”  is  up  to  five  characters,  and  is  determined 
when  the  data  is  taken.  “Frequency”  is  an  integer,  starting  with  0,  that  indi¬ 
cates  the  frequency  index  (0  is  the  first  frequency,  1  is  the  second,  and  so 
on).  “Phase”  is  either  0  or  1  in  the  lab  data  files  and  00  or  90  in  the  con¬ 
vert  output  files,  and  represents  whether  the  data  are  in-phase  or  quadra¬ 
ture.  “Bit  Significance”  is  either  h  or  1,  and  indicates  whether  the  file  con¬ 
tains  most-significant  bits  or  least-significant  bits.  Lab  data  files  are  stored 
in  character  format  (one  byte  per  datum),  and  convert  output  files  are 
stored  in  single-precision  floating  point  binary  format 


# 

Storage  Method 

JC 

DF 

Size 

Time 

1. 

Raw  Jeffcode 

y 

n 

424 

35.1 

2. 

Differenced  Jeffcode 

y 

y 

252 

30.3 

3. 

Reordered  Jeffcode 

y 

n 

220 

25.9 

4. 

Differenced  Reordered  Jeffcode 

y 

y 

200 

24.3 

5. 

Reordered  2-byte  Integer 

n 

n 

276 

27.4 

6. 

Differenced  Reordered  2-byte  Integer 

n 

y 

208 

27.8 

7. 

Reordered  2-byte  Integer  Separate  Files 

n 

n 

216 

25.5 

8. 

Differenced  Reordered  2-byte  Integer 
Separate  Files 

n 

y 

176 

28.3 

9. 

1.5-byte  Character  No  Compression 

n 

n 

384 

2.6 

10. 

2-byte  Character  No  Compression 

n 

n 

500 

11. 

1.5-byte  Character  Differenced  Reordered 

n 

y 

168 

33.3 

12. 

1.5-byte  Character 

n 

n 

204 

32.3 

13. 

2-byte  Character  Differenced  in  X-Y  Reor¬ 
dered  Separate  Files 

n 

y 

140 

26.1 

*  Very  approximate 


Table  V-l.  Experimental  timings  and  disk  usage  values  for  various  storage 
methods. 


V-25 


Hie  column  labeled  “JC”  indicates  whether  or  not  the  method  uses  jeffcode ;  if  so, 
some  computation  is  required  to  extract  the  data.  Likewise,  the  column  “DF”  indicates 
whether  or  not  the  data  is  differenced.  Again,  a  “y”  in  this  column  indicates  that  addi¬ 
tional  computation  is  required  to  extract  the  original  data.  The  size  is  given  in  k-bytes 
used  by  the  archive  file.  The  time  column  tells  how  much  time  was  required  for  extrac¬ 
tion,  in  seconds.  The  time  column  is  very  approximate  since  the  time  depends  on  system 
load  and  other  factors.  “Jeffcode”  is  the  data  format  produced  by  the  laboratory  pro¬ 
gram.  “Reordered”  means  that  the  data  has  been  re-written  into  a  “normal”  ordering 
(see  above),  whose  rows  and  columns  correspond  to  X-  and  Y-indices  of  the  data.  “2- 
byte”  and  “1.5-byte”  indicate  how  many  bytes  are  used  to  store  one  datum.  “Integer” 
indicates  that  the  counts  of  the  A/D  converter  are  re-formatted  into  the  machine  represen¬ 
tation  for  an  integer.  “Character”  indicates  that  the  data  arc  stored  as-is  from  the  A/D 
converter,  bit-for-bit.  “Separate  Files”  indicates  that  the  most-significant  and  least- 
significant  bits  arc  stored  in  two  different  files.  In  terms  of  disk  storage  space  required, 
method  13  is  the  best.  In  terms  of  extraction  time,  methods  9  and  10  are  the  best.  Other 
methods,  such  as  2-byte  character,  might  be  more  convenient  for  accessing  and  using  the 
data.  The  new  version  of  software  uses  method  7. 


5.  SAMPLE  LABORATORY  DATA 

It  is  beyond  the  scope  of  this  report  to  present  all  500  plus  pages  of  the  Tons-o-Data 
booklet.  We  have  presented  some  of  these  data  in  previous  quarterly  reports,  and  a  rela¬ 
tively  complete  discussion  of  the  data  is  presented  in  the  final  report  for  the  NSWC  pro¬ 
ject,  mentioned  previously  in  this  chapter.  For  the  purpose  of  this  report,  a  few  key  sets 
of  data  are  presented  (the  data  presented  here  bears  special  significance  because  it  is  con¬ 
sidered  in  model  calculations  in  Chapter  ID).  The  images  are  presented  here  as  they 
appear  in  the  Tons-o-Data  booklet.  The  data  used  by  the  model  are  these  data,  rotated 
and  scaled  as  discussed  above. 


5  a.  Data  Database  Listing 

A  fairly  complete  summary  of  the  stored  data  is  presented  here.  The  database  list¬ 
ing  starting  on  the  following  page  was  produced  by  the  database  programs,  and 
represents  one  of  the  forms  of  output  that  a  user  can  request  when  searching  for  particular 
data.  Page  numbers  in  the  listing  refer  to  the  page  numbers  from  the  Tons-o-Data  book¬ 
let,  as  in  the  Sample  I  and  Sample  3  data  presented  in  this  chapter. 


*•« 

Dbwctory 

f— pit 

Oaalana 

Scmot 

Satire* 

Q 

Baa. 

Range 

Tati 

001 

1907*0(1 9/mpU.n.' 

A 

nflacud  anisotropy 
ZTufeap 

aaaeOO 

small  loop 

good 

0.1  01 

50016000 

anisotropy 

003 

1907*0(1 9*nplb  jef 

B 

mfloctad  anisotropy 
ZTufeap 

maeOO 

small  loop 

good 

0.1  0.1 

50016000 

anisrtropy 

003 

1907*0(1 9*mpU  tran 

A 

anisotropy  mnnninad 
ZTufeap  tenon 

tnacOO 

small  loop 

good 

0.1  0.1 

S0016000 

anisotropy 

007 

1907*0(1 9*mplb.tnn 

B 

aniaouopy  nanatnnod 
ZTufeap  lamon 

maeOO 

small  loop 

good 

0.1  0.1 

50016000 

anisotropy 

009 

19C7M*c01 

B 

flaws  loan  SBIROS 
ZTufeap  impact 

macOS 

arhip 

good 

0.1  0.1 

100025000 

flaws 

Oil 

1907/doclO 

aiuniman 

reflected  aniaooopy 
ahaninian  Z7ufeap 

mac08 

1  inch  loop 

fair 

0.1  0.1 

1000  25000 

anisotropy 

aatiuauan 

013 

1907/decll 

ahaninum 

reflected  aniaooopy 
ahan  ininn  Z7ufeep 

maeOO 

linch  loop 

good 

0.1  01 

1000  25000 

anisotropy 

015 

19*7/Uacl4 

aluminum 

acflectad  aniaooopy 
ahaninian  Z7ufctp 

maeOO 

linch  loop 

good 

0.1  0.1 

SO  1250 

anisotropy 

017 

1917M«c02*amplel 

1 

flaw*  ZTufeap  monoaen- 
aor  aaurated 

maeOO 

arhip 

good 

0.1  0.1 

1000  25000 

flaws 

019 

1987/dec02A*mpIeA_B 

AB 

flaws  Iowa  foil  Z7ufeap 
monoaouar 

mac08 

arhip 

good 

0.1  0.1 

1000  25000 

flaws 

021 

19t7A*c07 

I 

flan*  mmnttannv 
SBIROS  Z7ufcap 

maeOO 

whip 

good 

0.1  0.1 

1000  25000 

flaws 

023 

1907*1110 

1 

flaw*  ZTufeap  striping 

mac08 

bridge 

good 

0.1  0.1 

1009000 

bridge 

025 

19t7*dl4Ampll.ll0tnils 

1 

flawa  Z7ufcap  sniping 

maeOO 

bridge 

fair 

0.1  0.1 

1009000 

bridge 

027 

I9t7*tfl4/panl25mils 

pan  12 

flawy  Z7ufcap  sniping 

maeOO 

bridge 

fair 

0.1  0.1 

1009000 

bridge 

029 

1907/50115 

2 

flawa  Z7ufcap  striping 

maeOO 

bridge 

fair 

0.1  0.1 

1009000 

bridge 

031 

1907/jull6 

2 

flawa  Z7ufcap  atnping 

maeOO 

bridge 

fair 

0.1  0.1 

100  25000 

bridge 

033 

1907*1123 

1 

flawa  Z7ufctp  atnping 

maeOO 

bridge 

fair 

0.1  0.1 

100031000 

bridge 

035 

1907*il24 

1 

flawa  Z7ufc*p  atnping 

maeOO 

large  sleeve 

fair 

0.1  0.1 

1009000 

flaws 

037 

1907*1100 

1 

flawa  Z7ufeap  atnping 

mac08 

bridge 

fair 

0.1  0.1 

1009000 

bridge 

039 

1 9t7*il09*ranpD 

2 

flaw*  Z7ufctp  sniping 

maeOO 

bridge 

poor 

0.1  0.1 

1009000 

bridge 

0«1 

1907*dO9*aduide.tMmpl 

1 

flaw*  Z7ufctp  sniping 

maeOO 

bridge 

fair 

0.1  0.1 

1009000 

bridge 

M3 

1907*al2/anpie5.1oc|n 

5 

aniaooopy  ZTufeap 
lemon  uanoninod 

macOg 

small  loop 

good 

0.1  0.1 

1009000 

anisotropy 

M5 

1907/junl  5/Bnpl5_tniao_toptide 

5 

reflected  aniaouopy 
27ufc*p 

maeOO 

small  loop 

good 

0.1  0.1 

1009000 

anisotropy 

047 

1907/joal  5/amp]2_aniso_topaide 

2 

teflaaad  aniaooopy 
Z7u/cap  sniping  baJus 

maeOO 

small  loop 

poor 

0.1  0.1 

1009000 

anisotropy 

049 

1907*mI5/loap_OfUy 

anisotropy  ZTufeap 

maeOO 

small  loop 

good 

0.1  0.1 

1009000 

anisotropy 

051 

1907*m02/uniple3 

3 

flaara  Z7ufcap  striping 

maeOO 

large  sleeve 

good 

0.1  0.1 

1009000 

flaws 

053 

19t7*m02Aample.pa4 

P*4 

flaws  Z7afcap  striping 

maeOO 

Urge  sleeve 

good 

0.1  0.1 

1009000 

flaws 

055 

1907*m26/loopjcf 

aniaooopy  ZTufeap 
reflected 

maeOO 

small  loop 

good 

0.1  0.1 

1009000 

anisotropy 

057 

1907*n26*RtplS.tm 

5 

•nuecropy  transmitted 
ZTufeap 

maeOO 

small  loop 

good 

0.1  0.1 

1009000 

anisotropy 

059 

1907*ia20*mplSjcf 

5 

anisotropy  ZTufeap 
aflectad 

macOS 

small  loop 

good 

0.1  0.1 

1009000 

anisotropy 
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Mnctary 

Palm 

Stator 

Source 

Q 

Hoc. 

Rwtgr 

Tool 

061 

lOSTftatBAempM 

4 

anieotzopy  2.7ufcop 
tftormtud  fourfold 

mocOI 

•mall  loop 

good 

0.1  0.1 

1009000 

anisotropy 

063 

19S7^un09£enl2 

pul  2 

omocBepy  ZTufcop 
momioad  fourfold 

mocOS 

mall  loop 

good 

0.1  0.1 

1009000 

anisotropy 

063 

Wn/junSOAcBl 

1 

Bom  ZTufcop  otziping 
btluo 

mocOS 

junk  sleeve 

poor 

0.1  0.1 

1009000 

flaws 

067 

19l7/]Ba30AaM2 

1 

Bom  Z7ufcop  otziping 

bolozo 

mocOS 

junk  sleeve 

poor 

0.1  0.1 

1009000 

flaws 

069 

WtUfm/M 

•moooopy  2.7ufcop 

mocOS 

small  loop 

good 

0.1  0.1 

1009000 

anisotropy 

071 

19*7Aneyl9 

1 

Bom  2-7  uf cop  otziping 

mocOS 

large  sleeve 

good 

0.1  0.1 

100  5000 

flaws 

073 

19t7Aney20*tonp2o 

2 

flow*  27ufcop  otziping 

mocOS 

large  sleeve 

poor 

0.1  01 

100  5000 

flaws 

075 

19t7Anay2(Vbmp2b 

2 

Bom  ZTufcop  (taping 

mocOS 

largesleeve 

poor 

0.)  0.1 

1005000 

flaws 

077 

1  m/rnt  y2 1  ifemp2o 

2 

Bom  ZTufcop  otziping 
bolozo 

mocOS 

large  sleeve 

poor 

0.1  0.1 

1005000 

flaws 

079 

19S7Aney21Amp2b 

2 

Bovo  Z7ufcop  (taping 

bolozo 

mocOS 

lsrge  sleeve 

poor 

0.1  0.1 

100  5000 

flaws 

Oil 

19l7Aney25Anoy25.no 

2 

omoooopy  ZTufcop 
reflected 

mocOS 

large  loop 

good 

0.1  0.1 

1005000 

anisotropy 

013 

1 9t7An*y26/BnpLPi  1 

P‘1 

onioccropy  Z7ufcop 
reflected  fourfold 

mocOS 

large  loop 

good 

0.1  0.1 

1005000 

anisotropy 

013 

19t7Aney26A>nplPl  2 

pl  2 

■meotropy  Z7ufcop 
reflected 

mocOS 

large  loop 

good 

0.1  0.1 

1005000 

anisotropy 

017 

19S7Aney27Ample4r 

4 

onioauopy  ZTufcop 
reflected  fourfold 

mocOS 

small  loop 

good 

0.1  0.1 

1005000 

anisotropy 

019 

1 9*7Aney27A«nplo4z 

4 

onicetropy  trenenriticd 
Z7ufcop  fourfold 

mocOS 

■mall  loop 

good 

0.1  0.1 

1005000 

anisotropy 

091 

19t7Anay27^xnl2z 

pin  12 

onioatrapy  tmurniued 
ZTufcop  fourfold 

mocOS 

small  loop 

good 

0.1  0.1 

1005000 

anisotropy 

093 

19S7Aney27/unpl4.eniio 

4 

ogioouopy  27ufcop 
reflected  fourfold 

mocOI 

small  loop 

good 

0.1  0.1 

1005000 

anisotropy 

095 

19*7Ane»2l/jpil_2 

1*1 

oniootrapy  Z7ufcop 
tronomiued  fourfold 

mocOS 

small  loop 

good 

0.1  0.1 

1009000 

anisotropy 

097 

19t7Anoy2l/pilc 

P‘1 

onioouopy  ZTufcop 
truunutiod  fourfold 

mocOS 

small  loop 

good 

0.1  0.1 

1005000 

anisotropy 

099 

19l7Aney2S/omplZemoo 

2 

oaioouupy  Z7ufcop 
trreunuuad  lemon 

mocOS 

small  loop 

good 

0.1  0.1 

1009000 

anisotropy 

10, 

19t7Aiovll 

AB 

Bom  tom  Z7ufeop 
oaiping 

mocOS 

whip 

foil 

0.1  0.1 

1000  25000 

flaws 

103 

19t7/dec22Ancnflcrl 

1 

Bom  ZTufcop  oaiping 

mocOS 

monster 

foit 

0.1  0.1 

1000  25000 

flaws 

105 

19l7Mov02/M_whipl  J2J 

A 

tom  Z7ufcop  reaping 

mocOS 

whip 

good 

0.1  0.1 

1000  25000 

flaws 

107 

19l7Aiec22Aieipl 

1 

Bom  Z7ufeop  oaiping 

mocOS 

harp 

good 

0.1  0.1 

100025000 

flaws 

109 

1 9I7too*02/M_whip_22 5 

A 

tom  ZTufcop  otziping 

mocOS 

whip 

good 

0.1  0.1 

100025000 

flaws 

111 

19tt^on06 

oolezl 

ooloicetl  monooonoor 
Bom  ZTufcop  oaiping 
crecko 

HW1 

whip 

foil 

0.05  0.05 

100025000 

crack 

113 

1 9*7Ax>v02%*_whip_«5 

A 

tom  ZTufcop  oaiping 

mocOS 

whip 

good 

0.1  0.1 

100025000 

flaws 

115 

19S7Am»20 

AB 

Bom  tom  foil  Z7ufcop 
monooonoor 

mocOS 

whip 

good 

0.1  0.1 

100025000 

flaws 

117 

19t7Awv03AmplA 

A 

Bom  oubaoction 

ZTufcop  oaiping 

mocOS 

whip 

foil 

0.1  0.1 

1000  25000 

flaws 

! 

I 

i 

! 
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*ZL 

Mndary 

lusepie 

Source 

Q 

Rea. 

Range 

Tsai 

119 

19g7*ov03taqplB 

B 

Bam  eubtnetion 
Z7uicep  renping 

mecOS 

whip 

fair 

0.1  0.1 

1000  25000 

flaws 

121 

19g7*ovOS 

AB 

Bearttowt  Z7ufcap 
reripang 

aacOt 

whip 

fair 

0.1  0.1 

100025000 

flaws 

123 

19*7*ov06 

arienuoai  27ufctp  foil 

macOS 

arhip 

good 

0.1  0.1 

1000  25000 

orientation 

12S 

19f7Axs!3 

2 

flam  Z7ufcep  atriping 

maeOg 

arhip 

good 

0.1  0.1 

1000  30200 

flaws 

127 

19t7Ax>15 

1 

flam  Z7ufc*p  striping 

macOg 

arhip 

good 

0.1  0.1 

1000  25000 

flaws 

129 

19*7*eaOB*efro.ierl 

1 

flam  ZTufcep  striping 

macOg 

window 

fair 

0.1 0.1 

100025000 

flaws 

131 

19S7*ct23 

B 

flam  tom  ZTufeap 
striping  impact 

macOg 

aleevc 

fair 

0.1  0.1 

1000  25000 

flaws 

133 

1917*029 

B 

flam  tom  Z7ufcap 
amping  impact 

macOS 

arhip 

fair 

0.1  0.1 

1000  25000 

flaws 

13S 

19t7*eril3/dcfreeterl 

1 

flam  Z7ufcep  (taping 

macOg 

arindow 

fair 

0.1  0.1 

1000  25000 

flaws 

137 

19g7*c«03/tb_try_0 

B 

flam  tom  Z7ufcap 
sniping  impact 

macOg 

arhip 

good 

0.1  0.1 

1000  25000 

flaws 

139 

19ggfrnll 

1 

Sam  ZTufeap  monoaen- 
aor 

HW1 

ringl 

good 

0.1  0.1 

1000  25000 

flaws 

HI 

1917*030 

A 

tom  Z7ufcap  auiping 

macOg 

arhip 

good 

0.1  0.1 

100025000 

flaws 

143 

19g7Aepl  S/TMR_2 

B 

amaatropy  ZTufeap 
reflected  lemon 

macOg 

small  loop 

good 

0.1  0,1 

SOOluOOO 

anisotropy 

145 

19t7*cpl5/TMX_l 

B 

aniaoiopy  Z7ufcap 
reflected  loncn 

macOg 

small  loop 

good 

0.1  0.1 

50016000 

anisotropy 

147 

19g7Aep23/tb_tryl 

B 

flam  tom  Z7ufeap 
striping  impact 

macOg 

large  sleeve 

fair 

0.1  0.1 

50016000 

flaws 

149 

19g7/tep23/tb_try2J> 

B 

flam  tom  Z7ufcap 
auiping  impact 

macOg 

large  sleeve 

good 

0.1  0.1 

50016000 

flaws 

151 

19S7/aep23fcb_try3_90 

B 

flam  tom  Z7nfcap 
atriping  impact 

macOg 

large  sleeve 

fair 

0.1  0.1 

50016000 

flaws 

153 

19g7/Mp2gftbtryl_90 

B 

flam  tom  ZTnfcap 
atriping  impact 

macOg 

large  sleeve 

fair 

0.1  0.1 

50016000 

flaws 

135 

19t7/acp29/a«lryl0 

A 

tom  Z7ufcap  auiping 

aoacOI 

large  sleeve 

fair 

0.1  0.1 

50016000 

flaws 

157 

19f7A«p3CVWUu 

1 

flam  Z7ufcap 

macOg 

arhip 

good 

0.1  0.1 

100025000 

flaws 

139 

1 9l7*ep3GWef meter 

1 

flam  SBOttt  Z7ufcap 
oriping 

macOg 

arindow 

fair 

0.1  0.1 

1000  25000 

flaws 

161 

19t7*pi23*np3jid 

3 

flam  Z7ufoap  atriping 

macOg 

large  sleeve 

poor 

0.1  0.1 

1005000 

flaws 

163 

I9f7/epr23*npl2 

2 

Sam  Z7ufcap  atriping 

macOg 

large  sleeve 

poor 

0.1  0.1 

1005000 

flaws 

163 

19t7/ap23/umpk3.di»g 

3 

flam  Z7ufeap  atriping 

macOg 

large  sleeve 

good 

0.1  0.1 

1005000 

flaws 

167 

19l7/aprl3 

2 

flam  Z7ufcap  atriping 
balaza 

macOg 

small  sleeve 

poor 

0.1  0.1 

1005000 

flaws 

169 

19t7*prl6 

3 

flam  Z7ufeap  atriping 

macOg 

small  sleeve 

fair 

0.1  0.1 

1005000 

flaws 

171 

19t7/apr22*nplZ90 

2 

flam  Z7ufeep  atriping 

macOg 

small  sleeve 

poor 

0.1  0.1 

1005000 

flaws 

173 

19t7/epr22/nnpI3.90 

3 

flam  Z7ufcap  atriping 

macOg 

small  sleeve 

fair 

0.1  0.1 

1005000 

flaws 

173 

19l7*pr30 

1 

flam  SBDUt  ZTufeap 

macOg 

small  sleeve 

good 

0.1  0.1 

1005000 

flaws 

177 

1917/tpfOI 

2 

Bam  Z7ufeap  atriping 

macOg 

small  sleeve 

good 

0.1  0.1 

1005000 

flaws 

179 

19g7/ept09/50tnil 

ttl»c0g 

flam  Z7ufcap  atriping 

macOg 

small  sleeve 

good 

0.05  005 

1005000 

flaw 

111 

19g7/tptO9/100mil 

2 

Bam  Z7ufcap  atriping 

macOg 

small  sleeve 

good 

0.1  0.1 

1005000 

flaws 
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Q 
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1(3 

19*7Auar25 

1 

Bam  2.7ufcap  sniping 

macOS 

small  sleeve 

good 

0.1  0.1 

100  5C30 

flaws 

US 

19t7Aa*r26 

1 

Sam  2.7ufcap  sniping 
aatun tod 

macOS 

small  sleeve 

good 

0.1  0.1 

1005000 

flaws 

ir 

19t7ta*r27/nnpl2 

2 

Sam  2.7ufcap  sniping 

msc08 

small  sleeve 

poor 

0.1  0.1 

100  5000 

flaw* 

119 

19C7Au£7tenpl  1  disg 

1 

Bam  2.7ufcap  sniping 

macOS 

small  sleeve 

good 

0.1  0.1 

100  5000 

flaws 

191 

19C7Aaai27/tanpllback 

1 

Sam  2-7ufcap  (taping 
saturated 

mac08 

small  sleeve 

good 

0.1  0.1 

100  5000 

flaw* 

193 

19*7An»yll 

1 

Sam  2.7ufcap  (taping 

*wt«« 

macOS 

large  sleeve 

fair 

0.1  0.1 

1005000 

flaw* 

195 

19t7Auy07 

bar 

Sam  2.7ufcap  (taping 
impact  baiazi 

mac08 

big  sleeve 

poor 

0.1  0.1 

1005000 

flaws 

197 

19ISQael2 

ccpprl 

traces  monoeensor 
ITufcap  saturated 

HW1 

hngl 

fair 

0.02  0.02 

10025000 

peboard 

199 

19SI/janl3/salaroelll 

■olarl 

aolaroell  manoaeatcr 
cracks  saturated  2.7ufeap 
Sam 

HW2 

nng2 

fair 

0.01  0.01 

100  30000 

crack 

301 

19tS/jsnl3/pcbosrd 

pcboardl 

Iraoaa  monceensor 
2.7ufcap  u  tun  led 

HW2 

ring? 

good 

0.02  0.02 

200  25000 

peboard 

303 

191l/>nl3/iolaroel]l,2 

■olarl 

cacka  Sam  monoacnaor 
aolaiccU  tatunied 
27ufcap 

HW2 

ring2 

fair 

0.020.02 

1000  50000 

crack 

303 

19SS/j4nl4/sclsrcell2 

aolsr2 

mcnoswtsor  solarcell 
Mtuntod  2.7ufcap 

HW2 

ring2 

fair 

0.02  0.02 

1000  50000 

crack 

307 

19tl/3tnMAaUroeQ3 

aoUr3 

Bam  endu  monoeensor 
aolaroell  2.7ufcap 

HW2 

ring2 

good 

0.020.02 

1000  50000 

crack 

309 

19S8/jinl4Aolareell2.2 

aolaz2 

monoaensor  aolaroell 
2.7ufcap 

HW2 

rmg2 

fair 

0.020.02 

1000  50000 

crack 

311 

19tS4u>14Mucdl4 

aolar4 

Sam  monooenaor  solar- 
cell  2.7ufcap  clicks 

HW2 

nng2 

fair 

0.020.02 

1000  50000 

crack 

213 

19SS/jtnl4/unplel 

1 

Sam  monoaensor 
weaves  2.7ufcap 

HW2 

ting2 

good 

0.1  0.1 

50050000 

flaws 

215 

19U/jsnl5/pcboard2 

oopprl 

traces  monoaensor 
ITufcsp 

HW2 

ring2 

good 

0.01  0.01 

50050000 

trace 

217 

19tt5an2S 

1 

SamO.lufcspmonoaen- 
aor  weaves 

HW1 

HW2 

good 

0.1  0.1 

50050000 

flaws 

219 

19lt/jan29 

copi2 

traces  monoaensor 
O.lufcap 

HW2 

nng2 

good 

0.0198  0.02 

50050000 

peboard 

221 

19SS/JuI9^>c{x>*id 

oopprl 

tnoes  monoaensor 
O.lufcap 

HW2 

nng2 

good 

0.1  0.1 

50050000 

peboard 

223 

19t  Vjan21  /oopprl 

1 

traces  monoaensor 
(Xlufcep 

HW7 

ring6 

fair 

0.1  0.1 

50050000 

peboard 

223 

19tS/jin22 

aolu4 

monoaansor  aolaroell 
O.lufcap 

HW5 

HW6 

poor 

0.1  0.1 

20000  50000 

crack 

227 

lWWoblO 

1 

Bam  monoaensor 

O.lufcap 

HW1. 

ring3 

fair 

0.1  0.1 

500  50000 

flaws 

229 

19*Meb06 

1 

Bam  monoaensor 

Ollufcap 

macOl 

macOl 

good 

0.1  0.1 

1006000 

flaws 

231 

1 91  Web  1 1 

1 

Sam  monoaensor 

Olufcap  weaves 

HW3 

ringd 

good 

0.1  0.1 

50050000 

flaws 

233 

198l/jan30 

2 

Sam  O.lufcap  monoacn- 
eor 

HW3 

nng2 

good 

0.1  0.1 

50050000 

flaws 

235 

19S8/jan20 

1 

Sam  O.lufcap  moncscn- 
sormaves 

HW5 

HW6 

good 

0.1  0.1 

50050000 

flaw 

237 

HMiail94npltI 

1 

Bam  tnonoawtaor 

O.lufcap  weaves 

HW2 

hng2 

good 

0.1  0.1 

100050000 

flaw 

239 

19tt/fabl3 

aolirS 

monosotaor  (olaicell 
O.lufcap 

HW3 

nng4 

good 

0.10.1 

1000  50000 

crack 
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Senior 

Source 

Q 

Rea. 
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Teat 

Ml 

19M/M>14/hotamaU6 

aolu6 

aolamell  monoeenear 
fliwi  O.lufcap  cndu 

HW3 

ring4 

good 

0.1  0.1 

1000  50000 

crack 

M3 

19M/r«bl4/hoUKcin 

•oUr7 

aolemeU  monoeenear 
Din  O.lufcap  cndu 

HW3 

ring4 

good 

0.1  0.1 

1000  50000 

crack 

MS 

I9*WebI  5  Aolamell  1 

tolerS 

aolemeU  mcstaceniar 

0. 1  ufcap 

HW3 

nn|4 

good 

0.1  0.1 

100050000 

crack 

M7 

19M/r«i>17^npbb 

B 

flew*  manotauor 

O.lufcip  impact  low 

HW3 

rinf4 

fair 

0.1  0.1 

1000  50000 

flaw 

M9 

im/!cbl7fcnptb 

B 

fltw*  monoaenaor 

01  ufcap  impact  low* 

HW3 

ring4 

fair 

0.1  0.1 

1000  50000 

flaw 

2S1 

lWWebl7/tanpl2 

2 

flaw*  monoaenaor 

O.lufcap 

HW3 

ring4 

poor 

0.1  01 

100050000 

flaw 

253 

19U/r<fcl«Aol*ncUfih 

aolub 

aolimcll  monaaenaar 
flaw*  O.lufcap  cndu 

HW3 

nng4 

good 

0.1  0.1 

1000  50000 

crack 

255 

19SMebl6Aalcmell6 

joijob 

aolatceU  roanoacnacr 
flaw*  O.lufcip  cndu 

HW3 

ringd 

good 

0.1  0.1 

500  50000 

crack 

257 

1988/fd>15/ccpper3 

coppe3 

tract*  monoaenaor 
0.1ufcap 

HW3 

ring4 

good 

0.1  0.1 

10040000 

pc  board 

259 

19*M«b22/Smplb22 

B 

flaw*  monoaensor  towi 
O.lufcap  impact 

HW3 

ring4 

good 

0.1  0.1 

100050000 

flaw 

Ml 

19»Meb25 

1 

flaw*  mono*en*or 
O.lufcap  taunted 

HW3 

ring! 

good 

0.1  0.1 

500  50000 

flaw 

M3 

19«*(5ui07 

pctenaor 

tnen  monoeemor 
0.1ufcap 

HW1 

whip 

fair 

0.00958  0.01 

10030000 

trace 

265 

19I8Anu04Moem 

1 

flaw*  monoaenaor 

0.1ufcap 

HW3 

ringS 

good 

020.2 

50020000 

flaw 

M7 

19SI/mtiO»/ilow 

1 

flaw*  monoaonior 

O.lufcap 

HW3 

hngS 

good 

0202 

50020000 

flaw 

269 

19SSAnu21 

1 

flaw*  mcaaaenaor 

Ol  1  ufcap  tee  weave* 

HW3 

nngl 

good 

0.1  0.1 

1000  50000 

flaw 

271 

I988Antr22Aee2 

2 

flaw*  monoaencor 

0.1  ufcap  tee  taunted 

HW3 

nngl 

fair 

0.1  0.1 

100050000 

flaw 

273 

198I/mirZ2Aecb 

B 

monoaenaor  0. 1  ufcap  tee 
tow*  impact 

HW3 

ringl 

0.1  0.1 

1000  50000 

flaw 

275 

1988Anet23 

1 

flaw*  monoaenaor 
(Llufcap  tee  weave* 
taunted 

HW3 

nngl 

good 

0.1  0.1 

100050000 

flaw 

277 

1988/eptOl 

eolarl 

aolamell  erteka 
monoaenaor  O.lufcap  tee 

flaw*  taunted 

HW3 

nngl 

fair 

0.025  0.025 

50050000 

crack 

279 

1988/apiOd 

ooprlO 

meet  monoaenaor 

0.1  ufcap  toe 

HW3 

nngl 

good 

0.05  0.05 

50050000 

trace 

211 

19S8/aprl5 

AB 

foil  monoaenaor  O.lufcap 
low*  flaw* 

macO! 

whip 

good 

0.075  0.075 

5009000 

target 

213 

19*V«prll 

aolarl 

aolarcell  cndu 
monotwiaor  0. 1  ufcap 
Saw*  u  turned 

macOl 

whip 

fair 

0.05  0.05 

50025000 

crack 

2SS 

19tt/aprl9 

eolarl 

ha  lan  enefca  aolamell 
monoaenaor  0.1  ufcap 
flaw*  taunted 

macOl 

whip 

fair 

0.05  0.05 

500  25000 

crack 

2 17 

19ll/ipr20 

aolu9 

cndu  motatemor  tolar- 
cell  incite  O.lufcap 

flaw* 

HW3 

ring20 

poor 

0.05  0.05 

1000  50000 

crack 

2*9 

mt/wpm 

aoIat9 

cndu  mono* wear 
incite  aolamell  O.lufcap 
lie  we 

HW3 

nr  (20 

fair 

0.05  aos 

100050000 

crack 

291 

19tt/*p«23/nel  .1 

net 

tow*  taunted  aiufcap 
monoaenaor 

macOl 

whip 

fair 

0.1  0.1 

50050000 

flaw 

V-31 
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Ttat 

»3 

19SIAp*2S/mel  .2 

net 

•owe  monoeentor 
Olufctpbtltrj 

mtcOI 

whip 

fair 

0.1  0.1 

50050000 

flaw 

295 

19*l/tpr25/ree2 

ae2 

lowt  mcnotauor 

O  lufetp  returned 

mtcOl 

whip 

fair 

0.1  0.1 

50050000 

flaw 

297 

19MApr2S/ree3 

ac3 

impact  towt  returned 
mcnoemtcr  0.  lufetp 

mtcOI 

whip 

fair 

0.1  0.1 

50050000 

flaw 

299 

19t*/apr25/r»ed 

reed 

lowt  maootcntcr 
OLlufcep  btltzs 

mtcOI 

whip 

fair 

0.1  0.1 

500  50000 

flaw 

301 

19*lApr25/r»c5 

B«S 

imnmonnaior 

O.lufctp 

mtcOI 

whip 

good 

0.1  0.1 

30050000 

flaw 

303 

19tlApr2S 

aolai9 

cacki  man  aeon*  or 
O.lufctp  toltrcell  flew* 

jrt>2 

nng21 

good 

0.05  0.05 

50050000 

crack 

30$ 

19flAnay02/ncd 

m*4 

anisotropy  roonoMrtfar 

0.  lufetp  trantmitted 
fourfold 

jrt>2nt 

loop21 

fair 

0.20.2 

50050000 

anisotropy 

307 

19»8/roay02/rae5 

nc5 

tniaeuepy  bkxkmmoi 

0.  lufetp  transmitted 
lemon  fourfold 

jrt>2m 

loop21 

fair 

0.20.2 

50050000 

anisotropy 

309 

19tSAney02Aolar 

aolu9 

CBck*  monotone  or 
O.lufctp  toltrcell  fliws 

jrb2ns 

ring21 

good 

0.05  0.05 

50050000 

crack 

311 

19SI/mty09/umpleA 

A 

•miauopy  O.lufctp 
moootentor  tmumiaed 
lemon 

mtcOI 

loop4 

good 

0.1  0.1 

50050000 

anisotropy 

313 

19tMmay09Auga 

AB 

0.  lufetp  monotensor  foil 
mutated  lows  flew* 

mtcOI 

loop4 

fair 

0.075  0.075 

500  50000 

target 

31S 

1 911/mtyl  QAtmpleA_B 

AB 

O.lufctp  manoeontor  foil 
low*  flew * 

mtcOI 

loopd 

fair 

0.075  0.075 

500  50000 

target 

317 

19SSAntyl 1  Aee2 

ae2 

moaoeentor  O.lufctp 
transmitted  tmeeuepy 
fourfold 

mtcOl 

loop4 

good 

0.1  0.1 

50050000 

anisotropy 

319 

19llAneyl  l/in»o_A 

A 

moBoecneor  0.  lufetp 
reflected  tnitouopy 
lemon  unturned 

mtcOI 

loopd 

fair 

0.1  0.1 

50050000 

anisotropy 

321 

19tlAnayl  1  /amao.AB 

AB 

tnonoe  enter  0.  lufetp 
BBUminod  miiaunpy 
teluBied  lemon 

mtcOl 

loop4 

fair 

0.1  0.1 

50050000 

anisotropy 

323 

19UAuyl2/nel 

reel 

monoeonear  0.1ufctp 
Orwumtund  anisotropy 
fourfold 

mecOl 

loopd 

good 

0.1  0.1 

50050000 

anisotropy 

32S 

19f  SAnayl  1  A«3 

nc3 

nanoeoBsor  0.  lufetp 
trenanisod  anisotropy 

tt Untied 

mtcOl 

loopd 

fair 

0.1  0.1 

50050000 

anisotropy 

327 

191Atotyl2M  _reflccled 

A 

moneuuor  0.  lufetp 
reflected  anisotropy 
returned  lemon 

mtcOl 

loopd 

fair 

0.1  01 

50050000 

anisotropy 

329 

19tSAntylOAarget_X 

AB 

manotantor  O.lufctp 
fltw*  fen  returned  tows 

mtcOl 

loopd 

fair 

0.1  0.1 

50050000 

flaw 

331 

19MAnayl2/Bnpl4 

4 

monotmtor  0.  lufetp 
transmitted  tnitouopy 
fourfold 

mtcOl 

loopd 

0.1  0.1 

50050000 

•miauopy 

333 

1 911/mtyl  3/*ample5 

3 

monatauoe  O.lufctp 
trtntmined  tnitouopy 
returned  lemon 

mtcOl 

loopd 

fair 

0.1  0.1 

50050000 

anisotropy 

333 

IMMntyld/rempleA.B 

AB 

monotauor  0.  lufetp 
tnntmiued  tnitouopy 
fourfold 

mtcOl 

loopd 

good 

0.1  0.1 

50050000 

anisotropy 

337 

19tflAnayl7Aeel 

nel 

monotmtor  0.  lufetp 
imptet  flaw*  returned 

HW9 

ling? 

poor 

0.1  0.1 

50050000 

flaw 
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339 

19UAneyl7Me45.1 

tael 

menoaanaot  O.lufcep 
impact  flews  beleie 

warn  led  levs 

HW9 

nngl 

poor 

0.1  0.1 

S00S0000 

Aaw 

341 

19tlAiuyl7/ae43.2 

ael 

roped  flan  tom 

HW9 

ring? 

fair 

0.1  0.1 

50050000 

flaw 

343 

19MAnayl7/biasl 

solar  1 7 

moooswisor  eolaiccll 
Cllulcep  bleat 

HW9 

ring? 

good 

0.1  0.1 

30050000 

blast 

343 

19Uteayl7Aias2 

solar)  7 

0.1ufcap  blast 

HW9 

ring7 

good 

0.1  0.1 

30050000 

blast 

347 

19UAnayl7AU>3 

solar  17 

raonaeouor  eolarcell 

0. lufeap  blest 

HW9 

ring7 

good 

0.1  0.1 

30050000 

blast 

349 

lWttnayl  17*oLu9-ll 

eolat9 

mnpoeroeoc  eoleiccll 

0.  Ill/cap  cache  Aiwa 

HW9 

ring? 

good 

0.1  0.1 

50050000 

crack 

331 

19MAnay  1  IMul  -4 

aolerl 

(Llufcep  cmcke  flaw* 

HW9 

ring7 

fair 

0.1  0.1 

50050000 

crack 

353 

19((Aney26/figStest 

1 

Hen  mopoeeaieot 

0. 1  ufeep  weevea 

HW3 

figS 

good 

0.1  0.1 

50050000 

flaw 

■ 

1988AMy25AempJeeb 

AB 

anisotropy  transmitted 
fourfold  monoaoisor 
O.lufcep 

macOl 

loopS 

good 

0.1  0.1 

30050000 

anisotropy 

357 

19((Anay26Aacl  clean 

nel 

anisotropy  transmitted 
monoeataor  O.lufcep 
fourfold  ae tun  led 

macOl 

loops 

fair 

0.1  0.1 

50050000 

anisotropy 

339 

1988AMy26A»elflew23J 

ael 

anisotropy  transmuted 
mooeaansor  0.  lufeap 
fourfold  flam  impact 

macOl 

loops 

fair 

0.1  0.1 

50050000 

anisotropy 

361 

1988/junOl 

I 

tnenceenear  a  lufeep 
weaves 

HW10 

ring7 

good 

0.1  0.1 

50050000 

flaw 

363 

1988/ym03/cipl 

1 

moatoemeor  O.lufcep 
saturated  flaws 

eapl 

eapl 

fair 

0.1  0.1 

500  50000 

flaw 

363 

1988/)un03/caplr 

1 

nueirameor  O.lufcep 
flaws 

capl 

capl 

fair 

0.1  0.1 

50050000 

flaw 

.  367 

19(S/jtin03/cop2 

1 

tnonoeanaor  0.  lufeep 
flaws 

cap2 

cap2 

fair 

0.1  0.1 

30050000 

flaw 

369 

19((/jun03/cep5 

1 

moateeroaof  0.  lufeep 
flaws 

cap5 

eapS 

fair 

0.1  0.1 

50050000 

flaw 

" 

19tt/jun06/v«mcn 

1 

nsooasausor  O.lufcep 
Htunted  flaws 

2(0 

335 

fair 

0.1  0.1 

50050000 

flaw 

373 

19(SJ5un06fcep2s 

1 

■nonoeeBsor  a  lufeep 
flaws 

cep2 

cap2 

fair 

0.1  0.1 

300  50000 

flaw 

375 

19tl/fun06/cep6 

1 

monosaneor  O.lufcep 
flaws 

cap6 

eap6 

fair 

0.1  0.1 

500  50000 

flaw 

377 

19tt/junQ7/v«nMn280 

1 

moneaeosor  0.  lufeap 
flaws 

2(0 

whip 

fair 

0.1  0.1 

30050000 

flaw 

|  379 

19((/jun07/vum<ii335 

1 

flaws 

335 

whip 

fair 

0.1  0.1 

50050000 

flaw 

311 

t9(lfjun09/'eightcheanels 

1 

0.  lufeap  flours  sen  pang 

mseOS 

heap 

fair 

0.1  0.1 

30050000 

flaw 

1 

19t8/junl3/h»ip3lch 

1 

O.lufcep  atripinfl  flows 

mecOS 

harp 

fair 

0.1  0.1 

30050000 

flaw 

315 

19((/junl3/eolet 

aolu5 

crsdte  solaroeU 
monosmsos  0.  lufeap 
bloat  flours 

HW13 

Iing7 

good 

0.0$  0.05 

30050000 

crack 

317 

19My>al4A]uod_iing 

1 

flaws  weaves  striping 
O.lufcep  aetuntioc 
be  lets 

HW9 

ring7 

fair 

0.1  0.1 

30050000 

flaw 

" 

19Hfrml4/RAEl 

ael 

flaws  impact  sniping 
so  tun  ti  on  O.lufcep 

balsa 

HW9 

ring? 

fair 

0.1  0.1 

30050000 

flaw 

391 

19(fc5url4Ajued_ting2 

ael 

So  are  impact  sniping 
setunuon  O.lufcep 
balezs 

HW9 

ring7 

fair 

0.1 0.1 

30050000 

flaw 

V-33 


hr 

Mraciary 

ra 

Faaluraa 

Scaur 

Source 

Q 

Has. 

Rauft 

Teal 

393 

I9*l/jun09/rwxpped0„  1 

i 

flaws  striping  0.  lufcap 

macOS 

hup 

fair 

0.1  0.1 

50050000 

flaw 

395 

19Uftal(flastf 

i 

balazs  striping  0.  lufcap 
flaws 

HW14 

m i|7 

fair 

0.1  01 

$0050000 

flaw 

397 

19tt^»21AeA6 

i 

flaws  striping  0.  lufcap 

macOg 

hup 

fair 

0.1  0.1 

$00  50000 

flaw 

399 

19M0un21  /drive_hesd 

i 

flaws  0.  lufcap  monoacn* 
•or 

Aimk 

good 

0.1  0.1 

$05000 

flaw 

401 

19tVjuo2Via«lbolLcmcap6 

tael 

flaws  0. 1  ufcap  monoacn - 
•or 

cap6 

csp6 

fair 

0.1  0.1 

$00  50000 

flaw 

403 

19Mfiua34/Mmnic_di«b 

cm 

flaws  0.1  ufcap  monoaen- 
aor 

eap6 

cap6 

fair 

o.os  aos 

$00  50000 

flaw 

40$ 

19M/jua22/iolar 

solu6 

flaws  0. 1  ufcap  monoacn* 
aoraolarccll 

<4t«lf 

A'tlf 

good 

0.0$  0.05 

50  5000 

flaw 

407 

19*l/jun23/nel_c«p6 

nel 

flaws  0.1  ufcap  monoacn  - 
■or 

cap6 

cap6 

fair 

0.1  0.1 

50050000 

flaw 

409 

19tI/junZ3/ncl  disk 

nel 

flaws  0.1  ufcap  monoaen- 
■or 

disk 

disk 

fair 

0.1  0.1 

505000 

flaw 

411 

19tl/jun2g/recl  _4ch*nneli 

r*el 

flaws  0.1  ufcap  striping 

HW14 

ring! 

fair 

0.1  0.1 

50050000 

flaw 

413 

194t/jun2k/xxs_4chxnnd* 

B4S 

flaws  0.1  ufcap  striping 

HW14 

rin|7 

fair 

0.1  0.1 

50050000 

flaw 

413 

19tt/jull  1/kanplel  _weaves 

1 

flaws  0.1  ufcap  monoaen* 
•or  weaves  saturated 

HW2 

hngl 

fair 

0.05  0.0S 

500  50000 

flaw 

417 

1918/jull  lAww_filier_caps 

1 

flaws  0.01  ufcap  weaves 

HW2 

ring! 

good 

o.os  o.os 

50050000 

flaw 

monos  ensor 

419 

19tS/jull  lA>o_filier_c«p« 

1 

flaws  0.1  ufcap  monoacn* 

HW2 

ringl 

good 

0.05  0.05 

500  50000 

flaw 

•or  weaves 

421 

1911/jull  5/umpleb  Jiw2 

b 

flaws  0.1  ufcap  tows 

HW2 

ringl 

fair 

0.05  0.05 

50050000 

flaw 

monmmaor 

423 

I9SI/JulI  9/Mmplcb_new_flsw 

b 

flaws  0.1  ufcap  tows 

HW2 

ringl 

good 

aos  aos 

50050000 

flaw 

monoaanaor 

42$ 

1911/jull  9/t«mpleb_new_fl*w_bsck 

b 

flaws  0.1  ufcap  tows 

HW2 

ringl 

good 

0.05  0.05 

50050000 

flaw 

monoaansor 

427 

19»l/jul20Mup_ova_itmpleb 

b 

flaws  0.1  ufcap  tows 

HW2 

arhip 

fair 

0.05  0.03 

50050000 

flaw 

moammtor 

429 

194lflul2flYwhip_over_b*cl 

b 

flaws  O.lufcap  tows 

HW2 

whip 

fair 

aos  o.os 

50050000 

flaw 

monosensor 

431 

19>l/jul21/casba_nn2 

carbs 

flaws  O.lufcap  weaves 

HW2 

HWg 

fair 

0.05  0.05 

S0050000 

flaw 

monosaiaor 

433 

19St/jul2Qfcarba_nml 

cuba 

moacMHor  .lufcap 

HW2 

ringl 

fair 

0.0$  0.0$ 

50050000 

flaw 

43$ 

19U/jul21/c*rtx_nn3 

euba 

moacaauor  .lvfcap 
balazs 

HW2 

whip2 

poor 

0.0$  0.0$ 

50050000 

flaw 

437 

19*l/juI22/carba_i»hip2 

caiba 

man  osmsor  balazs 

HW2 

whip2 

poor 

0.050.05 

50050000 

flaw 

439 

19it/jul22/caib<  cap9 

Cuba 

monosmsor  balazi 

ctp9 

cap9 

poor 

0.0$  aos 

50050000 

flaw 

441 

1 9SS/jul22/cub*_Ioe 

cuba 

monosmsor  balm 

HW3 

ring22 

fair 

0.05  0.05 

100050000 

weave 

443 

1 944/jul22/cubb_tee 

cubb 

monasaisor  balazs 

HW3 

ring22 

fair 

0.05  0.05 

100050000 

flaw 

445 

19U4ul2S/bubb_«nra* 

cubb 

weaves  O.lufcap 
mnwaaisos  balus 

HW2 

ringl 

fair 

0.0$  0.05 

50050000 

flaw 

447 

19l*/jul28/esrbb_foJd 

cubb 

flaws  O.lufcap  meuatn- 
acr  balazs 

disk 

disk 

poor 

0.0$  0.0$ 

505000 

flaw 

449 

HWj»B4te.qBid4 

ua 

flaws  0.1  ufcap  mcooaai- 
aos  balazs  tows 

HW2 

ring] 

fair 

0.1  0.1 

50050000 

flaw 

431 

19tt£iD6/kss_dlam5 

XU 

flaws  0. 1  ufcap  moooaer- 
•oi  balazs  lews 

HW2 

ringl 

fair 

0.1  0.1 

50050000 

flaw 

4S3 

19Mfril2S/caiba_iveaves 

cuba 

monoaanaor  0.  lufcap 
balazs  waavas 

HW2 

ringl 

fair 

.050.050 

50050000 

weave 

V-34 


Underj 

Sample 

Faalara* 

Sana  or 

Scum 

■Uttgc 

Teat 

455 

195SftiB2*otu_iM 

*ol*i9 

_ .,^r,r 

HW3 

ring22 

good 

.025.025 

1000  50000 

crack 

cock*  BamaolarccU 

457 

!9S56uC3Mtr 

*olu5 

manoi  0  1  ufcjp 

HW2 

nag! 

good 

.025  .025 

100050000 

flaw 

cock*  Haw*  aolucell 

459 

19950ul23AeUr2 

•olu5 

Saw  mat  nun*  or 

0.1ufcap  (oluccU 

disk 

fair 

.025  .025 

50050000 

crack 

461 

19SS^ul26boiu3 

aolar7 

halar*  moooaanaor 

0. 1  ufcap  flaw*  crack* 
aolucell 

disk 

diak 

fair 

.025.025 

50050000 

crack 

463 

inVjaWMdii 

kUj5 

cndu  flawt  0.  lufcap 

HW1S 

nng23 

good 

.025.025 

1000  50000 

crack 

465 

lMS/juOT/koUy 

•oUr5 

flaw  .lufcap  moooaensor 
aolaiceU 

HW15 

ring23 

good 

.025.025 

1000  50000 

flaw 

467 

195S^ul27/xu 

Ul 

moooaanaor  tow*  .lufcap 

HW1S 

ring23 

fair 

.099.099 

1000  50000 

flaw 

469 

19S«frtl29/cub9 

cub* 

awnoatnaor  .lufcap 

HW18 

ring23 

fair 

0.100  0.100 

50050000 

weave 

wcavaa 

471 

19SS/jul29/carba 

cub* 

monoautaor  .lufcap 
waavca  balais 

HW19 

ring24 

fair 

0.1000.100 

50050000 

weave 

473 

I9SI/aug01/earba_dct*iled_weave* 

cub* 

weave*  .lufcap 

HW19 

bag24 

good 

.050.050 

50050000 

weave 

moooaanaor 

475 

19SS/aug01fcaibb_giDove_groove«idc 

cubb 

monoaanaar  .lufcap 

HW19 

rtng24 

fair 

.050 .050 

500  50000 

flaw 

477 

l9SS/nig01At>lar_i] 

*olu5 

Bam  .lufcap  monaacn- 
•oraolarccll 

HW19 

nng24 

good 

.025.025 

50050000 

flaw 

479 

19SS/bugOSfeiUr_* 

*olu5 

Sam  .lufcap  moaoaan- 
aor  aolucell 

HW19 

ring24 

good 

.025.025 

50050000 

flaw 

451 

19SS/aug08/cc_fiber joint 

cubb 

mnctoannam  .lufcap 

Sam 

HW19 

nng24 

fair 

.050 .050 

1000  50000 

flaw 

453 

19St/aug09/be_on_iop 

b 

Sam  moooaenaor 

HW19 

ring24 

fair 

.050.050 

50050000 

flaw 

455 

19SI/*ug09Aiecn_boticm 

b 

Sam  tom  tncnoacnacr 

HW19 

nr  *24 

fair 

.050 .050 

50050000 

flaw 

457 

1955/augl  0/B_top 

b 

Sam  tom  moooaenaor 

HW8 

ship 

good 

.050 .050 

50050000 

flaw 

459 

19SSAuglO/B_boaom 

b 

Sam  tom  BMBocaoaor 

HW8 

whip 

fair 

.050.050 

50050000 

flaw 

491 

1955/iugl 1  ffl.tep  jong 

b 

Sam  tom  moooaenaor 

HWS 

ringS 

good 

.050.050 

50050000 

flaw 

493 

1955/iugll/B  bottom  ring 

b 

Bam  tom  monoseoaor 

HW5 

ringS 

good 

.050.050 

50050000 

flaw 

495 

!946^wl5^cbotrdl 

tracts  moooaanaor 

0.  lufcap 

HW2 

ring2 

.010.010 

50050000 

flaw 

497 

19SS/apr26fruelb 

ssel 

moots  manrO-lufcap 

macOl 

whip 

0.1000.100 

500  50000 

499 

19S5/*pr27Aie3 

s«e3 

moooaanaor  0.  lufcap 

jtbl 

ring21 

0.1000.100 

500  50000 

501 

19SSApr27A»e3* 

2M 3 

nutnoautaarO.  lufcap 

JM 

ring21 

0.1000.100 

50050000 

503 

1955/junOlAwlO 

1 

monoacaaorO.  lufcap 

HW10 

ring7 

.099.099 

50050000 

505 

1955/julOlfcarrgtlel  jtew.f.ctp* 

1 

moooaanaor  0.  lufcap 

HW2 

ring2 

.099.099 

50050000 

507 

1955/juM7«l»m4 

Uf 

moatoaanaor  0.  lufcap 

macOl 

lccpS 

0.100  0.100 

S0050000 

aubtracboo  Iraiaaninod 

509 

19SS/jul07Ml*m5 

Ul 

moooaanaor  Ot  lufcap 
aubtracboo  tranaraisod 

macOl 

loops 

0.1000.100 

50050000 

511 

19IS/juXr7Aju*d] 

Ul 

moots  to  aor  0.1ufcap 
aubtracboo  tranunined 

macOl 

loops 

0.1000.100 

50050000 

513 

1955/juXT7A|u*<12 

US 

menoaanacr  O.lufcap 

macOl 

loops 

0.1000.100 

50050000 

aubtracboo  tranamioad 
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Range 


Test 


SIS 

19M0olO7/l9Uid3 

UI 

subtraction  transmitted 

macOl 

loopS 

0.100  0.100 

50050000 

517 

19tVjul07A]ued4 

XU 

manoMnier  0.  lufeap 
subsection  transmitted 

mscO! 

loops 

0.1000.100 

50050000 

619 

!9Si^ull  lfoo_uaalum_cape 

1 

mannstsisnr  0.  lufeap 

HW2 

ringl 

.060.050 

50050000 

521 

19tt^ull3Mlam4 

XU 

imae not  0.  lufeap 
aubtraction  transmitted 

macOl 

loops 

.050.050 

500  50000 

S33 

19M0uU3Mlam5 

XU 

mcnoscnjor  0.  lufeap 
subtraction  uanmuned 

macOl 

loopS 

.050.060 

50050000 

523 

19tt0ull3^uadl 

XU 

subtraction  nararnined 

macOl 

loops 

.060.050 

50050000 

527 

19*l^ull3*ju»d2 

XU 

monosmaor  0.  lufeap 
subtraction  transmitted 

macOl 

loopS 

.050 .050 

S0050000 

529 

19U/jull3A)uad3 

XU 

tnanatcnsorO.  lufeap 
subtraction  transmuted 

macOl 

loopS 

.050.050 

50050000 

S31 

19SS^ull3/quad4 

XU 

monosaiaor  0.  lufeap 
subtraction  transmitted 

macOl 

loopS 

.050.050 

50050000 

533 

19*Vju 115M5nol 

XU 

monosensor  0.  lufeap 
subtraction  transmitted 

macOl 

loopS 

0.1000.100 

50050000 

535 

19M/jull5M5no2 

XU 

monaststsnr  0.  lufeap 
subtraction  transmitted 

macOl 

loops 

0.100  0.100 

50050000 

537 

19tVjalltAUsn>4 

XU 

monoamaor  0.  lufeap 
subtraction  transmitted 

macOl 

loops 

0.1000.100 

50050000 

539 

19*VjullS/dlam5 

XU 

monoamaor  0.  lufeap 
subtraction  traaaninod 

macOl 

loopS 

0.1000.100 

50050000 

541 

19Sl^ul2«i^uadl 

XU 

subtraction  tranamitted 

macOl 

loopS 

0.1000.100 

50050000 

543 

1911/jull  Mquad2 

XU 

subtraction  transmitted 

macOl 

loopS 

0.1000.100 

50050000 

545 

1911/jull  tA]usd3 

XU 

tnonoaansor  O.lufcap 
subtraction  tranamitted 

macOl 

loopS 

0.1000.100 

50050000 

547 

19ttflulll*|uad4 

XU 

aubtraction  trams  timed 

macOl 

loops 

0.1000.100 

50050000 

549 

19tl/jul27/dil!/dUm2 

XU 

monoacnaorO.  lufeap 
subtiactioo  trsnamised 

HW1I 

loopS 

0.1000.100 

100050000 

551 

198t/jul27/diS/<Qam« 

XU 

menoacnaor  0.  lufeap 
subtraction  tranamitted 

HW1* 

loopS 

0.1000.100 

100050000 

553 

19tl/jul27/difl>quidl 

XU 

monosenaor  0.  lufeap 
subtraction  tranamitted 

HW18 

loops 

0.100  0.100 

100050000 

555 

199t^ul27  Miff /qua  d2 

XU 

monoamaor  0.  lufeap 
aubtraction  uanninod 

HW18 

loopS 

0.1000.100 

100050000 

557 

19tt/)ul27/difr/quad3 

XU 

monossnsor  0.  lufeap 
subtraction  tranamitted 

HW1S 

loop5 

0.100  0.100 

100050000 

559 

XU 

monoamaor  0.  lufeap 
subtraction  gammed 

HWlg 

loop5 

0.1000.100 

100050000 

S«1 

19tl/jul21/dlam2 

XU 

monosensor  0.  lufeap 
subtraction  transmitted 

HW4 

loops 

0.1000.100 

500  50000 

563 

191|/jul21/dlamt 

XU 

monoacnaorO.  lufeap 
subtraction  transmitted 

HW4 

loopS 

0.1000.100 

50050000 

V-36 


Mrodarj 

tanple 

Ft* taros 

Saneor 

torn 

Range 

565 

l9ttfril2SAquadl 

Ul 

BOMMUor  0. 1  ufcap 
•uhtrarbnn  transmitted 

HW4 

loop! 

0.100  0100 

500  50000 

567 

19ttfiuOtt|uad2 

lit 

Bncun tot  0. 1  ufcap 
•ubtractioa  tranarained 

HW4 

loupS 

0.1000.100 

50050000 

569 

196^ndAp«d3 

XU 

•uboaoiau  tranamined 

HW4 

loopS 

0.1000.100 

500  50000 

S71 

19U/jul2IA{u«14 

XU 

subtraction  transmitted 

HW4 

loop! 

0.1000.100 

500  50000 

573 

XU 

subtraction  transmitted 

HW4 

loopS 

0.1000.100 

50050000 

575 

19ttfcug0tMiffMlamS 

XU 

subtraction  transmitted 

HW4 

loopS 

0.1000.100 

soosoooo 

577 

19f  1/aufOtAdiffAquid  1 

XU 

moocaensor  0.1  ufcap 
subtraction  tranarninrd 

HW4 

loopS 

0.1000.100 

500  50000 

579 

191tAsuj0t/difl>qusd2 

XU 

■MioMuor  0. 1  ufcap 
subtraction  tranamined 

HW4 

loopS 

0.1000.100 

soosoooo 

511 

1988/su|08/<bfl/quad3 

XAS 

•ubtnction  transmitted 

HW4 

loopS 

0.1000.100 

soosoooo 

513 

198 1/augOl/difT Aqua  d4 

XU 

iwminor  0.1  ufcap 
•ubtractian  tranamined 

HW4 

loopS 

0.1000.100 

50050000 

SIS 

1 98t/auf09/diffAllam2 

XU 

monaaaiaor  O.lufcap 
subtraction  transmitted 

HW8 

loopS 

0.100  0.100 

soosoooo 

517 

I9*l/»uf09/*ff/dlim8 

XU 

■ubtractian  transmitted 

HW8 

loops 

0.100  0.100 

soosoooo 

519 

1988/aug09AdiJIAjuadl 

XU 

monoMnaor  O.lufcap 
subtraction  transmitted 

HW8 

loopS 

0.1000.100 

soosoooo 

591 

19tl/>uf09MilTA]uad2 

XU 

monosouor  0. 1  ufcap 
■ubtractian  transmitted 

HW8 

loops 

0.1000.100 

soosoooo 

593 

1988/auf09/difTAquad3 

XU 

mannaananr  O.lufcap 
■ubtractian  tnnaniaed 

HW8 

loopS 

0.100  0.100 

soosoooo 

595 

19tt/auf09/difl)lquad4 

XU 

•ubtractian  tmuaninad 

HWt 

loopS 

0.1000.100 

50050000 

597 

19tt/auglOMifiMlara2 

XU 

manoMnaor  0. 1  ufcap 
•ubtractian  trananined 

HW8 

loopS 

0.1000.100 

50050000 

599 

IMtfcqglQMifl/dlamS 

XU 

tnonoacnaorOLlufcap 
•ubtraction  transmitted 

HW8 

loopS 

0.1000.100 

soosoooo 

601 

19flfcuflQAdifi>]usdl 

XU 

mnmmnanr  0. 1  ufcap 
•ubtractian  tranamined 

HWt 

loop 5 

0.1000.100 

50050000 

603 

1 988/tugl  OAdifTAquad  2 

lift 

•ubtractian  tranamined 

HWt 

loops 

0.1000.100 

soosoooo 

605 

1988fau»10AiifT/qu4d3 

XU 

monaaaiaor  O.lufctp 
•ubtractian  traasninod 

HWt 

loopS 

0.1000.100 

50050000 

607 

198l/aulI<WbfrAqu.d4 

lift 

•ubtractian  tranamined 

HWt 

loopS 

0.1000.100 

50050000 

609 

19f  l/au|l  l/diffMlam} 

HI 

mooouruor  0.1  ufcap 
•ubtractian  tranamined 

HWt 

loopS 

0.1000.100 

soosoooo 

611 

1988/aujl  lAKfCMlamt 

lift 

monaransor  0.1  ufcap 
■ubtractian  tranamined 

HWt 

loops 

0.1000.100 

soosoooo 

V-37 


Source 


butt 


(13 

19MMifl  1/difDqutdl 

(15 

19M/tujll/difIAq(utd2 

(17 

19(t*u(l  l/difl>qutd3 

(19 

19MMifllMifIA)utd4 

(21 

19Mbugl2M4oi( 

(23 

ifWhuglMSny 

(25 

(27 

19tl/tuI12>q2m( 

(29 

19tt/nj|12ik]3m| 

(31 

19tt/tu|12A)4mg 

(33 

19ll/iu*20/iolii 

(35 

19»Vtu*22/tu 

(37 

19U/tu|23/xu2 

(39 

19((/tug30/dltm 

(41 

19I*/tu*3 1  /xtsmo 

(43 

19IIAep01/iu02 

(45 

19U/upli/xu07 

(47 

I9tt/upl2/xu08 

(49 

19U/upl3fxu09 

(51 

mi^tpivtuio 

(S3 

196Utp26Mltm2 

(55 

19MAtp26/dltm5 

(57 

19Ufep2touqI 

(59 

19IUtp26/xuq2 

(61 

19MAop26/xttq3 

(63 

19U/bcp2*xuq4 

(65 

19U/wp2I/dltm2 

(67 

19tMtp2*MkmS 

in 

XU 

1U 

IK 

Ul 

Ui 

Ul 


•olu5 


dUm 

xu 

XU 

xu 

xu 

xu 

xu 

xu 

xu 

xu 

xu 

xu 

xu 


utxufUon  trxnninitd 


r  Olufctp 
■thtmtior  inmioad 

ubtracuon  sxmuuad 

atbtncben  trottmiaed 


r  0.1ufctp 
ubutcbor  uxtuniiMd 


r  Olufctp 
•uhmction  tnn&niucd 


rO.lufetp 
tubtnction  mumined 


i  Olufctp 
tubmctjoc  paunuued 


r0.1ufctp 
tufamction  transmiucd 


r  O.lufcip 
tubcrtcuar  trusiuned 


rO.lufetp 
tubtrtcuco  tmnninnd 


rO.lufetp 
tubtnction  tmwniued 

tubtnction  tnnsnined 


rO.lufetp 
•ubtrxctkn  titntminod 


rO.lufetp 
tubtncoaa  tnmtniuod 

tubtnetiee  tztntntiaed 

■ibtnction  tnnmiued 


rO.lufetp 
tubtnction  tnrwnincd 


HWS 

leopS 

0.100  0.100 

HWI 

loopS 

0.100  0.100 

HWI 

loop5 

0.1000.100 

HWI 

loopS 

0.1000.100 

HWI  9 

rin(24 

0.1000.100 

HWI  9 

an(24 

0.100  0.100 

HW19 

hng24 

0.100  0.100 

HW19 

rin(24 

0.1000.100 

HWI  9 

ring24 

0.1000.100 

HW19 

ringjd 

0.100  0.100 

HW20 

ring25 

.010.010 

HW20 

Dag25 

.019  .020 

HW20 

hn|25 

.019 .020 

HW20 

ting25 

.050 .050 

HW20 

zm$2S 

.041  .040 

HW20 

ring 25 

.040  040 

tut 

.0(2.012 

ton 
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5  b.  Sample  1  Data 


Laboratory  data  considered  in  Chapter  III  was  taken  from  Sample  1,  a  satin  weave 
sample  with  twelve  flat-bottom  drilled  holes  of  two  different  diameters  and  six  different 
depths.  A  complete  discussion  of  the  modeling  and  inversion  of  these  data  appears  in 
Chapter  HI.  Here,  images  from  the  Tons-o-Data  booklet  are  presented  before  any  rotat¬ 
ing  or  scaling,  as  outlined  above,  is  done  to  the  data. 
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Figure  V-5.  Drawing  of  Sample  1,  satu  weave  graphite  epoxy  panel  with 
twelve  flat-bottom  holes  of  two  diameters:  0.5"  and  0.25".  Depths  of  the 
holes  are  as  indicated. 


The  layout  of  Sample  1  appears  in  Figure  V-5.  Depths  of  the  holes  are  shown  in  inches. 
The  overall  size  of  the  sample  is  6"  by  6".  The  larger  holes  are  0.5"  diameter;  the  smaller 
holes  are  0.25"  diameter.  All  holes  in  this  sample  are  flat-bottom  machined.  Pages  from 
data  profiles  are  presented  in  Figure  V-6. 


5  c  Sample  3  Data 

Sample  3  is  a  satin  weave  sample  of  graphite  with  one  0.25"  diameter  flaw  in  the 
center.  Thus,  it  is  a  convenient  comparison  with  model  calculations:  it  is  easy  to  describe 
the  flaw  to  the  model,  and  relatively  easy  to  interpret  the  results  since  they  are  known  to 
come  from  a  single  flaw  rather  than  a  number  of  flaws.  Some  tests  were  done  with  Sam¬ 
ple  J,  both  right-side-up  and  up-side-down.  The  topside  images  are  presented  here 
because  they  represent  a  convenient  measurement  to  be  considered  for  model  calcula¬ 
tions. 
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Figure  V-6.  Data  from  Tons-o-Data  booklet  used  in  the  model  for  recon¬ 
struction  of  flaws.  See  previous  figure  for  layout  of  actual  flaws.  Data  here 
are  “raw”  laboratory  data,  without  normalization  and  phase  angle  rotation. 
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Figure  V-7.  Data  from  Sample  3 ,  a  satin  weave  with  a  single  flat-bottom 
machined  hole  in  the  center.  These  data  are  convenient  for  reconstruction 
because  of  the  relative  simplicity  of  the  flaw  geometry.  The  actual  flaw  is 
approximately  0.1"  deep  and  0.25"  in  diameter.  Data  here  are  “raw”  la¬ 
boratory  data,  without  normalization  and  phase  angle  rotation. 
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